DUU 

HAVA 

|40NT>_fil'> f 



jplE SCrtOOW 



NAVAL POSTGRADUATE SCHOOL 

Monterey, California 




THESIS 

THERMAL ANALYSIS AND DESIGN OF AIR COOLED 
ELECTRONIC CIRCUIT BOARDS USING 
A DESKTOP COMPUTER 

by 

Ricky Allen Foltz 
June 1980 

Thesis Advisor: Matthew D. Kelleher 

Approved for public release; distribution unlimited 



T 1 96584 



UNCLASSIFIED 



SECURITY CLASSIFICATION OF THIS page n»,m Bnffd) 



REPORT DOCUMENTATION PACE 


READ INSTRUCTIONS 
BEFORE COMPLETING FORM 


REPORT HUMIC* 


2. GOVT ACCESSION NO. 


». RECIPIENT'S CATALOG NUMBER 


4 TITLE fond SulHttf) 

Thermal Analysis and Design of Air Cooled 
Electronic Circuit Boards Using a Desktop 
Computer 


S. TYRE OF REPORT A PERIOD COVERED 

Engineer's Thesis, 

June 1980 


• ■ PERFORMING ORG. REPORT NUMBER 


7. AuTMOWtj 

Ricky Allen Foltz 


t. CONTRACT OR GRANT NUMBER|«J 


• PERFORMING ORGANIZATION NAME ANQ ADDRESS 

Naval Postgraduate School 
Monterey, California 93940 


10. PROGRAM ELEMENT. PROJECT. TASK 
AREA A WORK UNIT NUMBERS 


II CONTROLLING OFFICE NAME and ADDRESS 

Naval Postgraduate School 
Monterey, California 93940 


12. report DATE 

June 1980 


is. number of pages 

135 pages 


14 MONITORING AGENCY name A AOCRESS/II rflt/aronf tram 1 Controlllni 01(1 co) 

Naval Postgraduate School 
Monterey, California 93940 


«S. SECURITY CLASS, (oi thle report) 

Unclassified 


is«. OECLASSIF* CATION/ DOWNGRADING 
SCHEDULE 


i* Distribution statement (oi thi « 

Approved for public release; distribution unlimited 


17 DISTRIBUTION STATEMENT fat t»m «k*ir*cl entered tn sided 30. It dllloronl teom Report) 


<s supplementary notes 


Electronics cooling, interactive computer program, Thermal Network 
Analysis 


A thermal design procedure for air cooled electronic circuit boards 
has been developed for the Hewlett-Packard Model 9845 desktop com- 
puter. The system of interactive programs, called THERMELEX, 
performs thermal analysis of printed circuit boards to predict 
either the junction temperatures for given power dissipation levels 
or the maximum power levels for given junction temperature limits. 
The system includes the following features: totally interactive 

with all input in question and answer format, simple data 



DD , jj"“i 1473 eoit.onommovmiioiioliti UNCLASSIFIED 

(Page 1) S/N 0102-014. AA01 



1 



SECURITY CLASSIFICATION OF THIS FAOI Omi» Bmorod) 



UNCLASSIFIED 

t«eu*«Tv cl *»» i>tc atiq- ow Twit n»*m t»« ( 



verification and correction capabilities, ability to store and 
retrieve circuit board descriptive data totally under program 
control, wide variety of output formats including tabular and 
graphical. By using internal selection of heat transfer corre- 
lations, the THERMELEX system depends only on input of physical 
parameters for thermal predictions. 



DD Forr^ 1473 
S/N oV&2 -014-6601 



2 UNCLASSIFIED 



lieuUlTV CWAMtHCATlON 0* ▼Nlf PAflt Dmtm 



Approved for public release; distribution unlimited. 



Thermal Analysis and Design of Air Cooled Electronic 
Circuit Boards Using a Desk Top Computer 



by 

Ricky Allen ^oltz 

Lieutenant Commander, United States Navy 
B . S . , Stanford University, 1971 



Submitted in partial fulfillment of the 
requirements for the degree of 



MASTER OF SCIENCE IN MECHANICAL ENGINEERING 

and 

MECHANICAL ENGINEER 



from the 

NAVAL POSTGRADUATE SCHOOL 
June 1980 



■"vn IOtt;T 



Hjsm 



ABSTRACT 

A thermal design procedure for air cooled electronic 
circuit boards has been developed for the Hewlett-Packard 
ilodel 9845 desktop computer. The system of interactive 
programs, called THERMELEX, performs thermal analysis of 
printed circuit boards to predict either the junction tem- 
peratures for given power dissipation levels or the maximum 
power levels for given junction temperature limits. The 
system includes the following features: totally interactive 

with all input in question and answer format, simple data 
verification and correction capabilities, ability to store 
and retrieve circuit board descriptive data totally under 
program control, wide variety of output formats including 
tabular and graphical. By using internal selection of heat 
transfer correlations, the THERMELEX system depends only on 
input of physical parameters for thermal predictions. 
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NONMENCLATURE 



Aair Average area for cooling air flow 

Agap Area of bottom of component 

Rl Area for conduction in either a CPU 

or a circuit board lead 

Apara Area of the component experiencing 
parallel flow 

Areg Area of the region on the circuit 
board 

Aspin Surface area of the component lead 

Astag Area of the component experiencing 
stagnation flow 

Axpin Crossectional area of the component 
lead 

Bdh Height of circuit board (perpen- 
dicular to air flow) 

Bdl Length of circuit board (parallel 

to air flow) 

Cpair Heat capacity of the air 

CPU Conduction Path Unit (size 

defined by user) 

Dgap Distance between the component 
bottom and the circuit board 

Dh Hydraulic diameter 

DIP Dual Inline Package 

Epsb Emissivity of the circuit board 

Epse Emissivity of the component 

Fair Volumetric air flow rate 



[m 2 ] 

[m 2 ] 

[m 2 ] 

[m 2 ] 

[m 2 ] 

[m 2 ] 

[m 2 ] 

[m 2 ] 

[m 2 ] 

[m 2 ] 

[J/Kg-degK] 



[m 2 ] 

[ra 2 ] 



[m 2 /sec] 
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Fr 


Roughness factor 




Gnu 


Kinematic viscosity of air 


[m 2 /sec] 


Gz 


Graetz number 




Havg 


Average height of the components 
present on the board 


[m] 


Hb 


Heat transfer coefficient for the 
circuit board 


[Watt/m 2 -degK] 


Hduct 


Heat transfer coefficient for the 
duct formed by the circuit boards 
and equipment case 


[Watt/m 2 -degK] 


He 


Height of the component 


[m ] 


Hpara 


Heat transfer coefficient for areas 
receiving parallel air flow 


[Watt/m 2 -degK] 


Hr 


Dimension of region in the vertical 
direction or vertical distance between 
finite difference nodes [m ] 


Hstag 


Heat transfer coefficient for areas 
receiving stagnation air flow 


[Watt/m 2 -degK] 


Kair 


Thermal conductivity of cooling air 


[Watt/m-degK] 


Kb 


Thermal conductivity of the circuit 
board 


[Watt/m-degK] 


Kl 


Thermal conductivity of the board 
leads 


[Watt/m-degK] 


Kpin 


Thermal conductivity of the compo- 
nent lead 


[Watt/m-degK] 


Le 


Length of the component 


[m ] 


Lpin 


Length of component lead 


[m ] 


Lr 


Dimension of region in the hori- 
zontal direction or horizontal 
distance between nodes 


[m ] 


Mair 


Mass flow rate of cooling air 


[Kg/sec] 


Navg 


Average number of component in a 
vertical column perpendicular to 
the air flow direction 
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Npin 

Nxr 

Nyr 

Perim 

Pow 

Pr 

Q 

Oe-b 

Q i 

°2 

Q 3 

R 

Re 

Rb- 

conv 

Rb-hor 

Rb-ver 

Re- 

conv 



Number of leads on a component 

Number of regions in the air flow 
of horizontal direction 

Number of regions perpendicular to 
the air flow direction or vertical 
direction 

Wetted perimeter of the air duct 

Energy dissipated in the component 

Prandtl number 

Rate of heat transfer 

Total rate of heat transfer from 
the component to the circuit board 

Total rate of heat transfer out the 
bottom of a circuit board region 

Total rate of heat transfer out the 
right side of a circuit board region 

Total rate of heat transfer out the 
top of a circuit board region 

Total rate of heat transfer out the 
left side of a circuit board region 

Thermal resistance for heat 
transfer 

Reynolds number 

Total thermal resistance for convec- 
tion from circuit board surface to 
the cooling air flow 

Thermal resistance for conduction 
in the circuit board between nodes 
in the horizontal direction 

Thermal resistance for conduction 
in the circuit board between nodes 
in the vertical direction 

Total thermal resistance for con- 
vection from component surface to 
the cooling air flow 



[m ] 

[Watts] 

[Watts] 

[Watts] 

[Watts] 

[Watts] 

[Watts] 

[Watts] 

[deg K/Watt] 

[degK/Watt] 

[degK/Watt] 

[degK/Watt] 

[degK/Watt] 
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[degK/Watt] 



Rgap- Conductive thermal resistance for 
cond the air gap between the component 
bottom and the circuit board 

Rgap- Effective thermal resistance for 
rad radiation from the component bottom 
to the circuit board 

Rl Total conductive resistance between 

finite difference nodes (centers of 
region) 

Rpara Thermal resistance for heat transfer 
from area receiving parallel air 
flow 

Rstag Thermal resistance for heat transfer 
from area receiving stagnation air 
flow 

Rtop-r Effective thermal resistance for 
radiation from the top of the 
component to next board 

Rtote- Total thermal resistance between 
b the component bottom and the 

circuit board 

Tair Local air temperature 

Te Temperature of the component surface 

Tb Temperature of the board surface 

Thb Thickness of the circuit board 

Tj Temperature of the component 

junction 

Vair Velocity of the cooling air 

Wavg Average width of the components 
present on the board 

We Width of the individual components 

X Direction parallel to cooling air 

flow also referred to as the hori- 
zontal direction due to orientation 
of circuit board picture on screen 
of computer 



[degK/Watt] 

[degK/Watt] 

[degK/Watt] 

[degK/Watt] 

[degK/Watt] 

[degK/Watt] 
[degC, degK] 
[degC, degK] 
[degC, degK] 
[m ] 

[degC , degK] 
[m/sec] 

[m] 

[m ] 
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Xi Distance from the entrance of the [m] 

cooling air 

Y Direction perpendicular to cooling 

air flow also referred to as the 
vertical direction due to orienta- 
tion of circuit board picture on 
screen of computer 

Zb Distance between circuit boards [m , 



mm ] 
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I . INTRODUCTION 



A. BACKGROUND 

Electronic components generally convert a significant 
fraction of the input power into internal joulian heating. 

When the components are large, widely separated, and air is 
permitted to circulate freely among them, this heat generated 
within the component is transferred to the environment through 
natural convection. Indeed, millions of pieces of home elec- 
tronic equipment have operated reliably for years depending 
only on natural circulation for cooling. The major emphasis 
over the last 15 years has been a continuing effort toward 
a reduction in physical size of components and increased 
component density within electronic equipment. The military 
is particularly in need of smaller, more complex, yet relia- 
ble equipment that must often be sealed from an extremely 
hostile environment. This results in the need for more 
efficient cooling methods. 

The age of large scale integration (LSI) is here. 

Although the power dissipated in each active junction within 
a component has greatly decreased, the large number of heat 
sources on each chip has created serious heat removal 
problems. There is a genuine need for electronic designers 
to not only be aware of the problems associated with higher 
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temperatures, but they must also be prepared to solve those 
problems. References [1] through [10] are a sampling of 
the many sources that indicate the increased emphasis on 
cooling problems throughout the electronics community. 

The need to operate electronic equipment with maximum 
junction temperatures below those levels that result in 
failure is well recognized, but even moderately high temper- 
atures in electronic components result in progressive 
deterioration and reduced reliability. A generally accepted 
thumb rule is: for every 10 deg C increase in junction 

temperature, the lifetime of a component will decrease by 
one-half. [3, 7] The increasing costs associated with each 
failure make it imperative to address reliability during 
the equipment design phase and to provide sufficient cooling 
to maintain temperatures as low as practical. 

There are many methods of removing the heat dissipated 
within electronic components. These methods include the 
natural air convection previously mentioned, thermo-electric 
devices, heat pipes, cold plates and even complex refrigera- 
tion systems which use pumped liquid coolants [8]. The 
complex systems required for these solutions can lead 
to lower reliability through increased probability of 
failures in the cooling system. While many of these methods 
can allow extremely high power densities for specific appli- 
cations, the most widely used method for cooling of components 
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mounted on printed circuit boards is forced air cooling. 

Air is readily available, abundant, non-corrosive, non-toxic, 
non-flammable, dielectric, and is easily pumped using readily 
available fans and blowers. For those situations where there 
is a need to seal the equipment, interior cooling air is often 
circulated through heat exchangers. 

All too often the problem of cooling the electronic com- 
ponents is attacked after the circuit boards have been designed 
or even produced and assembled [9], In some cases, the need 
for increased cooling is recognized only after repeated 
failures in service have resulted in unhappy users and ex- 
cessive repair costs. At best both the electronic design and 
the thermal design progress concurrently but are worked on by 
separate design groups. These groups may have conflicting 
range goals that can interfere with the important long range 
goal of reliability. 

Many of the efforts to directly involve the electrical 
designer in the thermal problems are directed towards overall 
system cooling. Reference [11 ] details a thermal design 
program called VENTBOX. VENTBOX treats a cabinet enclosure 
with entire circuit boards modeled as distributed heat 
sources. This program fails to address the problems of in- 
dividual components and, therefore, is unable to analyze 
the circuit board. 

Electronic circuit analysis programs are often utilized 
in the thermal analysis of circuit boards [12 and 13] . This 



15 



technique requires the development of an equivalent thermal 
circuit, usually by a packaging engineer, and thus the thermal 
design is removed from the hands of the electronic designer. 
Reference [14] details a thermal analysis program for circuit 
boards that is much more accessible to an electronic designer 
since the inputs are physical dimensions and types of compo- 
nents rather than equivalent thermal network parameters. 

All these programs share a major weakness in that they 
depend on a large general purpose computer facility. This 
not only can result in excessively long turn-around times 
in batch processing, but also the expense of computer time 
may become a factor. In addition, the input data typically 
consists of long lists of numbers that must be in the proper 
format with the correct option selection codes. Likewise 
the output also consists of even longer lists of numbers 
with the key information hidden in their midst. 

One solution to automated thermal design and analysis 
of avionics systems is being developed as a joint project 
by the Air Force Flight Dynamics Laboratory at Wright 
Patterson Air Force Base and Boeing Aerospace Company. The 
ITAD (Integrated Thermal Avionics Design) system is expected 
to include many existing analysis programs and a large ever 
growing data base containing reliability parameters. It 
is anticipated that ITAD will be accessed by user through 
both batch methods and interactively through color graphics 
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terminals. The scope of this project is enormous and the 
benefits to the system designer will be many, provided he 
has access to the large computer at Wright Patterson Air 
Force Base [15, 16]. 

Another solution on a much smaller scale is to create 
a program for a compact, stand alone desktop computer that 
is easily accessible to the designer of electronic circuit 
boards. Such a system, if interactive and "friendly" to 
the casual user, would allow rapid evaluation of various 
circuit board designs at the conceptual level. Rejection 
of those designs with poor thermal characteristics could 
then occur before they leave the drafting table saving both 
future efforts and dollars. 

B. OBJECTIVES 

The main objective of this thesis was to develop an 
interactive thermal analysis program utilizing the Hewlett- 
Packard 9845 desktop computer. It was considered important 
to include the following features: 

1. Simplicity of operation: All data input and control 

of program flow are done in a question and answer format in- 
cluding specific instructions with each question. 

2. Graphical data checking: Display circuit board 

replicas on the screen both for verification of input and 
to provide a more useful form of output. 
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3. Data correction capabilities: Use interactive methods 

to allow correction of portions of the data without the need 
to repeat all the input. 

4. Automatic storage and retrival of data: Circuit board 

descriptions should be written to and read from mass storage 
devices using simple questions and answers rather than requir- 
ing specific knowledge of operations of the devices. 

5. Analysis methods hidden: Various empirical heat 

transfer correlations should be used based on the physical 
descriptions and the user should be relieved of the need to 
make decisions concerning the details of the heat transfer 
analysis methods. 

6. Sensitivity analysis: Allow automatic parameter 

changes to investigate the effects on the thermal performance 
of these changes with plotted data. 
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II. DESCRIPTION OF THE THERMELEX SYSTEM 



A. GENERAL 

1 . The Computer 

The name THERMELEX will be used to refer to the 
system of interactive programs created to perform thermal 
analysis of air cooled electronic circuit boards. The 
programs are written for the Hewlett-Packard Model 9845 
desktop computer using the Hewlett-Packard extended version 
of the basic language. A sophisticated operating system 
hides the complexity of the 9845 from the user and provides 
protection from his mistakes yet provides the power and 
flexibility needed for application programs such as the 
THERMELEX system. 

The computer used to create and debug THERMELEX is 
the 9845A with 64 k bytes of read/write memory (option 203) , 
internal printer (option 500) and graphics package. In 
addition, dual tape drives and dual floppy disks were utilized 
in the development of the system. These additional mass 
storage devices provided considerable increase in the con- 
venience and speed of storage and retrieval operations but 
THERMELEX is designed to operate with only the standard tape 
drive. Recent advances in the internal electronics have 
resulted in this model being superseded by the 9345B model 
with larger memory capabilities and several other features. 
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THERMELEX is available in either an A-version or B-version 
for use in the respective model of the 9845. 

Figure 1 shows the 9845A but the 9845B is identical 
in physical appearance except for the name tag. The screen 
at the top will display the messages from the system to the 
user, what is typed by the user and the output from the 
THERMELEX programs. This output may be printed on the screen 
in what is known as the alpha mode or it may be presented as 
pictures and words from the graphics mode. Output of either 
mode may also be produced on the thermal printer above the 
keyboard in the inclined area. Directly underneath the 
screen there are four pull-out reference cards that explain 
error messages and other operating conditions that may occur. 
At the extreme upper right corner of the sloping section is 
the standard tape transport (:T15) for mass storage of 
programs and data. 

The keyboard contains not only a set of standard 
typewriter keys but also several other groups of keys that 
are important to THERMELEX. The numeric keys at the lower 
right allow easy entry of numbers as an alternative to 
those in the alphanumeric group of keys and allows numeric 
calculations to be performed even during execution. Between 
these groups are the gold colored program control keys. All 
responses to questions are followed by pressing the CONT key 
at the bottom of this column of keys. The group of special 
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function keys in the upper right are defineable by the user 
or from program control. They are also used in THERMELEX 
as priority interrupts to simplify the input of data. This 
use of these keys is further explained later. The final key 
of interest to the THERMELEX system user is the AUTOST key 
in the lower right corner of the EDIT/SYSTEM FUNCTIONS group. 
This key allows automatic loading and execution of the first 
program in the THERMELEX system. Explanation for its use 
is found in Appendix 1, the User Instructions. For further 
explanations of the features of the Hewlett-Packard model 
9845 computer, see Ref [17] supplied with the computer. 

2. THERMELEX System 

Due to the memory limitations of the 9845A, THERMELEX 
is divided into three major programs which are generally 
brought into case from mass storage under program control as 
they are needed. However, each program is designed to stand 
alone which can allow the experienced user to bypass some of 
the questions and answers required to help the inexperienced 
users. AUTOST is the first program of the three; here several 
pages of user instructions (see Appendix 1) may be printed, 
the special function keys are defined and the user is directed 
along a path to follow towards the other programs. In BOARDS , 
the circuit board descriptive data is entered, verified and 
stored on any one of the mass storage devices. The user may 
also elect to produce a printed copy of the input data for 



21 



his records. The THERML program performs the thermal analy- 



sis and provides output in various tabular or graphical forms. 
These three major programs will be described in more detail 
in the pages that follow. 

There are also several smaller files in the system. 
BDSKEY and STDKEY alter the definitions of the special 
function keys. TNAMES contains a listing of all variable 
names and program section names with explanations to aid in 
any future modifications of the THERMELEX system. DEMO-P 
and DEMO-T are data files containing example circuit board 
descriptive data to allow demonstrations of the system and 
to help the new user become familiar with the capabilities. 

B. AUTOST 

1. General 

AUTOST is the driver program for the system. This 
name allows the program to be automatically loaded and exe- 
cuted if the AUTOST key is latched down and if the tape 
containing the THERMELEX system is in the standard right- 
hand tape drive (:T15) when the main power switch is turned 
on. Regardless of how the program is loaded, the first 
question will be concerning the location of THERMELEX. The 
program will ask which mass storage device contains the 
system since this is necessary for correct program control. 

At this point, the program defines the default mass storage 
device using the "MASS STORAGE IS " " command, and the 
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default mass storage device should not be altered while using 
THERMELEX. 

2 . Instructions 

AUTOST will provide a printed set of user instruc- 
tions either on the screen in short segments or on the thermal 
printer in 8.5 by 11 inch pages. These instructions present 
an overall system view and are intended to familiarize the in- 
experienced user and supplement the extensive instructions 
and prompts that are presented on the screen in all sections 
of THERMELEX. 

3 . Special Function Keys 

AUTOST redefines the special function keys as required 
for the system and can provide a paper key-code overlay as a 
substitute for the plastic model (Hp Part #7120-6164) . Figure 
3 shows both styles. The paper style may be used as a guide 
for filling in the appropriate blanks on the plastic version. 

4 . Option List 

Finally, AUTOST provides a list of options for the 
user to select from. He may load either of the other two 
programs in THERMELEX or produce a complete copy of the 
programs in the system through selection of the appropriate 
number from the menu of options. This technique is employed 
throughout the system whenever possible because of the sim- 
plicity involved in entering a single number over other 
methods of option selection. 
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5 . The Copy Option 

The copy option is included as a built-in feature 
since there are many separate files in THERMELEX and a 
separate command is required for each file to be copied. 

Any mass storage device may serve as the destination for 
the system; however, when copying to tape, it is suggested 
that it be blank due to the number of open blocks required. 
Following completion of this option, program control 
returns to option list. It is expected that the first 
option (keyboard entry of circuit board description) is the 
most likely to be used and the description that follows 
assumes this path. 

C . BOARDS 

1 . General 

The major purpose of the BOARDS program in the 
THERMELEX system is to provide a simple method for input 
of the descriptive data needed in THERML for the creation 
of the thermal model and the analysis. This data includes 
the geometric and material properties of the circuit board 
as well as what components are on the circuit board. These 
components are limited to DIP'S (Dual Inline Packages) and 
chip carriers (FLAT Packs) . The thermal model used for the 
analysis of the board is a finite difference scheme and the 
physical location of the components is approximated by the 
locations of the finite regions created by the user. 
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The descriptive data may be entered from the keyboard 
in response to simple questions; or, it may be retrieved from 
mass storage data files previously recorded using BOARDS . 

The instructions included as Appendix 1 provide specific 
explanations of how to use this program. 

2 . Checking of Input Data 

All input data is displayed using either the graphics 
capabilities of the 9845 or with printed lists produced on 
the screen or paper. The user may, therefore, easily verify 
the validity of the data he has entered. This feature is 
included whenever possible throughout THERMELEX. 

3 . Correction of Input Data 

There are many opportunities to change the descriptive 
data using the same question and answer format of interactive 
programming. At no time is it necessary for the user to have 
extensive knowledge of the program's internal workings or the 
machine level commands performed within the program. As a 
further means of correcting past errors, the backup feature 
built into the THERMELEX system makes it extremely easy to 
return to previous questions. This feature is accessed by 
depressing the special function key (Ko) in place of a 
response to a request for input. Program control will jump 
back to the previous question to allow re-entry of data. 
Repeated use of Ko allows backing up to any desired point 
and resumption of program flow from that point. 
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4 . The Circuit Board 



To describe the geometry of the circuit board under 
investigation, the user must enter the length, width, thickness 
and the thermal conductivity of the circuit board. With no 
input, the thermal conductivity will be assigned a default 
value of 0.29 Watts/m-degC , a representative value for the 
bonded glass laminates generally used for construction of 
circuit boards. The length parameter is defined to be in the 
direction of air flow regardless of which dimension is the 
largest. 

In order to create the finite difference thermal 
model, the board must be divided into logical regions by 
placing any number of equally spaced vertical lines and any 
number of equally spaced horizontal lines on the board up to 
a total of 50 regions. More regions would be possible in the 
9845B due to the larger memory capacity but this would 
require some program changes. Each region will either contain 
a component assumed to be centered in the region or will be 
empty. Since there may be regions with no components, there 
is no unique set of descriptive data for a given circuit 
board design. The user is encouraged to try different com- 
binations of horizontal or x regions and vertical or y 
regions to describe a circuit board. Figure 4 shows two 
such possibilities for a board with six 14 pin DIP'S. Both 
the six region model and the 42 region model describe the 
same board. Experience is helpful in making decisions as 
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to how to divide the board and some boards may not fit into 
the constraints imposed by THERMELEX, but many will. 

5 . The Components 

Each of the defined regions will be empty or contain 
one of the following components: 14, 16, 24, 40 pin DIP'S 

with either vertical or horizontal orientation and 16, 24, 

40, 64 pin chip carriers. These components will be assured 
to be centered in the region (with the exception of the 
40 pin DIP which will be assumed to occupy two regions) . 

The input of component type for each region is via the 
special function keys. A special paper key-code overlay for 
use in this component input section may be provided if a 
plastic version is not available, see Figure 5. As a 
flashing cross appears in each of the regions, the appro- 
priate special function key is depressed. This defines the 
type of component, length, width, orientation and draws the 
component on the screen for visual verification of correc- 
tions. After all regions have been defined, corrections 
are possible through removal and replacement of components 
using the same special function keys. 

6 . Thermal Conduction Paths 

Each of the defined regions on the board may exchange 
energy with the adjacent regions via conduction through the 
board itself and any thermally conductive material added to 
the board such as electrical leads. In addition, many circuit 
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board designs involve metal conduction rails that provide 
both mechanical support and a method to transport the heat 
from the componets. There are provisions in BOARDS to model 
either the electrical leads or the conduction rails. For 
circuit boards without conduction rails, the user specifies 
the average lead width (mm) , thickness and thermal conductivity 
of the lead material. These geometric parameters are used 
to determine the thermal resistance of a single lead connect- 
ing the region centers or nodes. 

When conduction rails are included on the circuit 
board, the effects of the electrical leads are ignored. To 
allow modeling of rails that may have different widths, the 
concept of a conduction path unit (CPU) is introduced. A 
CPU is defined to be . 1 mm in width but the user specifies 
the thickness and conductivity of the material. These para- 
meters are used to determine the thermal resistance of a 
single CPU connecting the region centers. 

With the thermal resistance of single CPU or elec- 
trical lead determined, the user need only specify the number 
of such resistances that connect the circuit board regions. 

For example, if a conduction rail is 5.8 mm in width, it may 
be modeled as containing 58 CPUs. While the concept of 
CPUs is totally artificial, this concept does allow modeling 
of the Navy Standard Electronic Module (SEM) and Improved 
Standard Electronic Module (ISEM) for those situations when 
these modules are used in forced air cooled systems. 
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With the physical description of the circuit board 
now complete, the user may elect to produce a picture of 
the circuit board from a dump of the graphics to the internal 
printer. Pressing special function Key 3 will produce a 
picture such as Figure 6. This picture may be used as a 
final verification of the validity of the circuit board 
description that has been entered. The remaining descrip- 
tive data to be entered concerns the component parameters. 

7 . Temperature or Power Input 

The user must specify either the average power to be 
dissipated in each component or the maximum junction tempera- 
ture. When power levels are specified, the steady state 
junction temperatures will be calculated in THERML . When 
junction temperatures are specified, the maximum allowable 
component power levels will be calculated. In either case, 
the user must also specify a case to junction thermal resis- 
tance for each component. This is usually supplied by the 
manufacturer and provides the means in the thermal model to 
link the component power level and junction temperature to 
the case surface temperature. Figure 7 shows the data 
summary sheet provided after the data has been entered and 
verified . 

8 . Storage of Data 

Although it is possible to load THERML and perform 
the thermal analysis directly since all data is passed 
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through a common block, it is strongly urged that the circuit 
board description be stored on mass storage. A few simple 
answers allows the storage to tape or disk for retrieval at 
a later time. This data file may be accessed, verified, 
altered and recorded back to mass storage using BOARDS to 
investigate the effects that changes in design have on the 
temperatures or power levels obtained in THERML . Both 
Appendix 1 and the program provide easy to follow instructions 
for retrieval of previously recorded descriptive data files. 

Totally under program control, the third major program 
in the THERMELEX system, THERML will be loaded at the direction 
of the user and execution started at the proper entry point. 

D. THERML 

1 . ■ General 

Recall from the previous discussion that BOARDS es- 
tablishes all the parameters that describe the circuit board. 
While the major purpose of THERML is to analyze the thermal 
performance, it must first establish the parameters that 
describe how the circuit board is installed and cooled, such 
as the spacing between boards, volumetric air flow per board 
and inlet temperature of the cooling air. Using the total 
set of board and installation parameters, the thermal model 
is set up, solved and the results presented. 

The approach taken in the thermal analysis is to 
construct an approximate thermal network [7] for every path 
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of heat transfer from the components and the circuit board. 
The thermal resistance (R) for each path is calculated from 
the set of parameters and if a temperature difference (AT) 
for that path is known, the expression: 




will determine the rate of heat transfer (Q) for that path. 

When the rate of heat for a given path is known, the expression 

AT = (Q) x (R) 

will determine the temperature difference. 

An overall heat balance for the components and the 
circuit board, with all the heat that is generated going into 
the cooling air stream, is used to determine either the steady 
state junction temperatures (when component power levels are 
specified) or the component power levels (when the junction 
temperatures are specified) . 

2 . Air Flow - Thermal Model 

The cooling air flow is assumed to come from an in- 
finite heat sink which supplies air at a constant inlet 
temperature. All the heat dissipated within the components 
is assumed to enter the air stream with no heat conducted 
into the card guides or electrical connectors. These guides 
and the mechanical support sections of the connectors are 
typically made of plastic with high resistance to heat 
transfer and this assumption leads to conservative estimates. 
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As the air travels from inlet to outlet, it is assumed 
to remove heat from both the component surfaces and the surface 
of the circuit board. It is further assumed to travel in air 
lanes defined by the height of each region perpendicular to 
air flow and not mix until the outlet. As the air removes 
heat from each region, the local air temperature (Tair) will 
increase. The process is described by the general equation: 



Tair at X + AX = Tair at X + 



Q added in AX 
(Mair) x (Cpair ) 



where Mair = Mass flow rate of air [kg/sec] 

Cpair = Heat capacity of air [J/kg-degk] 

AX = length of a region in the air flow direction 
This process results in the temperature of the cooling air 
stream being modeled as a series of steps as the air travels 
from inlet to outlet within an air flow lane. The local air 
temperature and the convective resistance determine the local 
convective heat transfer. 

The air velocity is determined from the physical 
description entered by the user and this air velocity is 
used to determine the convective resistances. Recall that 
circuit board spacing (Zb) , board height (Hb) and volumetric 
air flow rate (Fair) are part of the descriptive data set 
previously entered. With no components present, the air 
velocity (Vair) would be defined by: 



Vair = 



Fair 

(Zb) x (Hb) 
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However, the components tend to reduce the air flow area 
(Aair) by blockage. The average blockage area is determined 
by calculating the average number of components in a vertical 
column of regions perpendicular to the air flow. In addition, 
the average width (Wavg) and height (Havg) of the components 
is determined and the average air flow area calculated using: 



Aair = (Zb) x (Hb) - (Navg) x (Havg) x (Wavg) 



This area is used to determine the air velocity from: 
Fair 



Vair 



Aair 



In addition the average area is used to determine the wetted 
perimenter (Perim) of the air duct formed by the boards 
according to: 



Perim = 2Hb + 2Zb + 2 (Navg) x (Havg) 



The hydraulic diameter (Dh) is therefore: 
4 Aair 



Dh = 



Perim 



which is also used in calculations of some of the heat 
transfer coefficients. 

3. Components - Thermal Model 

All heat generated within the components must be 
transferred away through conduction to the board, radiation 
to the board and to some radiation sink, and through con- 
vection to the cooling air. Figure 8 shows each of the 
thermal networks for heat transfer. 
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DIP components are usually mounted by soldering the 
electrical leads to the circuit board electrical conductors 
(Figure 9) . These connector pins also act as low resistance 
thermal conductors between the component and the circuit 
board. The total thermal resistance from the component to 
the circuit board may be calculated from the cross sectional 
area of each pin (Axpin) , the length of the pin (Lpin) , the 
thermal conductivity of the pin material (Kpin) and the 
number of pins (Npin) using the expression: 

Rpin ~ "(Kpin) x (Apin ) x (Npin) 

There is also a small gap between the bottom of the 
component case and the circuit board. Figure 9 shows the 
mounting for DIP cases where this gap is on the order of 
one mm; for chip carriers, this distance is smaller yet. 

Due to the small distance (Dgap) involved, it is assumed 
that no convection occurs in this gap and that the conduc- 
tive resistance of the air gap (Rgap-cond) may be calculated 
from the expression: 



Rgap-cond = (Kair ) V (Agap) 

where Agap represents the surface area of the bottom of the 
component and Kair represents the thermal conductivity of 
the air. 

The component case will also radiate energy to the 



circuit board. Assuming the gap to act as two parallel 



plates of equal areas with emissivities of Epsb and Epse, 
the heat transfer rate may be calculated (Ref [18] ) from: 



where : 



Q = 



(Sig) x (Agap) x (Te 4 - Tb 4 ) 

~ 1 , + 

Epsb Epse 



Sig = Stefan-Boltzman constant 

(5.67 x 10 -8 . Watts — 
m - degK 4 

Te = Component surface Temperature (degK) 

Tb = Circuit board Temperature (degK) 

This radiation term may be simplified by expanding (Te 4 - Tb 4 ) 
in a Taylor series about Te and retaining only the linear 
portion of the series. When this is done, an effective gap 
radiation resistance (Rgap-rad) may be calculated from: 



Rgap-rad = Epse + Epsb - (Epse) (Epsb) 

(r) x (Sig) x (Epse) x (Epsb) x (Te 3 ) 

Since the component case temperature (Te) is an unknown, this 

resistance will be recalculated as the solution progresses. 

These three resistances (Rpin, Rgap-cond, Rgap-rad) 
may be combined into a total resistance between the component 
and the circuit board (Rtote-b) . When the component is a 
chip carrier, the air gap and lead length are assumed to be 
one-tenth that of the DIP case [13] . 

Although the radiation heat loss from the component 
is extremely small and normally neglected in thermal analysis 
of electronic circuit boards, it is included here for completeness. 
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The component will radiate to the other components and the 
back of the adjacent circuit board. For the purposes of the 
thermal model under discussion, it is assumed that radiation 
is from the top of the component to the adjacent board and 
that this adjacent board temperature is the same as the board 
temperature below the component. It is further assumed that 
the area of the circuit board is much greater than that of 
the component. Using the same linearization technique pre- 
viously discussed results in the expression: 

Rtop-r = 1 - E P se =- 

(4) x (Sig) x (Epse) x (Ae) x (Tej 

for the effective thermal resistance for radiation from the 

top of the component to the adjacent circuit board. 

Convection heat transfer to the air, the final heat 
loss path from the component, depends on the heat transfer 
coefficient and the area for that mode of heat transfer. 

The component is assumed to experience two separate modes of 
convection. Those portions of the component that are perpen- 
dicular to the air stream are assumed to experience a stagna- 
tion form of air flow while the top and sides see parallel 
air flow. Recall that air flow is defined as being from left 
to right as required for data entry in the BOARDS program. 
When component types are entered as being horizontal or 
vertical using the special function keys, the component 
length (Le) and width (We) parameters are set such that Le 
is along the air flow direction and We is perpendicular to 
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the air flow. The standard correlation for plane stagnation 
flow may be written as: Ref [19] 



Hstag = (.57) x (Kair) 



(Pr* *) x ( 



Vair 



(We) x (Gnu) 



where Hstag = Heat transfer coefficient for areas receiving 
stagnation air flow 
Pr = Prandtl number for air 
Gnu = Kinematic viscosity for air 
Using the component height (He) , the area for stagnation heat 
flow (Astag) may be calculated from: 



Astag = (2) x (We) x (He) 

The portions of the component that are parallel to 
the air flow (top and sides) are assumed to experience the 
same heat transfer coefficient that would occur in a smooth 
duct modified by a roughness factor (Fr) to account for the 
presence of the components acting to increase this heat 
transfer coefficient (Hpara) . Hpara depends on the Reynolds 
number (Re) determined by: 

_ (Vair) x (Dh) 

Re GHU 

In the entrance region of the duct formed by the circuit 
boards, up to a distance of approximately ten hydraulic 
diameters, the flow is assumed to be laminar. Reference [20] 
contains the following correlation for the heat transfer 
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coefficient in the entrance region of a smooth duct: 

r (—) ij l 

Hduct = (- 664 ) x (Kair) * (Gz) x (l+(7.3) x [ Gz } I 
(1.1) x (Dh) Pr J 

where Gz is the Graetz number defined as: 

_ (Re) x (Pr) x (Dh) 

Gz xl 

and Xi is the distance from the entrance of the cooling air 
flow. 

The value of the heat transfer coefficient predicted 
from the above equation has been found to be low from compari- 
sons to experimental data of Ref [12] and Ref [14]. The 
expression : 

Fr _ i + 5 x (Perim - (2) x (Bdh) - (2) x (Zb)) 

Perim 



has been created to adjust the predicted heat transfer co- 
efficient for the rough duct. The expression reduces to one 
when no components are present and is less than two for all 
reasonable board constructions. The resulting heat transfer 
coefficient that is used for the parallel sides (Hpara) is 
therefore : 



Hpara = (Hduct) x (Fr) 

For positions beyond the entrance length the flow 
may be laminar or turbulent and the appropriate heat transfer 
correlation must be chosen. The transition from laminar to 
turbulent is assumed to occur at a Reynolds number of 1000 
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since the components act as turbulence promoters for the 
air flow. For laminar flow, the expression: 



Hduct = < 5 ~ 4) * (Kair) 

Dh 

is used and the result modified with the roughness factor 
(Fr) to obtain Hpara. For turbulent flow the Dittus-Boelter 
relationship is used: 

tt (.023) x (Re' 8 ) x (Pr’ 4 ) x (Kair) 

Hpara m 

Heat is convected not only off the top and sides of 
the components but also from the surface of the leads or 
pins. Since the thermal conductivity of the pin material 
is so high, their surface is assumed to have the same temper- 
atures as the surface of the component. The effective com- 
ponent parallel side area for heat transfer is, therefore, 
determined from: 

Apara = (Le) x (We) + 2 (Le) x (He) + (Aspin) x (Npin) 



where Aspin is the surface area of the pins that experience 
parallel air flow. The thermal resistances for convection 
from the component may then be calculated from: 

R P ara = (Hpara) x (Apara) and ' 

„ _ 1 

Rstag (Hstag) x (Astag) 
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These two resistances may then be combined as parallel re- 
sistances to give a total convective thermal resistance from 
the component to the air (Re-conv) of 



Re-conv 



(Rpara) x (Rstag) 
(Rpara) + (Rstag) 



Returning now to the basic premise that all the heat 
produced within the component must be transferred through one 
of the heat paths illustrated in Figure 8, the heat balance 
equation for the Ith component is: 



Pow (I) = Qto board + Qto air + Qrad off top 



utilizing the general expression for heat flow as a function 
of temperature difference. This may be expressed as: 



Pow (I) 



Te(I) - Tb ( I) . 
RToTe-b ( I ) 



Te(I) - Tair(I) 
Re-conv (I) 



Te(I) - Tb (I) 
Rtop - R ( I ) 



This expression may be rearranged and solved for the compo- 
nent temperature (Te) . The resulting expression is used to 
eliminate Te from the final set of equations that are solved 
for the set of board temperatures. 

4 . The Circuit Board - Thermal Model 

The circuit board has been divided into regions as 
determined by the user while entering the circuit board 
description in program BOARDS . Each of the regions will 
lose or gain heat by the normal processes of conduction, 
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convection and radiation. In this thermal model, radiation 
from the circuit board regions is neglected due to the lower 
temperatures of the circuit boards. The board temperature 
(Tb) is assumed to be uniform within a given region for 
purposes of convective heat transfer to the cooling air 
stream. For purposes of calculating the heat conducted 
between regions, the temperature differences are assumed to 
exist over the distance between the centers of the regions. 
These region centers are the nodes in the finite difference 
model used to analyze the circuit board. 

The surface of the board forms part of the air duct 
previously mentioned in the discussion concerning the calcu- 
lation of the convective heat transfer coefficient from the 
parallel sides of the components (Hpara) . The heat transfer 
coefficient of the board varies with distance from the air 
centered in the regions, the convective heat transfer co- 
efficient for the board (Hb) is assumed to be equal to Hpara. 
The area for heat transfer from each region (Areg) includes 
both sides but excludes any area under the component (Ae) 
since the air is assumed not to flow in the small gap between 
the component and the circuit board. The area (Areg) is cal- 
culated from the equation: 

Area = (2) x (Bdl) x (Bdh) _ 

Areg Nreg 

where Bdl = Length of circuit board 

Bdh = Height of circuit board 
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The convective resistance from the circuit board region 
(Rb-conv) is therefore: 



Rb-conv - -y-rrr— r r 

(Hb) x (Areg) 

The conduction of heat within the circuit board 
material and any added conductive material is assumed to 
occur between the nodes located at the centers of each 
region. For a given region, this heat flow is assumed to 
occur only out the four sides of the region. The distance 
between nodes in the horizontal or air flow direction (Lr) 
is calculated from the circuit board length (Bdl) and the 
number of regions in the horizontal direction (Nxr) . The 
expression used is: 



Lr 



Bdl 

Nxr 



Likewise, the vertical spacing betv/een nodes (Hr) is calcu- 
lated as: 



Hr = 



Bdh 

Nyr 



These distances are used with the circuit board thickness 
(Thb) and the circuit board material thermal conductivity 
(kb) to determine the base board conductive resistances in 
both the vertical (Rb-ver) and horizontal direction (Rb-hor) . 
The expressions used are: 

Hr 



Rb-ver = 



(kb) x (Lr) x (Thb) 
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and 



° r Tkb) x (Hr) x (Thb) 

The conductive resistances in the material added to 
the circuit board is treated in the same manner since these 
materials are also assumed to connect the nodes. The con- 
ductivity of the material (kl) and physical area (al) have 
been specified for either a conduction path unit (CPU) or 
electrical lead. These are used in the expressions: 



Rl-ver 



Hr 

(kl) x (Al) 



and 



Rl-hor 



Lr 

(kl) x (Al) 



to determine the thermal resistance for conduction in a 
single conduction path. This resistance is divided by the 
number of such paths to determine the total conductive re- 
sistance of the material added to the circuit board. 

This resistance is assumed to be in parallel with 
the conductive thermal resistance of the bare circuit board 
and a net conductive thermal resistance is calculated for 
each of the four directions using a normal product-over-sum 
formula for parallel resistances. Figure 10 shows how a 
typical region would thus be connected to the four adjacent 
regions. The values of these conductive thermal resistances 
are stored in a two dimensional array, Rl (I,J), where the 
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first index represents the region number and the second index 
is the direction number ( 1 , 2, 3, 4). These numbers represent 
bottom, right side, top, and left side, respectively. For 
example. Figure 10 shows how R1 (1,2) connects region number I 
and region number I + 1. For those regions on the edges of 
the circuit board, the conductive thermal resistance in the 
direction off the board are set to very high value due to the 
assumed adiabatic boundary conditions for all edges. The set 
of conductive thermal resistances, R1 (I,J). is used in the 
heat balance for the circuit board. 

With the component treated as the source of heat for 
the circuit board, one may again apply a simple heat balance 
for a region as follows: 

Qe-b = Qconv + Q x + Q 2 + Q 3 + Q 4 

where Qconv represents the heat transfer to the cooling air 
and through represent the heat conducted to the adja- 
cent regions as shown in Figure 10. This expression may be 
written for Region I as: 

Te(I) - Tb (I) = Tb (I) - Tair(I) Tb(I) - Tb(I + Nxr) 
Rtota-b Rb-conv Rl(I,i) 

. Tb (I) - Tb ( I + 1) , Tb (I ) - Tb(I-Nxr) 

(R1(I, 2) + R1(I, 3) 

. Tb (I) - Tb(I-j) 

+ — ri ( i ; 4 ) 

As previously stated, the heat balance equation for the component 
derived in Section 3 may be solved for Te(I) and that expression 
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used to eliminate Te(I) from the equation above. The only 
unknowns will then be the board temperatures. A heat balance 
for every region may be performed resulting in Nreg simultan- 
eous algebraic equations. The coefficients of these equations 
may then be collected into a matrix and solved using an LU 
decomposition method [21 ] . 

5 . The Analysis 

In performing the thermal analysis of the circuit 
board, there are two separate situations to be considered: 

1. Specified component power - determine steady state 
junction temperature; 2. Specified junction temperature - 
determine the maximum power allowable. For the first situa- 
tion a case temperature (Te) is assumed for each component 
since this is necessary to determine the effective radiation 
resistances. All thermal resistances and the local air 
temperatures are then calculated using the given variables. 

The resulting system of Nreg simultaneous equations is solved 
using LU decomposition for the temperature of each region of 
the circuit board Tb. 

This set of board region temperatures is used to 
determine the set of component case temperatures using the 
expression for Te determined from the heat balance of the 
component. These component temperatures are compared to 
those from the previous iteration. If the largest difference 
between any of the temperatures is less than the convergence 
criterion established by the user (default .1 degC), the 
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results are presented in the graphics mode on a fascimile 
of the circuit board as shown in Figure 11. 

In the display of the circuit board, each region 
contains the component type, junction temperature, power 
level and case temperature for the component in that region. 
Empty regions contain only the temperature of the circuit 
board. Those junction temperatures that are within 5% of 
the maximum are starred (**) for easy reference. This output 
is dumped to the internal thermal printer to provide a hard 
copy. 

If convergence has not been reached, the new case 
temperature is used to calculate new radiation thermal re- 
sistances and the new system of equations solved. Figure 12 
shows the intermediate display of all temperatures that is 
presented on the screen while the next iteration is in 
progress. If longer than 20 lines, the maximum display area 
for output on the screen, this display may be moved up or 
down with the display control keys at the center top of the 
keyboard (see Figure 2) . Convergence typically occurs in 
less than four iterations when solving for component junction 
temperatures . 

For the second situation with specified junction 
temperatures (Tj), a power level of .25 watts is assumed for 
each component. The component case temperature may then be 
determined from the expression: 



Te = Tj 



Pow 
Rj -c 
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where Rj-c represents the case to junction thermal resis- 
tance set by the user. As in the first situation, all the 
thermal resistances and the local air temperatures (Tair) 
are calculated. It is important to note that Tair depends 
on the assumed power levels since these values will change 
with each iteration. The resulting system of simultaneous 
equations is again solved for the board region temperatures 
and the component case temperatures. The resulting compo- 
nent power levels for each region are determined from the 
expression 



Pow 



Te - Tb 
Rtote-b 



This component power level for each component is compared to 
the results of the previous iteration or the assumed values 
for the first iteration. Convergence is assumed when the 
largest percentage difference from these comparisons is 
less than the convergence criterion established by the user 
(Default 1%) . 

The output is presented in the same manner as before 
except those power levels within 5% of the minimum are 
starred (**) in this situation (Fig. 13) . Since both the 
radiation thermal resistances and local air temperatures 
depend on the power levels, more iterations are required 
before convergence is reached. Typically less than six are 
sufficient . 
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6 . What Now Option List 



Following a complete cycle through THERML , there are 
several options available to the user. A different circuit 
board may be analyzed either by entering the descriptive data 
set from the keyboard or by retrieval from mass storage. The 
previously recorded data set may be read in using THERML with 
no data checking capabilities or BOARDS may be loaded under 
program control to allow data checking or changes to the 
circuit board. 

In addition the same circuit board may be reanalyzed 
with a new set of installation parameters, i.e. , board spacing, 
inlet air temperature and volumetric air flow rate per board. 

7 . Sensitivity Analysis 

Another option available allows the user to investi- 
gate the effects of various air flow rates. The user specifies 
a maximum air flow rate and five separate analyses are performed 
for each of five air flow rates up to the maximum specified. 
Figure 14 shows how the results of this sensitivity analysis 
are presented for the situation of specified power levels 
while Figure 15 is an output for the second situation of 
specified junction temperatures. The outputs from each 
analysis such as Figures 11 and 13 are not presented during 
the sensitivity analysis since the graphics mode is used but 
an output such as Figure 12 may be presented on the thermal 
printer if records of individual components are needed for 
each of the separate air flow rates. 
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A final option available to the user is termination. 
It is important to terminate the THERMELEX system under 
program control to insure normal key definitions are returned 
and graphics parameters correctly assigned. 
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111 * RESULTS AND CONCLUSIONS 



The THERMELEX system offers the designer of electronic 
circuit boards the means to predict the thermal performance 
of air cooled circuit boards and avoid the problems that 
often surface only after the equipment is in use. The 
system is easy to use yet provides valuable data in a 
variety of formats that can help the designer to make im- 
portant design decisions regarding circuit board layout and/or 
cooling air flow parameters. 

The THERMELEX system has been tested with a wide variety 
of circuit boards to insure that various combinations of the 
possible components and empty regions will create valid 
descriptive data sets and reasonable results. In all test 
cases the results have been satisfactory when compared to 
expected results. For example, higher component power 
levels result in higher junction temperatures and decreasing 
air flow produces higher temperatures with all else the same. 
Several tests were made that used the results of a power 
level prediction as input to the same circuit board to 
insure that predicted junction temperatures were the same 
as those specified for the original test. 

In addition to the above tests for general validity, 
direct comparisons to experimentally determined component 
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case temperatures reported in reference (14) were made. 

The circuit board consisted of 25 equally spaced 14 pin DIPs 
and is depicted in Figure 6. Three different air flow rates 
and four different component power levels were used. The 
results of these experiments are shown in Figure 16. These 
same flow rates and power levels were used as inputs for 
THERMELEX and the predicted maximum case temperatures are 
also included in Figure 16. At the lower power levels, the 
agreement with the experimental results is encouraging but 
for the higher power levels there are significant differences. 
In particular, the predicted maximum case surface temperature 
shows a much stronger dependence on the air flow rate than 
the experimental results would indicate. It is believed 
that several effects neglected in the the thermal model 
become significant for low flow rates and higher component 
power levels. For low air flow rates, the effects of natural 
convection heat removal become more important, thus decreas- 
ing the actual surface temperature. In addition, at higher 
temperatures, the conduction into the electrical connectors 
and mechanical supports will also tend to hold the surface 
temperatures lower for the experimental results. Further 
work is needed to resolve the differences between the 
experimental results and those predicted by THERMELEX. 
Particular attention towards refinement of the thermal model 
is required. 
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IV. RECOMMENDATIONS FOR FUTURE WORK 



While the present version of THERMELEX can be a valuable 
tool for predicting thermal performance of electronic circuit 
boards , improvements and refinements would be useful in 
several areas. The first recommendation would be for more 
experimental verification with particular emphasis towards 
developing better empirical heat transfer relationships for 
use in the present thermal model. The thermal model should 
also be expanded to include both the effects of natural con- 
vection and boundary conditions other than adiabatic. In 
particular the user should be able to specify a constant 
temperature for one or more of the physical circuit board 
boundaries. This would allow the modeling of installations 
that include metal card guides or cold plates. 

The present model does not include the interaction that 
can occur between circuit boards mounted closely together. 
These effects could be included in the radiation sink 
temperatures "seen" by a board or included as effects on the 
local air temperature for the bottom and top of the circuit 
board . 

Additional improvements could also be made in order to 
increase the number of circuit board designs that may be 
analyzed with THERMELEX. The limited set of components 
could be increased to include DIPs with different numbers 
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of pins and other flat pack case styles. Discrete components 
such as resistors, capacitors and switches are also able to 
transfer heat to the air stream and could also be included. 
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Figure 1 . The Hewlett-Packard 9845 desktop computer 
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Figure 3: General Purpose Keycode Overlays 
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Alternate Region Definition for Single Circuit Board 
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igure 4. 
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Figure 5: Component Definition Keycode Overlay 
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Graphics Dump of Circuit Board 
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THE ABOVE BATh IS FOR NELC/NF'GS TEST BOARD 
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Figure 7. Data List of Circuit Board 
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Figure 8 Thermal Network of Components and Circuit board 
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Figure 9: DIP Mounting on Circuit Board 
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CONDUCTIVE THERMAL RESISTANCE BETWEEN REGIONS 
Figure 10. Thermal Network of the Bare Circuit Board 
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OUTPUT DATA FOR NEZLOMPGS TEET BOARD 
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.600 W 




. 600 H 




33.7 C 
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.600 H 
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CIRCUIT BOARD DESCRIPTION IS STORED UNDER TEST 6 



gure 11. Final OUTPUT DATA from THERTiL . (Temperature Solution) 
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BRTA FOR NELC/NPGS TEST BOARD 
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igure 12. Intermediate OUTPUT Data from THERML 
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OUTPUT CflTfi FOR NELC/J1PCS TEST BOARD 
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Figure 13: Final Output Data From Therml (power solution) 
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Figure 14. Sensitivity Plot: TJUNC vs Air Flow 
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Figure 15. Sensitivity Plot: Power vs Air Flow 
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MAXIMUM CASE TEMP (deg C) 




Figure 16. Comparison of THERMELEX Predicted Case Temperatures 
with Experimental Results from Ref. 14 
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The THERM EL EX system is designed to perform thermal analysis of air cooled 
sctronic circuit boards. The system consists of three major programs: 
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The screen will now shift from the alpha mode to the graphics mode and draw 
•' ® « P * y board w i t h t h e i n. ag 1 n ar y r e g i o n s m ar k e d o f f w i t h dotted 1 1 n e s . This 
cture will remain on the screen for only 3 sec after it is drawn but it may 
returned to the screen using 12 to allow further time for checking. 
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When -all regions hawe been specified, the screen will return to the alpha 
mode and ask if any corrections are needed. Use the special functions keys 
as required to either return to gr aph i c al pk. a or produce a hard copy. Enter 
t h e n u m b e r s o f t he c o n n e c t i n g r e g 1 o n s t h a t h a<< e 1 n c o r r e c t v a lues for t h e r. u m bers 
of leads or CPU's. hs before, CGNT with no entry signals the lack of further 
correct tons. The physical description of the circuit board is now complete. 
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be permanent 1 y 1 ost . I nf orm at. i on c one e rn i ng t he m i n i mum dat. s f i 1 
p r e s e n t e d b y t ki e p r o g r am t o k. e 1 p i n t ki i s. d e c i s l o n . 



«- >y o O ~ Ml 

•li x Cl* x - *i' 

3 t> - ,ti 

•fl «-• - 3, 3 3 3-. 

Ml 3 — *3 hi -H - . 

O *0 *TI '3. X i.i 

— X 'i* O a* 



O m i i* '1 T) ►' 3 

•i> —r — IZ 3 Z 

4- ■« Cl 3 v 3 <i> 

Ml w - W Ml 

3 . i_ >11 • O 3 Ty o 

•'l O _Z "I X 1*1 rr> X 

_ <3 x 3 *rt -• 

X o .« ^ Ti 

UJ IJ ■— "I x 

z: • a* <y i> mi >it 

LL '•! TJ _ *1* 'll Ml • J 

U J. i. S. e "l 

X o 3 o x ■ T i «it 

i. iTl 

ill ill ill ill Hi _ _ 

3 3' x O — ' J i. 

— 3- « <Tt W — -il = 



■ii a* ■!» o 3 •!' o — 

w £ L L IJ *1 

•li ill 3 n' 3 

X u 3 ' i* x X - — 

x O — 'J 3 — 

3 4. jj 
*> — i. Ill 3 O 3 

3 a* 3 O 

o '•* i< o — 3 

X «- H> X — »-• 3 HI 

U> 3 X x ij x ^ 

3 X «- 3s ''y 3 M l 

o c <T1 3 M l *1* ; 

3 •- CJi S >n 



3' 3 
f_ n> z O X 

3 Ml 111 O ■«- 

Ml o 3 M l •_ -n H I 

ItT tj o ij iy 

Z a* Z ‘U *■ 'j Z X 

w 3 <n jZ 2* v Z 3 

i. in _i. *" ° = 

O <*■ - n 3* 'li . 

Z O •» 3 h- x z •_ 

3 -z 111 X '« - 

'll Ml *1' I'll O *3 — *-■ 

M I M l 3 O — 3 . 

3 C 'I' u 'l* a* X • 

•<- — Z Ml 3s X X 



— o mi x n . o 

— 3 — X 'll "l X 'll 

3 3- —iv Z 

w Z 3 *tl 3 i) 

Z 0 0 3 3 3 

<0 3 — 3* m i O 

Ml Z 'll — Hi _ 'll 

3 ' 1 ' j j "i Z 3. Mi 



TTI ill iTl |-| 

•n - Jj 

'll jj JU 



I — ' I Z 'll O 

Ml UL I 3 - i_ IJ 

■X .XI _ '!■ ^ 'll 

• W - - .rt 

'll JZ - T) Z- 

m - X 3- O 

— 3 3 - "I _J | 

3 "I 3 IJ Ml UL | 

o z 'I- i' in :-s ujj 

i.i o "O •— jZ *— i 



01 



•li u. ill ^ 

*. O 3 



w- LU 

^ r •'! Ill _J 

~ JZ Ui 



rn 'X 'll • - >n = UJ 

: I 3 'II i u I 

O X M 31 IJ h- 

jE - Z Z 0 u. Hi 
«- Ml .71 O 3 Cli X 



IJ — in Ml 3 

•3 3 'Ti - >n 

ITl ITl z -■ •- 

•i' n* o >n c 3> 

■i- i3 J 3 O O 



•Tt — 



•3 

■li — 



*- 1' *> X 

JZ — o TI 



O h- >il O LU 

w - 3 LU 

>Tl • %- - 

Z t. Ml 

'Tt Hi hi tj Ml 

— TI X HI 'll 

'3l Z •-> Z Z 

’•■■■ O '!• X 
■i' 3 «-> 

■3 O '3 ’3 



xj - — 

LU V|_ 'll XI 

t— I iTl 111 HI = 

*-'■ 'Ti X ^ 

•Z 'Tl ». 

— Tl i-j 3. Uj 
■:o| -n — 
3 -3 O 5 I 
o 'j H 

' — "3 Z'| c I.'*' 

x s_ iXI in in 

3 u. X 

= IJ X 3* X 

.fl r v o n .3. 

Z Z Z mm 

3« _ ' I' X = 

X mi " o|ir 
<v O X 'I J 

4-2- Mil 
X Z 3* 'i* 

o "Z .. 



•n hi ij hi o hi 

X " I Cn H' 3 X 

K 3 i. TMI j 



fl ft jj •" *- 

x - 3 

i. O 'Tl 

O £ 1' Z 

X ' i> Mi h» 

'I' mi Z >v *- 

£m 3 0 X 3 

■i' £ 

«-> O O KJ 



: 'ii <n 

iy 'I' 

: to — 



z •- z ** 

- £ -ii 3* — 3 

'Tl O iy hi X 

Z Z Z Ml 

•n x i_i >z ■— "i 

*□ Ml ' 1 * 0 i 

M- X — T) X' 

•y in iy «- 3' i_ X 

3 mi X 'i' |T t 

— ♦- I s Z i j 

Cl *— -3 3 * X " U O 

■- — O L* O 3 

i. lit IJ 111 l.i'l 

"1 Z 3 Z ^ 'rt • 

H' ij ill *:• X 'rt iy 

T> 3 Z 'rt £ — 



C* X Z 



■ M- 3 
'rt X Ml 'i' 



— X 'J LI 



X — UI 
— Ill *T 

U 3 CO 



Ml *- — Cm- O 

•3 'rt — iy *3 

- X - 3 'li 

fl j 3 o 3' 'i' 

x X 3 X 

3 mi *3 = i i.i 

O 'I' Z O X - 

o — irt -j u 

X O 

x 'Tl X U 0 £ 

•rt - — x 

j- 3 u — 'rt . 



“"'I i- J 

TJ 'I' 



' X Ij x «-• 



O TI ' ll 
- iy £ 
*_ !_ irt 



t_ |TI X O O O H* *3 irt X Ml o 

i I O 3 "3 x '3 ’3 'i' iy 'rt x 

mi X 3 Ty 3 • 3 — X Ml 'I' 

.rt .rt x «rt tj ■■! ui o mi iy «-■ hi mi "i jt 

£ 'y ZT x i'y il' >- •— X Z 'rt 3 'j 

£ *3 3 3' Ml iTl i j ij 0 'rt 

hi irt rn x <rt — 3 «-• — mi *— *- 3 - X 

w 3 3 O 'ii X -rt ■— j- hi 7J "3 ^ 

.rt 3 'ii X iy — hi in ij z 'J Z 

■— 'li £ x *3 ij ... z j . X Hi 3 'rt iy ij 

X --irt XX 'i' TJ 'y 

x 3 3 x 3 : • 

3 3i — — 'rt 'y H' — — — 3 '!' 

X 'ii O — — T- irt Mi 3 'rt 'Tl 

rj . — ij O £ irt irt ij TJ 

X Z — *3 X mi Mi mi 3 3 'it 

x j; 0 mi 3 X •— a* 'rt *3 

£ - 0 0 X 3. 

3' — *3 • X x — m* £ 3' O 

3 •— 3 Ml '!• X 3> >— •!' 3 'J 

■4- x ill 'J *3 3 Ty •<- ■>- Ty 

3 O - 'rt O 3 3 •- x 3 

•j 'rt X X £ • •— Ml 3 U 

. — x X j . X 3 Ml s O 'll 'j 
>— iTl Tj ifl ij (Tl ij iy IJ j -r. x X 

O Tl i 3 X 3 '!• 'll 3' X «-• IJ 

X 31 3 3 3 X' 'J 3 



i*J X U_ 

iy 31 *D 

X — irt 

3 ^ 



IJ Tj 

O 3 O 



'J 



'rt 

£ irt O X O ci X X 



75 





76 



The final results art written onto a picture of the circuit board with 
each region containing component type, junction temperature, power level and 
cast t e f- p e r a t u re. E n. p t y r e g 1 o n s c o n t a i n o n 1 y t h e t e rn p e r a t u r e of t he circui t 
board. Those junction temperatures within 5-. of the maxi muni and those power 
1 e v e 1 s w i f h i n 5 *. o f t h e rn i n i rn u m ar e s t ar r e d O * > for e a s y r e f e r e n c e as t r o nb 1 e 
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APPENDIX B 
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AUTOST DRIVER PROGRAM FOR 

thermelex 

A SYSTEM OF PROGRAMS FOR THE HP 9345 
TO PERFORM THERMAL ANALYSIS OF ELECTRONIC CIRCUIT BOARDS 



PREPARED AT THE NAVAL POSTGRADUATE SCHOOL MONTEREY 
R. A. FOLTZ LCBR USN 



■ f 4 -**.*»* 






♦ | | | | | AUTOST | | | | | JUNE 38 * 

* * 

PRINTER IS 0 

PRINT C H R I *. 2 7 ' .■&"?< 1 3 4 T " 

DIM S-TC60] 

COM Map 
PRINTER IS 16 
GOSUB Setst o 
GOSUB Ir.it 

IF Ana l = " BACK _UP" THEN GOTO 248 
GOTO Start 



290 St arte: GOSUB Error 

300 Start: St - a NOTE TO USER ABOUT SPECIAL FUNCTION KEYS" 

310 ON ERROR GOTO Start-: 

320 GOSUB Page-head 

330 PRINT " The following ussr-keyi will be in effect unless otherwise note 
d. " 

340 PRINT "In general these key def f i n i t. i ons are valid on 1 ■ 1 when the program 
i a uai 1 1 ng " 

358 PRINT "for an input < i e Run light in lower right corner ia off;. The Back 
up option" 

360 PRINT "CkO; returns the program to a previous question or control point." 
370 PRINT 

330 PRINT "This should be used when the user decides he h-s made some error on 



i r.put 
390 
4O0 
410 

6164)“ 



t hat ’ 
PRINT 
PRINT 
PRINT 



"car. be corrected by re-entering the data." 

"If you do not have a plastic key code overlav filled in <hp # 7120- 



428 

438 



PRINT "a paper one will be provided." 
PRINT L I N < 1 ) 



440 


PRINT 


" 


KEY 


' k 0 ' 


will 


back up to the previous question." 




450 


PRINT 


" 


KEY 


'k 1 ' 


will 


display print, ed output. 


(EXIT GRAPHICS’* 


460 


PRINT 


" 


KEY 


'k2' 


wi 1 1 


d i sp 1 ay p 1 ot » ed out put . 


(GRAPHICS)" 


470 


PRINT 


" 


KEY 


'k 3' 


will 


pr i nt c opy of pi ot t ed out put . 


(DUMP GRAPHICS’* 


480 


PRINT 


" 


KEY 


' k 4 ' 


will 


disable internal printer. 


(PRINTER 


IS 16.* 


490 


DISF " 


P r e s a 


C0NT 


WHEN 


i READY FOR MORE EXPLANATIONS" 






500 


PAUSE 
















510 


PRINT 


" 


KEY 


'kZ' 


will 


enable internal printer. 


(PRINTER 


IS 8>" 


520 


PRINT 


" 


KEY 


' K6 


will 


rewind left hand tape. 


(REWIND 


: T 1 4 > 


530 


PRINT 


" 


KEY 


' K 7 ' 


will 


rewind right hand tape. 


(REWIND 


: T 1 5 > 



( 1 > 
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AUTOST 



DRIVER FOP THEPMELEX 



546 PRINT 
550 PRINT 
560 PRINT 
570 PRINT 
CONT - 
530 PRINT 
ion in" 



KEY "k 14 ■'will answer Yes." 

KEY 'Kl5'wi11 answer No." 

LIN (2) 

"In general when answering questions or 
"with no other entry will assign either 



590 PRINT "parentheses (the first option in the 1 
600 PRINT LIN( 1) 

610 PRINTER IS 16 

620 COSUB Overlay 

630 IF ftnsl= " B ACK_JJP " THEN 260 

640 DISP "LOADING^ EY DEFF I N I T I GNS " 

650 IF M ap = 0 THEN LORD KEY " EDSKEY " 

660 WRIT 1000 

670 GOTO 690 



selecting an option pressing 
the default or the first opt 
st. > as the desired option." 



6S0 Optionse: GOSUE Error 

690 Options: Si="MRIN PROGRAM OPTIONS" 

700 ON ERROR GOTO Optionse 
710 GOSUE Page he ad 

720 PRINT " The THEPMELEX program package consists of thre 
one of" 

730 PPINT "which is in core now providing these instructions, 
re cal led" 

740 PRINT " B 0 R R D S and T H E P H L . Both have the capability to read circuit to 
ard data" 



major sections, 
The other two a 



750 



PRINT "from 



Ma 



Storage Device <. eg Tape or Disk) but only BOARDS can: 



760 

770 

730 

790 

SCO 



PPINT "<a> Input circui* beard des c r i pt i : from keyboard...": 

PRINT "(b) Allow graphical check ingot descriptions.."; 

PRINT "(c) Make corrections to aescri pt ions as required. . . "J 
PRINT "(d) Record descriptions to mass storage..." 

PRINT L I N 1 ) , " THEFML performs the thermal analysis and output of results 



810 
820 
830 
E" 
840 
GE" • 
850 
860 
870 
880 
890 



PRINT L IN*, 1 ) , "YOUR OPTIONS R P E : " ; 

PRINT " 1. LORD BORPDS TO INFUT DESCRIPTION FROM k E v BOARD . " 

PRINT LIN(ly,SPA(20>, "2. LORD BORRDS TO READ DESCRIPTION OFF MRSS : 

PRINT LIN( 1 ■ , SPR(20 s "3. LORD THEFML TO PERU DESCRIPTIONS OFF MASS 
WITH IMMEDIATE ANALYSIS." 

PRINT L I N ( 1 > , SPA ( 20 , " 4 . COPY ALL PROGRAMS " 

An « f = " 1 " 

INPUT "YOUR CHOICE *“ < 1 . 2 . 3 , 4 > " , Ansi 
IF Rns-?= " BACK UP" THEN Start 
Ans = INT (VAL ■ Ar.sf ) > 



900 IF ( Ans >0 > AMD 'Ans'5) THEN 930 

910 GOSUE Err in 

920 GOTO Options 

930 IF Ans= 1 THEN Map =2 

940 IF Hns = 2 THEN Map = 3 

950 ON Ans GOSUE’ Ge r. , Ge n , Therm 1 , Copy 

960 GOTO Options 

970 Gen: DISP "LOADING GENERAL CIRCUIT BOARD PROGRAM " 

980 LOAD "BOARDS", 1 
990 STOP 

1000 The r ml : DISP "WORKING LOADING THEFML " 

1010 LOAD "THEFML", 10 
1O20 STOP 



iTORAG 
ST ORA 



79 



AUTOST 



DRIVER FOR THEFMELEX 



1030 Copye:GOSU£ Error- 

1040 Copy: ! THIS SECTION FOR PROGRAM REPRODUCTION ONTO ANOTHER MASS STORAGE 

1050 ON ERROR GOTO Copue 

1030 $1 = "COPY I NG THERHELEX PROGRAMS" 

107O GOSUB Pagehead 

1030 PRINT " This section allows easy reproduction of THEFMELEX from one mass 
storage" 

1090 PRINT "device t.o another. The mass storage medium (taps or disk) must be i 
n i t i al i zed . " 

110O PRINT "It the medium intended as the destination is new or you wish to eras 
ae all" 

1110 PRINT "files then the use of option 2 will perform this task. Since the TH 
ERMELEX" 

1120 PRINT "system requires only appr ox i rn at 1 y 500 records.it is possible to use 
a medium" 

1130 PRINT "that already contains files to be retained; however, this may not al 
low" 

1140 PRINT "sufficient space for the data files containing the circuit board des 



criptions." 

1150 PRINT "Option 1 will record in the available spaces if possible but it is r 
ec ommende d " 



1130 PRINT "that an entire mass storage meriia be devoted to THERMELEX. " 
1170 PRINT LIN< 1 , "YOUR OPTIONS APE:"; 

1130 PRINT " 1. COPY ALL PROGRAMS WITHOUT INITIALIZATION." 

1190 PRINT LIN< 1 ) , 3PA<20>, "1. COPY ALL PROGRAMS AFTER INITIALIZATION." 
1200 Ans I = " 1 " 

1210 INPUT “YOUP CHOICE FROM ABOVE ( 1 or 2 ) " , A n s I 
1220 IF Ans 1= " BACK-UP " THEN Options 
1230 Ans= INT< VALv Ans I ■ > 

1240 IF < Ans= 1 ) OP (Ans=2> THEN 1270 



1250 GQSUB Err in 
1230 GOTO Copy 
1270 Ans I =Msus 1 I 

12S0 EDIT "WHAT MASS STORAGE DEVICE CONTAINS THE THEF.MELEX SYSTEM ( : T 1 5, : T 1 4 , : FS 
, ETC > " , Ansi 

1290 IF Ansf=" BACKJJP" THEN Copy 

1300 IF AnslC 1 , 1 3->": " THEN Ansl = " : " &Ansl 

1310 Msusll=Ansl 

1320 Ansi = " : T 1 4 " 

1330 EDIT "WHAT MASS STORAGE DEVICE CONTAINS THE DESTINATION MEDIUM CT14,:T15,: 
F8, ETC)", Ansi 

1340 IF Ans 1= " BACK _UP" THEN Copy 

1350 IF AnsIC 1 , 13' : " : " THEN Ansl= " : " ^Ansl 

1330 IF Ansi- Msusll THEN 1450 

1370 BEEF- 

1330 PRINT PAGE , L I Nt 1 0) , SPA ■; 1 O :• , "HOW CAN I COPY FROM " ; Hsus 1 1 ] " TO ";Ansl 
1390 PRINT SPA (10) , "TRY AGAIN" 

1400 WAIT 2000 

1410 BEEF- 

1420 WAIT 15O0 

143G BEEP 

1440 GOTO Copy 

1450 Msus2l = Ar.sl 

1430 IF Ar.s=l THEN Copy2 

1470 BEEP 

1430 PRINT PAGE, LIN 1 ' 10>, TAB< 15), " 

1490 PRINT 

1500 PRINT TAB'. 15) , “*,*,.-**♦* 4***-**- WARN I NG *■**♦■•■■■■****-*■••-** ■* " 
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1510 PRINT L I M\ 2 • 

1520 PRINT TREC 15-' , "THIS PROGRAM WILL ERASE ALJ_ FILES ON " :Msus2f 
1530 PRINT L I N < 2 1 , T A E ^ 1 5 , " 

1540 PRINT L I N <. 1 > , $ P A 10>,"USE i evQ TO ■ BACK UP IF THIS IS HOT YOUR DESIRE" 

1550 PRINT L I N ( 1 > , SPA < '20 • , "PRESS CONT IF THIS IS OK " 

1 560 WAIT 200 
1570 BEEP 
1530 INPUT Ansi 

1590 IF Ans i = " EACH _UP " THEN Copy 

1600 DISP "WORKING INITIALIZING ";Msus2l 

1610 INITIALIZE M s u £ 2 i 

1620 IF Msus2iC2, 2] . >"T" THEN Copy2 

1630 DISP "WORKING SETTING UP ";Miu-2l 

1640 CREATE " DUMMY *' ".Ms us 2i , 420 

1650 CREATE " ENDSPA " :s.Ns us 2i , 1 

1660 PURGE " D U M M Y " & M su s 2 f 

1670 Copy2 : ! IN THIS SECTION THE ACTUAL COPYING TAKES PLACE 
1630 RESTORE 1750 
1690 FOR 1=1 TO S 
1700 READ Hami 

1710 COPY Naftii&Msus 1 i TO Nami?Msus2i 

1720 DISP "COPYING "jNamei;" FROM "jMsusli;" TO ";Msus2i 
1730 WAIT 1500 
1740 NEXT I 

1750 DATA BDSKEY, STDr EY , DEMO-T . DEMO-P , AU TOST , BOARDS , THEPML , THAMES 
1760 RETURN 
1770 ' 

1780 Set st os- .* GOSUB Error 

1790 Set st c- : ! THIS SECTION DETERMINES WHERE THE THERMELX SYSTEM IS AND SETS 
1300 ' MASS STORAGE TO THAT PLACE 

1810 ON ERROR GOTO Setstoe 

1820 Si = " LOC AT I ON OF THE THERMELEX SYSTEM" 

1830 GOSUB Page he ad 

1840 PRINT " It is necessary that the mass storage media that contains the 
ERMELEX" 

1 S50 PRINT "system be located in the default mass storage device. Therefore, 
e default" 

1860 PRINT "mass storage will be converted by the program. This should not b 
1 1 e r e d " 

1370 PRINT "during the use of THERMELEX. Data files located on other mass st 
ge devices" 

1830 PRINT "may be accesed by appending the appropriate device code to the fi 
n ame . " 

1890 Ans i= " : T l 5 " 

1900 EDIT "CHANGE OP ENTER BELOW THE LOCATION OF THERMELEX ?" , Ansi 
1910 Ansi = UF'Ci< Ansi > 

1920 IF AnsiC 1 , 1 ] • >" : " THEN Ansi = " .* "$<Ansi 

1930 IF (Ans i C 2 , 2 ] = "F" > OP < Ansi C2, 2]="T" > OP t Ansi C 2 , 2 ] = " Y " > OP < Ansf C2. 2J = 

> THEN 1970 

1940 BEEP 

1950 GOSUB Err in 

1960 GOTO Set s t o 

1970 Msus 1 i=Ansi 

1930 MASS STORAGE IS hsus l I 

1990 RETURN 

2O0O Overlay: ! THIS SECTION PREPARES AN OVERLAY FOR THE USER 

201O Ansi= " Y" 

2O20 INPUT "DO YOU DESIRE A PAPER KEY-CODE OVERLAY < V OR 10? M .Rnsi 
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2030 IF vURCf 'final '*"H" > OP <flnsl*"BP.CK UP" > THEN RETURN 
2040 PRINTER IS O 

2050 PRINT SPAU5) , "TEAR OFF FOR GENERAL PURPOSE KEY CODE OVERLAY ",LIN<2> 

2060 PRINT " | Eac k up | EXIT | ENTER ] DUMP [~FRT I S | F R T IS | REWND | REWND | 

PLACE" 

2O70 PRINT "| ,ko> | GRAPHIC | GRAPHIC} GRAPHIC | 16 | 0 } : T 1 4 | :T15| 

ABOVE LEYS" 

2030 PRINT "TEAR HERE THEN PRESS CONT " 

2090 PRINT L I N < 1 ) 

2100 PAUSE 

2110 PRINT "j (13) | | | | j EDIT j YES | HO | 

PLACE" 

2120 PRINT "| | | | j | LINE | | | 

BELOW KEYS" 

2130 PRINT LINO) 

2140 PRINTER IS 16 

2150 RETURN ! END OF Over law 

2160 ! 

2170 P-ageheaa: ! THIS SECTION PLACES DESIRED HEADING ON A BLANK CRT 

2130 PRINT PAGE, TAB < 34-LEN < S-f ^ :• . " *■* " ; CHRI < 1 32) ; Si; CHP f ■. 1 23) J " *^",LIN<2) 

2190 RETURN 

2200 Err i n : ! THIS SECTION ALERTS THE USER TO AN ATTEMPT TO INPUT INVALID DATA 

2210 BEEP 

2220 PRINT PAGE 

2230 DISP “INPUT OUT OF RANGE TRY AGAIN" 

2240 WAIT 25O0 
2250 BEEP 
2260 WAIT 1500 
2270 BEEP 

2280 RETURN < END OF Err in 

2290 Error: ! 

23O0 PRINTER IS 16 

2310 PRINT LIN' 10;, TAB'20’ , "EFPOP NUMBER " ; EF.PN ; " HAS OCCUPED IN LINE "JERRL 
2320 DISP "PFESS CONTINUE WHEN READY TO RESUME PROGRAM FLOW" 

2330 PAUSE 
2340 ! 

2350 PagsertS! THIS SECTION BREAKS THE INSTRUCTIONS INTO PAGES FOR THE CRT 
2360 DISP "Press CONT when read 1 ; for more initrctions" 

2370 PRINT LlfHBl ar.h s> 

2380 PAUSE 
2390 RETURN 
2400 ! 

2419 P ageprt J ! THIS SECTION BREAKS THE PRINTED INSTRUCTIONS INTO PAGES 

2420 PRINT LIN<B1 anks ) 

2430 PRINT TAB< 35 ■ , "Pg. "JPagenum 
2440 PRINT LINv 1 ) , " 

2450 IF P agenum >5~THEN RETURN 

24 60 PRINT L I N ( 4 ) , TAB < 37-LEN< S* ) - 2 ) 5 CHP l < 132 ) ; 3 1 \ CHPf <1 28 .■ ; TAB <76) ; "Pg. " ; Pagenum 
+ 1 , L I N < 2 ) 

2470 RETURN 
2430 » 

2490 InstelGOSUB Error 

2500 Inst : i THIS SECTION PREPARES A SET OF WRITTEN INSTRUCTIONS FOP THE USER 
2510 ON ERROR GOTO Inste 

2520 $*=" INSTRUCTIONS FOP USE OF THEPMELEX" 

2530 GOSUB Pagehead 

2540 PRINT " fi written set of instructions can De prepared for the user that 
will help" 
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2550 PRINT "on*? to become familiar' with th? THERMELEX ij'item without actually us 
i n *3 the " 

2560 PRINT "the computer. These are intended to only supplement the set of dire 
c t i ons" 

2570 PRINT "given during the e ecution of the system of programs. These instruc 
t i ons will" 

2580 PRINT "normally be presented on the screen but if you desire a hard copy :*f 
the" 

2590 PRINT “ instruct ions, ENTER YPPT rather than Y. - 
2600 A n s X = " N " 

2610 INPUT "DO YOU DESIRE PPINTED INSTRUCTION'S* v N or V OP YPRT ) ? " , RnsX 
2620 IF flnsf = " EfiCh_UF " THEN RETURN 

2630 IF <UPCf 'finsf Cl , 13)»"Y"> OP » UPCf < Fins X '■ = " N " THEN 2660 
26-40 GOSUB Errir, 

2650 GOTO Inst 

2660 IF UPCf •' Rns X ^ = " N " THEN RETURN 
2670 PRINTER IS 16 
26S0 PRINT PAGE 

2690 IF UPCXCAnsX C 1 , 2] >="YP" THEN PRINTER IS 0 
2700 P= 1 6 

2710 IF UPCX'RnsXC 1 , 2] >="YP" THEN P=0 

2720 PRINT LIN<4> , TAB 1 , 37-LEN 1 . SX 1 2.-, “ ; CHRX <. 1 32 > : ST I CHF'f ( 1 23 '* ; " t*" f LINC2> 

2730 PRINT " The THERMELEX system is designed to perform thermal analysis of 
air cooled" 

2740 PRINT "electronic circuit boards. The system consists of three major progr 



Driver and system setup." 

Input and recording of data." 
Analysis and output of results," 



am s : " 

2750 PRINT 

2760 PRINT LI N< D , TAB* 1 7> , " l . AUTO ST 
2770 PRINT LIN'D, TAB* ID, "2. B OAF PS 
2730 PRINT LIN'D, TAB<17>, "3. THERML 
2790 PRINT L I N < 1 > 

2800 PRINT "The A U T 0 3 T program establishes special key deff ini t ions and acts a 
s the driver" 

2310 PRINT "for the other programs within the THERMELEX system. The key deffini 
t i ons are " 

2328 PRINT "valid throughout the system except when entering the specific comp on 



ent type" 

2330 PRINT "codes during execution of BOARDS . " 

2840 B 1 ank s = 0 

2850 IF P= 1 6 THEN GOSUB Pagecrt 
2860 PRINT 

2870 PRINT " This driver program is entered by placing the t ape in 
y tape" 

2S80 PRINT "transport <:T15) at the riqht side above the special keys. 
: ", LIN<1) 

2390 PRINT TAB<10>, “TYPES LOAD " ; CHR X ■. 34 > ; " AUTOST : T 1 5 " J CHF'f < 34 ; ; " , 1 
SS EXECUTE 
2900 PRINT 

2910 PRINT "If the 9345A is off, the driver may be loaded by insert, ing 
in ; T 1 5 , " 

2920 PRINT "latching the AUTO ST Key in the down position, and turning 



he primar 
You then 
Then Pre 



the tape 



the power 



on. " 

2930 PRINT "In either case the first question will be concerning the location of 
the " 

2940 PRINT "THERMELEX system. The program will convert the default mass storage 
devi ce" 

2950 PRINT "to what ever device contains the system. It is important to not c ha 
nge this" 
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29c-0 PRINT "while using THERMELEX. These instructions may be provided on paper 
or CRT and" 

2970 PRINT "a paper keycode overlay may be provided as a guide to the key deffin 
i t i ons" 

29S0 PRINT "that is loaded by h U T 0 9 T . R plastic over 1 ay is available from Hew 
1 e 1 1 - P ac k ar d " 

2990 PRINT "<part # 7120-6164) which mav be filled in using the paper overlay as 
a guide." 

3000 PRINT "The other proqrams are entered from the option selection listCMenu) 
in RUTOST . " 

3010 PRINT "Normally B 0 R P D 9 is next loaded to allow the input of circuit board 
d e s c r i p t i v e " 

3020 PRINT "data from the keyboard. The instructions that follow apply to this 
opt i on. " 

3030 PRINT 
3040 E 1 ank s = 4 

3050 IF P= 1 6 THEN GO SUB Pagecrt 

3060 PRINT " The user should have either the circuit board or a facsimile av a 
ilable" 

3070 PRINT "while he is using the EORPDS program. Orient the board such that 
the cooling" 

3030 PRINT "air flow is f rom the left to the right. Me asure the length < defined 
in the air" 

3090 PRINT "flow direction), width and thickness of the circuit board in mil lime 
ters (mm) . " 

3100 PRINT "It i j necessary to divide the board into evenly spaced imaginary reg 
ions in both" 

3110 PRINT "horizontal and vertical directions. The centers of each of these re 
g i o n s will" 

3120 PRINT "form * he nodes of a finite difference grid to be used in the thermal 
anal ysi s" 

3130 PRINT "of the electronic circuit board. The imaginary lines that are creat. 
ed divide" 

3140 PRINT "the circuit board into regions that will either be empty or will con 
tain one" 

3150 PRINT "component. Rn exception to this rule is the case of 40 pin DIP'S wh 
ich occupy" 

3160 PRINT "two regions due to their large size. Rn attempt to approximate the 
real circuit" 

3170 PRINT “board with a model requires some flexibility and imagination at this 
point. “ 

3180 PRINT 

3190 PRINT " The thermal model used for the analysis will assume that each co 
mponent is" 

32O0 PRINT "centered in the region. The following components can be included!" 
3210 PRINT 

3220 PRINT " DIPS (horizontal or vertical > . CHIP CRPPIEPS " 

3230 PRINT " 14, 16, 24, 40 pi ns 16,24,40,64 pins" 

3240 PRINT 
3250 Pagenurn=l 
3260 B 1 an h s = 4 

3270 IF P=0 THEN G0SUB Pageprt 
3280 IF F' = 1 6 THEN G0SUB Pagecrt 

3290 PRINT "The limit it ions imposed by restricting the components to the list ab 
ove will " 

3300 PRINT "prevent some circuit boards from being precisely modeled using THEPM 
ELEX but " 

3310 PRINT "are necessary to provide for the simple input o f the descriptive dat 
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3320 PRINT 

3330 PRINT " When BOARD'S is loaded ustng option #1 in BUT 03 T a control pa 
ramete-r is" 

3340 PRINT "assigned which causes the program to begin at. the section which asks 
Tor the" 

3350 PRINT "geometric parameters of the bare circuit board. You will be asked t 
o enter the" 

3350 PRINT "length of the board in millimeters. Enter the length using either t 
he number" 

3370 PRINT "pad at. the f a>' right of the keyboard or the numbers above the letter 
keys" 

3330 PRINT "followed by the C0NT Key. The next question will appear at the bo 
ttom of the" 

3390 PRINT "screen while the length that was just entered will be printed irnmedi 
at 1 y below" 

3400 PRINT "the instructions. Answer each question remembering to use the CQNT 
Key and" 

3410 PRINT "check the values on the screen. If a mistake is made press K0 and 
the program" 

3420 PRINT "will ask the question again. After the length, height, thickness, c 
onduc t i v i t y , " 

3430 PRINT "number of regions in the 'X' direction and the number of regions in 
the 'Y' " 

3440 PRINT "direction have been entered, the program will ask for a short (less 
than SO " 

3450 PRINT "characters > descriptive title. This title will become part of the d 
at a set for" 

3450 PRINT "the circuit board and will be used to identify the output. Two exam 
pies are! " 

3470 PRINT 

34S0 PRINT SPA(5) , " TE:T BOARD FOR PROJECT UMPTYFPhTS or DESCRIPTIVE TITLE 
41234. 9A . " 

3490 Elanks=0 

3500 IF P = 16 THEN G0SUB Paqecrt 
3510 IF P=0 THEN PRINT 

3520 PRINT " The screen will now shift from the alpha mode to the graphics mo 
de and draw" 

3530 PRINT "an empty board with the imaginary regions marked off with dotted 1 in 
es. This" 

3540 PRINT "picture will remain on the screen for only 3 sec after it is drawn b 
ut i t may" 

3550 PRINT "be returned to the screen using k2 to allow further time for check 
ing. " 

3550 PRINT "K1 will return the screen to the alpha mode for the ne.,t- section." 
3570 IF P=0 THEN PRINT 

3530 PRINT " The special function keys will now be used to input the contents 
of each of" 

3590 PRINT "the regions defined on the blank board*. This is the last chance to 
use the " 

360O PRINT "Back-up option <k0> if the board is not correct. A new Key code ove 
r 1 ay will" 

3610 PRINT "be provided if necessary and must replace the general purpose overt a 
y on the" 

3620 PRINT "keyboard. The keys defined in this section act *3 priority interrup 
t s and will" 

3630 PRINT "allow the screen to remain in the graphics mode as the components ar 
e def 1 ned. " 

3640 PRINT "A flashing cross will move to each region in sue ess ion and wait for 
a valid" 
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3650 PRINT "key to be pressed. Every region must be defined even if empty. It 
i s " 

3660 PRINT “important, to wait for the cursor to move to the next empty region be 
fore keys" 

3670 PRINT "are pressed. The bad -up option will NOT work while defining comp 
onents but" 

3680 PRINT "there will be an opportunity to make corrections in the next section 

3680 PRINT "When all regions have been defined, the screen will return to the al 
pha mode to" 

37O0 PRINT "ask if there are any corrections. Fit this point k2 will return the 
graph i cs" 

3710 PRINT "to the screen and kl will return the alpha mode to the screen." 

3720 B 1 ank s = 0 

3730 IF P= 1 6 ThEN G08UB Pagecrt 
3740 PRINT 

3750 PRINT "To make corrections, the numbers of ALL incorrect regions are ente 
red and the" 

3760 PRINT "components erased. When there are no more to be removed, pressing 
CONT with" 

3770 PRINT "no entry will allow all empty regions to be redefined. 8 flashing c 
ros s will" 

3780 PRINT "visit each empty region followed by further correct, ion opportunity a 
s needed. " 

3790 PRINT "When the board is correct, the general purpose key definitions are 
ret urned" 

38O0 PRINT " ar.d the normal key code overlay should be replaced. " 

3810 PRINT 
3320 Pagenum=2 
3830 B 1 ank s = 1 1 

3840 IF P=0 ThEN GOSUB Pageprt 

3850 PRINT " DIP components generally are produced in two different case sty! 
es, ceramic" 

3860 PRINT "or plastic. These are different sizes and this affects the thermal 
analysis" 

3370 PRINT "performed in T H E P H L . The next section provides several methods to 

specify the" 

38S0 PRINT "case styles for the DIP'S. These are: 1 > all pi as tic... 2) all ceram 

i c . . . " 

3890 PRINT "3) majority ceramic (user specify which are p 1 as t i 4 ) majority c 

er antic (user" 

3900 PRINT "specify which are c ►: r am i c ) . . . 5 ) user specify case style for each DIF- 
component . " 

3910 PRINT "In general, opt i on 3 or 4 require the user to input the number of re 
gions that" 

3920 PRINT "are either plastic or ceramic and are useful for specifying the styl 
e when" 

3930 PRINT "when only a few are different. When there is a nearly equal mix, op 
t i on 5 " 

3940 PRINT “requires a 1 for plastic style and a 0 for ceramic style. After all 
DIP'S are" 

3950 PRINT "specified, a list of region numbers is displayed on the screen (in 
alpha mode)" 

3960 PRINT "with the plastic cases in inverse video (numbers black on light back 
round) and" 

3970 PRINT "ceramic cases in normal video. Correct as needed and again no corre 
ct i ons i s " 

3980 PRINT "indicated by pressing CONT Key with no entry." 
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3990 IF F =0 THEN 4040 

4OO0 DISP "PRESS CONT WHEN READY FOR ANOTHER PAGE OF INSTRUCTIONS" 

4010 PRINT LIN(2> 

4020 PAUSE 
4030 GOTO 4O50 
4040 PRINT 

4050 PRINT " The next section concerns the thermally conductive paths that ar 
e added “ 

4030 PRINT "to the circuit board. Both the electrical conductors and any t. herrna 
1 c onduc t ion" 

4070 PRINT "paths or rails will transfer heat between regions. When any conduct, 
i on " 

4080 PRINT "rails are present, the effects of the electrical leads will be negli 
gable and" 

4030 PRINT "ignored. To aid in t. he modeling of the conduction rails it is neces- 
sary to" 

4100 PRINT "introduce the concept of a Conduction Path Unit (.CPU). Each of the 
rails is" 

4110 PRINT "modeled as containing an integer number CPU's. The width of a CPU i 
s defined " 

4 120 PRINT "to t,* .1 mm while the user must define t. he thickness (mm) and therm a 

1 " 

4130 PRINT "conductivity <Watts M-C) of the material that forms the conduction r 
ails." 

4140 PRINT "The thermal c onduc t i v i t i es of several aluminum alloys and copper are- 
presented" 

4150 PRINT "to aid the user. If there are no conduction rails present, the elec 
trical lead" 

4130 PRINT "description is entered in the much the same way; thickness (ram), ave 
rage width" 

4170 PRINT "(mm) and the conductivity ( Wat t s^M-C) . This will default to 334 Wat 
ts/M-C, the" 

4 139 PRINT "value for copper. This is the most coramom material for traces on cir- 
cuit boards. " 

4190 IF P=0 THEN 4240 
4200 PRINT LIN ( 2 > 

4210 DISP "pc-E 38 CONT WHEN READY FOP ANOTHER PAGE OF INSTRUCTIONS" 

4220 PAUSE 
4230 GOTO 4250 
4240 PRINT 

4250 PRINT " In the next section the numbers of leads or CPU's that cross the 
i nt er i or" 

4230 PRINT "boundaries between regions are entered. The display will alternate 
between" 

4270 PRINT "the graphics and the alpha mode with the flashing cross or cursor mo 
vir.g to the" 

4230 PRINT "appropriate location on the board to be defined. The graphics will 
remain on" 

4290 PRINT "the screen for two seconds and then return to the alpha mode for an 
i nput 1 i ne . " 

4300 PRINT "It will ask for the input of the number of leads or CPU's between tw 
o spec i f i c " 

4310 PRINT "regions; ie if the conduction rail is 3.7 mm wide, there are 37 CPU' 
s in that" 

4320 PRINT "path. Special function key 3 <k3) may be used to dump the graphics 
to the" 

4330 PRINT "printer for a hard copy of the graphics picture. It is important to 
account " 
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4340 PRINT " 
number. " 

4350 PRINT - 
ogram flow" 
4360 PRINT •' 
attempt to" 
4370 PRINT - 
ries have" 
4380 PRINT " 
4390 PPINT 
4400 PRINT " 
ransf ered " 
4410 PRINT " 
edges of" 
4420 PRINT - 
edges" 

4430 PRINT - 
bound ary . " 
4440 PRINT " 
the resul ts 
4450 PRINT " 
cuit boards" 
4460 PPINT " 
dissipated" 
4470 PRINT " 
4480 PRINT 
4490 Page rum 
4500 Blanks* 
4510 IF P=0 
4520 PRINT " 
the alpha" 
4530 PRINT " 
ns keys" 

4540 PRINT " 
y. Enter" 
4550 PRINT " 
the numbers 
4580 PRINT " 
of further" 
4570 PRINT " 
omp 1 et e . " 
4580 PRINT " 
tts) or the" 
4590 PRINT " 
4600 PRINT 
4610 PRINT " 
ec i f ied, " 
4620 PRINT " 
ature is" 
4630 PRINT " 
coriiponent " 
4640 PRINT " 

4650 PRINT " 
t i on" 

4660 PRINT " 
mponent " 

4670 PRINT " 
s the link" 



for the leads on both sides of the circuit, board when entering this 
As before follow all entries with the C0NT ke-. Because of the pr 
in this section, the Back_up option <k0> will not function. Do not 
make corrections or return to a previous section until after all ent 
been made. Rll corrections will be made in the next section." 

In the thermal model used by the THEPNELEX system, all heat is t 
to the cooling air stream. None of the heat is. transferred out the 
the circuit board. For this reason, ♦ he sides of the regions at the 
will be labeled with 0 to indicate that no leads or CPU's cross this 
When this heat loss from the board edges and connector is neglected, 
will be conservative. These assumptions mav be unacceptable for c in- 
used in installations where a significant fraction of the total heat 
is conducted awav from the board through card guides or connectors." 

THEN COSUB Pageprt 

When all regions have been specified, the screen will return to 
mode and ask if an 1 ’ corrections are needed. Use the special functio 
as required to either return to gr aph i c s' al ph a or produce a hard cop 
the numbers of the connsc*. ing’ regions that have incorrect values for 
of leads or CPU's. Rs before, C0NT with no entry signals the lack 
corrections. The physical description of the circuit board is now c 
Rll that remains is to specify either the component power levels <Wa 
junction temperatures for each component." 

When the maximum junction temperatures for the component s are sp 
the maximum steady state power level that will result in that temper 
calculated for each component. When the steady state power of each 
is specified, the steady state junction temperatures are calculated. 
In either situation it is also necessary to specify the case to junc 
thermal resistance Rj_c <Peg C'WattJ. The thermal model uses the co 
surface temperatures for heat transfer calculations and Pj_c provide 
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4680 PRINT "betus'sn rhs component power, junction t eriiper at ure and the surface t. e 
mper at ure . " 

4690 PRINT “If R i ■: 13 specified as zero <0> then the case surface t sniper at. ure w 
ill be - 

4?O0 PRINT "equal to the junction t eiriper at ure . " 

4710 PRINT 

4720 PRINT " After all component power levels or junction t emps rat ures have b 
e e n set," 

4730 PRINT "the complete list is displayed in the alpha mode to allow checking a 
nd" 

4740 PRINT "changes. Corrections are inserted by entering the incorrect region 
number then" 

4750 PRINT "entering the correct values in responce to the question asked. If d 
es i red " 

4760 PRINT "a hard copy of all data for the circuit board will be provided from 
the printer." 

4770 PRINT 

4780 PRINT " The next section allows recording of the circuit board descript i 
ve data " 

4730 PRINT "on to a mass storage media. This should be done. The recorded data 
file may be" 

4800 PRINT "accessed at some later time either in BOARD? for editing and chang 
es or" 

4310 PRINT "directly in THEFNL to perform the thermal analysis. BOARD? will 
create the " 

4320 PRINT "the necessary data file to record the description but if a data file 
of " 

4830 PRINT " sufficient size exists, it may be used. All old data in that file w 
ill" 

4840 PRINT "be permanently lost. Information concerning the minimum data file s 
i z e is" 

4850 PRINT "presented by the proqram to help in this decision. " 

4360 PRINT 

4370 PRINT " Built in error traps throughout the entire THERMELEX system will 
save the " 

4330 PRINT "user from loss of previous inputs in case of errors. I nf ormat i on ab 
out the" 

4390 PRINT "error is presented and using the error message guide on the plastic 
slide outs" 

4900 PRINT "below the screen, the user may be able to make corrections. When re- 
ady, as" 

4910 PRINT "signaled by CONT key, program flow returns to a previous point and 
oper at i on" 

4920 PRINT "resumes. The data recording section where file names are entered is 
p ar t i c u 1 ar 1 y " 

4930 PRINT "prone to spelling errors or entry of incorrect mass storage unit spe 
c i f i ers . " 

4940 PRINT "Be careful but remember th system will attempt to catch those inevit 
able errors" 

4950 PRINT "and act. on them before they become fatal." 

4960 PRINT 

4970 PRINT " The end of a cycle through BOftPD ? is now complete. As the page 
he adi ng" 

4930 PRINT "that will appear states, WHAT NOW ?. More work with BOARD? is pos 
3 i b 1 e . The " 

4990 PRINT "user may either input another circuit board description from the key 
board and" 

50O0 PRINT "record that to mass storage-; or, he may retrieve a different circuit 
board" 
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5010 PRINT "discript iv§ data set. from mass storage for checking and corrections. 

5020 PRINT "Another option is to perform the thermal analysis of the circuit boa 
rd using" 

5030 PRINT "the third program in the THERMELEX system, THERML. " 

5040 PRINT 
5050 Pagenum=4 
5060 E 1 an k s = 6 

507O IF P= 0 THEN G0SUB Pageprt 

5030 PRINT " Before an e*p!anation of t-he program flow in THERML . some e •: p 1 
anat ion of" 

5090 PRINT "the use of BOARDS ♦ o edit previously recorded data files is in ord 
er . This" 

5100 PRINT "option may be accesed through the use of RUTOST as explained on t. h 
e first page" 

5110 PRINT "o*' these instructions or the B 0 h R D S program may be directly enters 
d into the" 

5120 PRINT "computer with the command*. LORD " ; CHFf < 34 • : " BOARDS I T 1 5 , 1 " ; CHRS 3 

4>J" then Press EXECUTE " 

5130 PRINT 

5140 PRINT "Cf course the appropriate mass storage specifier should be used both 
here and" 

5150 PRINT "when entering the file name for any previously recored data file. S 
ee the " 

5160 PRINT "Operating and Programming Manual for the 9S45 if the above is not cl 
ear. " 

5170 PRINT 

5 ISO PRINT " Following the input of the file name that contains the descript iv 
e data, " 

5190 PRINT "*;-,e data file will be read and all the variables that are normally s 
ct from" 

5200 PRINT "the keyboard will be defined. The circuit beard will be drawn on th 
e screen" 

5210 PRINT "in graphics mode and changes to the specific components that occupy 
the regions" 

5220 PRINT "may be made. However; it. 
gions in" 

5230 PRINT "any way. This major revis 
rd i nput . " 

5240 PRINT "Correction to case styles is next and then the numbers of leads or C 
PITs " 

5250 PRINT "between regions will be labeled on to the graphics picture of the ci 
rcui t board" 

5260 PRINT "Changes will be allowed to the numbers but no» the physical sizes of 
the leads" 

5270 PRINT "or CPU's; this also requires a new board input. A list of power le 
ve 1 s , " 

5230 PRINT "junction temperature and case to junct'ion thermal resistances is pre 
sent ed " 

5290 PRINT "for checking and corrections as required. This edited descriptive d 
at a set may" 

5300 PRINT "then be recorded either back to the same data file or to a new file. 
5310 PRINT 

5320 PRINT " There are two de mors t r at i on data files included in the THERMELEX 
system. " 

5330 PRINT " D £ M 0 - T contains a typical circuit board description needed to dete 
miine the" 



s not possible to change the number of re 
on of the board must be done as a new boa 
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5340 PRINT "steady state junction temperature ^ i e component power levels are spe 
c i f i e d ) . " 

5350 PRINT “ PE H 0 - P contains a different circuit board for which maximum powers 
levels are “ 

5360 PRINT "determined in the thermal analysis. The use of either of these can 
be helpful" 

5370 PRINT “to familiarize the user with both 3 0 A P B 3 and the last program in t 
he THERELEX" 

5330 PRINT "system, THERNL . " 

5390 PRINT 

5400 PRINT " T H E R M L creates the mathematical model, performs the analysis of 

the circuit" 

5410 PRINT "board and produces the output. The variables used in the thermal mo 
del are set" 

5420 PRINT "using either a data file from mass storage or they are passed in a c 
ommon block" 

5430 PRINT "when T HEP ML is loaded from E 0 A R D 3 . There are no correction oppo 
rtunities" 

5440 PRINT "for the circuit board descriptive' variables in this program." 

5450 PRINT 

5460 PRINT " Assuming that THEPML has been loaded from the What Now option 
list in " 

5470 PRINT 11 E 0 ARDS , the first input is the inlet t empe r at ure < deg C"- of the coo 
ling air." 

5430 PRINT " After this is entered and printed on the screen, the volumetric air 
flow- 

5490 PRINT "ra f e <M A 3^jec) is requested. Since the same set of general purpose 
spec i al " 

5500 PRINT "function keys are valid in this program, the k eye ode- overlay should 
remain and" 

5510 PRINT " k 0 i Back-up' 1 may be used to return to a previous question when input 
errors are" 

5520 PRINT "made. The ne.;t input is the spacing between » he circuit boards for 
rac k mount s . " 

5536 PRINT This spacing and the board height determines the sire of the immagin 
ary air duct." 

5540 PRINT "containing the circuit board. When combined -nth the amount of air 
flow for" 

5550 PRINT "each circuit board the air velocity and heat transfer correlations a 
re" 

5560 PRINT "determined. Since all the power dissipated on the circuit board mus 
t be " 

5570 PRINT "removed by the air flow, it is important for these values to be as p 
rec i se " 

5530 PRINT "as possible. After the thermal analysis is complete, there will be 
opport un i t es " 

5590 PRINT "to alter these important air flow par amt er s and examine the effects 
on the" 

5600 PRINT "output." 

5610 PRINT 
5620 Pagenum=5 
5630 Blank s = 6 

5640 IF P=0 THEN COSUB Pageprt 

5650 PRINT " Eefore the analysis begins, the convergence criteria must be set 
The" 

5660 PRINT "analysis consists of assuming a solution then defining a set of sim-j 
t aneous" 

5670 PPINT "equations (one for each region) which are solved for either the powe 
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5630 PRIMT "or the junction t empe rat ures . The result* are compared to the previ 
ouj results" 

5690 PRINT "(the assumed values for the first ittjration) and if the largest dif 
f erenc e " 

5700 PRINT "between these values is less than a maximum set by the user, the res 
ults are" 

5710 PRINT "presented in graphics. If the convergence criteria is not met, the 
resul ts" 

5720 PRINT "are printed to the screen in alpha mode and another i iteration is en 
tered. “ 

5730 PRINT "The closer to zero one sets the convergence criteria, the longer the 
proc ess" 

5740 PRINT "takes. Default values for the convergence criteria are .5 Deg C or 
1 5c change in" 

5750 PRINT "power. These typically require 2 to 6 i tt.gr at ions before they are m 
et. Each" 

5760 PRINT "itteration ta*es 5 to 100 sec to perform depending on the number of 
reg i ons. " 

5770 PRINT "fi: the solution precedes the user may view the i nt e rmed i at e results 
to watch" 

5730 PRINT "the progress." 

5790 PRINT 

5300 PRINT " The final results are written onto a picture of the circuit, boar 
d with" 

5310 PRINT "each region containing component type, junction temperature, power 1 
evel and" 

5320 PRINT "case temperature. Empty, regions contain only the temperature of the 
circuit" 

5330 PRINT "board. Those junction temper attires within 5*. of the maximum and t ho 
se power" 

5340 PRINT "levels within 5*1 of the minimum are starred (•**•> for easy reference 
as trouble" 

5350 PFINT "spots. This graphics output is automatically dumped to the printer 
to insure" 

5860 PRINT "that a hard copy of the results exist." 

5370 PRINT 

5880 PRINT " This completes a cycle through THERNL and again the question i 
s What Now 

5890 PRINT " E C A F D 3 may be returned to core, to make changes to the circuit boa 
rd or" 

59O0 PRINT "to enter a new circuit board. Another analysis of the same circuit 
board" 

5910 PRINT "may be done with a new set of air flow paramters or, a new circuit b 
oard" 

5920 PRINT "descriptive data set may be read in from mass storage with THEPML . 
5938 PRINT 

5940 PRINT " fin additional option is sensitivity analysis. Plots of the ma-.i 
mum " 

5950 PPINT " Tjunc vs. Air Flow Rate or the minimum Power vs. Air Flow Rate may b 
e pr oduc e d . “ 

5960 PRINT "A maximum air flow rate is specified and five separate analyses are 
performed" 

5970 PPINT "and the results plotted and dumped to the printer. The maximum air 
flow rate" 

5930 PRINT "specified should be an integer multiple of five to make for better 1 
ook i ng axes . " 

5990 PRINT "It is also possible to produce hard copys of the printed results for 
each of the" 
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6000 PRINT "air flow rates. Upon completion of this sensitivity analysis, progr 
am flow" 

6010 PRINT "returns to the option list previously discussed." 

6010 PRINT 

6030 PRINT " Termination should be done under program control using the final 
opt ion." 

6040 PRINT "This will insure that the standard key definitions are returned arid 
that the" 

6050 PRINT "proper graphics parameters are set for the next user of the computer 

6060 PRINT LIH<2>,Tft£<32), "GOOD LUCK" 

607C Pagenum=6 
6030 B 1 an k s = 1 6 

6090 IF P=0 THEN GO SUB Paoeprt 
6100 PRINTER 13 16 
6110 RETURN 
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5 

10 

15 

20 

25 
30 
35 
40 
45 
50 
55 
60 
65 
70 
75 
80 
85 
90 
95 
100 
105 
1 10 
115 
120 
125 
130 
135 

140 St ar t_c- s r : OH Hap GOTO 01 dpi ct_correct , Bdpi ct. , 01 dpi ct 
145 DISP“"MAP= " {Map 
150 PAUSE 

155 Starts: COSUB Error 

160 Start: Sf= "THERMAL ANALYSIS PROGRAM FOR ELECTRONIC CIRCUIT BOARDS" 

165 ON ERROR GOTO Starts 
170 GOSUE Page he ad 

175 PRINT " Thii progr am allots the user t input a detailed description of an 
e 1 ec tron i c " 

ISO PRINT "circuit board and perform a thermal analysii to predict either the s 
teady state" 

185 PRINT "temperatures of the components if power levels are given or the comp 
onent power" 

190 PRINT "level that will result in the steady state temperature specified by 
the user . " 

195 PRINT "The description can come from a data file previously created by this 
progr am" 

200 PRINT "or can be entered at this time." 

205 A n s i = " N " 

210 INPUT "DO YOU DESIRE TO RETRIEVE A PREVIOUSLY STORED DESCRIPTION <N OR Y>" , 
fins* 

215 IF Ans -f = " E A C K_UP " THEN 160 

220 IF < UPC f ''flusfll, 1 3 .■ * " N " ) OP < UPCf C Rnsf Cl, 1 3 > = " Y" > THEN 235 
225 GOSUE Err in 
230 GOTO St art 

235 IF UPC f < Ans f C l , 1 ] > = " N " THEN Edpict 

240 GOTO Oldpict 

245 Oldpicte: GOSUE Error 

250 Oldpict: 1 THIS SECTION RETRIEVES A BOARD DESCRIPTION OFF A MASS STORAGE DEV 

ICE 

255 ON ERROR GOTO Oldpicte 



BOARDS 






DATA INPUT AND STORAGE FROGRAM FOR 

THERMELEX ***< 

A SYSTEM OF PROGRAMS FOR THE HP 9345 
TO PERFORM THERMAL ANALYSIS OF ELECTRONIC CIRCUIT BOARDS 

PREPARED AT THE NAVAL POSTGRADUATE SCHOOL MONTEREY CA 
R. A. FOLTZ LCDR USN 



Mill 



BOARDS 

* 4 -*r Hr 



Mill 



OPTION BASE 1 
PRINTER IS 16 

COM Map, Name i. Pi ct 1 bl f C503 . Case ft 50I , Bdl , Edh, Sea, Thi ck_b, Tem_*o 1 .Kb, A1 , K1 
COM SHORT Xt50 • . Yv50,' , Tj<5CO , Pcwt 50 > , P j_c <50> 

COM INTEGER I t vpe < 50 > , N 1 < 4 . 50 > , Nxr , Nyr , Nreg 
SHORT L e < 3 0 > , U e • 5 0 > 

INTEGER Np inf 59' 

DIM Sf C703 . Msu=lC223 
IF Map = 0 THEN St art 
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260 0 1 de =P I 

265 Si=" BOARD DESCRIPTION FROM MASS STORAGE DEVICE” 

2?0 GOSUB Page he ad 

275 PRINT " The data file containing the board description must have been st 
ored by " 

230 PRINT "this program. Enter the data file name below, be sure to include th 
e m as s " 

285 PRINT "storage unit, specifier if needed, eg : T14, : F3, : Y12 etc" 

290 EDIT "UNDER WHAT FILE NAME IS THE CIRCUIT BOARD FILED (input or change belo 
u) ", Name * 

29 5 IF N aide t~" BACK UP" THEN 160 
300 IF Name** "NULL" THEN Oldpict 
305 ASSIGN #5 TO Name* 

310 READ 45, 1 

315 READ *5 ; N ame T , P i ■: t 1 b 1 T , C as e * , Bd 1 , Bdh , Sc a, N;,r , Nyr , Mr eg , Th i c k_b , Kb , A 1 , Y. 1 , T em_s 
o 1 

320 FOR 1 = 1 TO Nr eg 

325 READ #5 ; 1 1. ype I > , I > , Y< I) , T j < I > , Row! I ) , R j_c < I > 

330 FOR J*1 TO 4 
335 READ 45; Nl f J, I > 

340 NEXT J 
345 NEXT I 

350 Oldpict correct: GOSUB Board_pict 

355 GRAPHICS 

360 FOR 1=1 TO Nr eg 

365 IF 1=1 THEN 385 

370 IF KNxr-l THEN 330 

375 IF < 1 1 y p e < I ) = 1 0 > AND It ype • I -N \r >= 1 0 ' THEN 400 

380 IF < I type (I =9) AND ■ I type I -1 >=9J THEN 400 

335 IF 1 1 y p e I ) = 0 THEN GOSUB Empty 

390 IF 1 1 v p e < I ) = 0 THEN 400 

395 ON I type (I) GOSUE Hd i p 1 4 , Vd i p 14 , Hd i p 16, 7d i p 1 € , Hd i p24 , Vd i p24 , Errol de , Erre 1 de 
, Hd i p40 , Vd i p40 . F 1 at 1 6 , F 1 at 24 , F 1 a* 40 . F 1 at 64 



400 


NEXT I 








405 


GOSUB Bd cha 








410 


IF Ans* = " BACK 


JJP" 


THEN 


01 dp 


415 


GOSUB Case_t.j 


P 






420 


IF AnsT = 'BACK. 


UP " 


THEN 


405 


425 


GOSUB Lsads_o 


Id 






430 


IF fins *=" BACK 


UP" 


THEN 


415 


435 


GOSUB Le ads_c 


ha 






440 


IF Ans*= " BACK 


JJP" 


THEN 


425 


445 


GOSUB Tempi n 








450 


IF A n s S = " B A C K 


JJP" 


THEN 


435 


455 


GOSUB Stow 








460 


IF fins-T = " BACK 


UP" 


THEN 


445 


465 


GOSUB What now 






470 


IF Ansf = " BACK 


_UP" 


THEN 


455 



475 GOTO St art _oue r 
480 Bdpicte.* GOSUB Error 

485 Edpi ct : ! THIS SECTION IS THE NORMAL ENTRY POINT FOR NEW BOARD INPUT 
490 ON ERROR GOTO Bdpicte 
495 0 1 de=0 

500 $*=" KEYBOARD INPUT OF CIRCUIT BOARD DESCRIPTION" 

505 GOSUB Pagehsad 

510 PRINT " This routine draws a blank circuit board and divides is into regi 
ons of " 

515 PRINT "interest as determined by the user. There are a maximum of 50 regions 
auai able " 
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520 PRINT "with on* component o'regi on. The circuit board 1 a then presented on th 
e " 

525 PRINT "CRT either normal =i:e or scaled up or down by some integer it' larger 
or smaller" 

530 FRINT "than 120mm by 170mm. The specific components are selected with the us 
er keys" 

535 PRINT “in the upper right corner of the keyboard. If you do not have a key 
code " 

540 PRINT "overlay for the keys the program will provide one for you. Informat i 
on to be" 

545 PRINT "used in the thermal analysis routine is entered in this section of th 
e progr am . " 

550 PRINT “You will ha--e the option of recording the board description onto mass 
st or age . " 

555 PRiNT 

560 PRINT "Conductivity of the board will default to .2942 W at t s/M-de gK with no 
entry. " 

565 PRINT LINC 1 > , "NOTE : LENGTH IS DEFINED TO EE IN THE DIRECTION OF AIR FLOW 

570 MAT 1 1 ope =ZER 
575 Ansf =" NULL " 

5 SO INPUT "WHAT IS THE LENGTH OF THE BOARD IN MILLIMETERS" , AnsS 
535 IF Ans$ = 1 BACK UP" THEN Start 
590 IF AnsX= " NULL " THEN 575 
595 Bd 1 =AE$C VAL < Ans S .- ’> 

600 PRINT LINC 1 >, TAB-. 10) , "LENGTH =" ;BdlJ" mm" 

605 Seal ex*INTCBdl 171.- + 1 

610 IF B d 1 < 1 7 1 x 2 THEN Se al ex = .5 

61? Ansi= " NULL " 

620 INPUT "WHAT IS THE HEIGHT OF THE BOAPD IN MILLIMETERS" , Ans $ 

625 IF Ans t = "BACK_UP" THEN Bdpict 
630 IF Ansf =" NULL" THEN 615 
6 35 Bdh = AES< VAL< Ans $ :- > 

640 PRINT LINC 1 TAE-- 10 -, "HEIGHT *" ;Bdh;" mm" 

645 Sc al ey= I NT < Bdhx 1 2 1 > «■ 1 

650 IF Bdh< 121 2 THEN Scaley=.5 

655 Sca=MAX<Scal*y. Scale,,) 1 SCALE FACTOR TO FILL CRT DISPLAY 
660 IF Sc a= 1 THEN 630 

665 IF Sc a< 1 THEN PRINT L I N < 2 ) , “ *4**N0TE***e DUE TO SMALL SIZE OF THE BOA 

RD CRT DISPLAY WILL BE 2X SIZE" 

670 IF S c a ' 1 THEN PRINT L I N < 2 > , " **»*NQTE* *** DUE TO LA r ‘GE SIZE OF THE BOA 

RD CRT DISPLAY WILL BE 1 " " ; Sc a; " S I ZE " 

675 An s S = " NULL " 

630 INPUT "WHAT IS THE NUMBER OF DIVISIONS --REGIONS) IN THE HORIZONTAL DIPECTIO 
N " , Flnsl 

6S5 IF A r. s T = " B A C K U P " THEN 615 
690 IF Ans i = " NULL"" THEN 675 
695 Nxr* ABS ( VAL < Ans t ' ) 

700 PRINT L INC 1 ) , TABc 1 S ) , " # Xreg =";Hxr 
705 Ans$=" NULL" 

710 INPUT "WHAT IS THE NUMBER OF DIVISIONS (REGIONS) IN THE VERTICAL DIRECTION" 

, Ansi 

715 IF Ansf = " BACK UP" THEN 675 
720 IF Ans T = “ NULL " THEN 705 
725 Nyr = ABS ( VAL < Ans * > ) 

730 PRINT LINC 1 ), TAB' 10) . Yreg =" ;Nyr 

735 Nreg=Nxr*Nyr ! Nreg = HUMBER OF REGIONS 

740 IF Nr e g<' 5 1 THEN 770 
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" , An 5 I 

= " ; Thickj; " mm" 

ONDUCTIVITY OF THE BOARD? < WATTS- "M -dag K>" 



745 BEEP 

750 PRINT PAGE, LIN< 15) , SPA<30> , "TOO MANY REGIONS" 

755 WAIT 1000 
760 GOSUB Err in 
765 GOTO B d p i c t 
770 Ans$="NULL" 

775 INPUT "WHAT IS THE THICKNESS OF THE BOARD f..m 
730 IF Ansl=" BACkJJP" THEN 705 
785 IF Ans 1 = " NULL " THEN 770 
7 90 Th i c k_b = RBS < VAL Ans S > > 

795 PRINT LI N(1 >, TAB< 10) . "Thick 
800 Ansi = " . 2942" 

80S EDIT "WHAT IS THE THERMAL C 
810 IF Ans *»" BACKJJP* THEN 755 
815 Kb=AB3 < VAL < Ansi } 

820 PRINT LIN( 1 ) , TABU0J- , "Cond £ *"JKb 

825 EDIT "INPUT A SHORT -.<50 CHARACTERS > WORD DESCRIPTION OF THE CIRCUIT 
HERE ",Pict1bl* 

830 IF Piet 1 bl «*" BACK JP" THEN 800 

835 IF LEN<Pictlbll-'v51 THEN 350 

840 GOSUB Ernn 

345 GOTO S25 

850 GOSUB Board_pict 

855 GOSUB Graphic ey 

860 IF Ansl="EACK UP" THEN Bdpict 

865 GOSUB Bd eh a 

370 IF Ansi = " B AC V JJP " THEN 850 

875 GOSUB base t 1 ■ p 

830 IF Ansi = " EAC K _UP " THEN 365 
985 GOSUB Leads 
890 IF flrtsl= " BACKUP" THEN 375 
895 GOTO 4S5 
900 Eoard_p ict« : 

90S Board_p i c I 



1 END OF Ed p l •: t. 



OSUB Error 

This SECTION PLACES 



A BLANK. BOARD ON THE SCREEN 



910 ON ERROR GOTO Bo ard_pi.it * 

915 PLOTTER IS "GRAPHICS" 

920 GRAPHICS 
925 MSCALE 0, 10 
930 CSIZE 3 
935 LORG 5 
940 MOVE 26,130 

94 5 LABEL " air flow > " 

950 Bdh=Bdh/3c a 

955 Bdl=£dl/Sca 
960 LORG 6 
965 LDIR PI/2 
970 MOVE 175,70 

975 IF Sc a> 1 THEN LABEL "THIS PICTURE IS 1 ✓ " &VAL1 < Sc a> & " SIZE" 

930 IF 3c a< 1 THEN LABEL "THIS PICTURE IS 2X SIZE" 

985 LORG 5 
990 LDIR 0 

995 MOVE 92-LEH C P i c t 1 b 1 1 'i 2 , 135 

10O0 LABEL USING "K" ; R l c »- 1 1 1 1 

1005 MOVE 0,O 

1010 DP AW O , B d h 

1015 DRAW Ed 1 , Bdh 

1020 DRAW B d 1 , u 

1025 DRAW 0,0 
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BOARDS 



DATA INPUT FOR THERMELE 



1030 LIME TYPE 3 

10 35 Lr = Ed1 .'N'T « LR = LENGTH OF EACH REGION 

1040 Hr = BdlvNyr » HP = HEIGHT OF EACH REGION 

1045 FOP 1=1 TO N r - 1 

1050 MOVE I * L r , 0 

1055 DRAW I * l_ r , B d h 

1060 NEXT I 

1065 FOR 1=1 TO Nvr-1 

1070 MOVE 0, I *Hr 

1075 DRAW Bdl , I*Hr 

1030 NEXT I 

1035 LINE TYPE 1 

1090 LORG 5 

1095 Nr =0 

1100 FOR J=Nyr TO 1 STEP -1 

1105 FOR 1 = 1 TO Nxr 

1110 Nr=Nr+l 

1115 X<Mr;=Lr*> 1-1 '*Lr • 2 

1 120 Y < Nr > = < J-l -Hr + a-t-Hr. 2 > 

1125 MOVE X C Nr > , V Nr .> 

1130 LABEL USING " K " ; VAL i 1 Nr > 

1135 NEXT I 

1140 NEXT J 
1145 Bdh=Bdh*3ca 
1150 Bd1=Bd1*Sca 
1155 WAIT 1500 
1160 EXIT GRAPHICS 
1165 RETURN 



1170 Graphkeye: GOSUB Error 

1175 Graphkey: ! HERE IS SUBROUTINE TO ALLOW INPUT OF BOARD DEC F: I FT I ON FROM KEYS 
DATA INPUT 



11 SO ON ERROR GOTO Graphkeve 
1135 AnsX="Y" 

1190 31= "KEY CODE OVERLAY" 

1195 GOSUB Page he ad 

1200 PRINT *’ This section allows the components to be defined and drawn on the 
screen" 

1205 PRINT “using the keys in the upper right corner of the keyboard. The previ 
ous key" 

1210 PRINT "deff ini t ions are not oalid while entering the components but will 
be " 

1215 PRINT "during later portions of this program. As the flashing cursor moves 
to each" 

1229 PRINT “region in numerical order press the coresponding key and wait for t. h 
e cursor to" 

1225 PRINT "appear at the next region. All regions must be defined even if with 
k 0= >EMPT Y . " 

1230 PRINT 

1235 PRINT " NOTE: kll- kl4 refer to chip carriers or flat packs." 

1240 PRINT 

1245 PRINT " If you do not have a plastic kev code overl ay filled in 'HP pa-'t 
# 7 1 20-6 1 64 > , I will make a paper one for you." 



1250 A n s S = " N " 

1255 INPUT "DO YOU NEED A PAPER OVERLAY <N or Y)?",AnsX 
1260 IF AnsT = " EACf _JJP " THEN RETURN 
1265 IF UPCt’.AnsXC 1 , 1 3> = "N" THEN 1320 
1270 PRINTER IS 0 

1275 PRINT SPA< 10> , "TEAR OFF FOR COMPONENT DEFF INITIQN KEY CODE OVERLAY" , LIN<2 
) 
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1230 

1235 

BOVE 

1290 

T 

1295 

1300 

1305 

BELOW 

1310 

1315 

1320 

1325 

1330 

1335 

1340 

1345 

1350 

1355 

1330 

1365 

1370 

1375 

1380 

1335 

13 30 
1395 
1400 
1405 
1410 
1415 
1420 
1425 
1430 

14 35 
1440 
1445 
1450 
1455 
1460 
1465 
1470 
1475 
14S0 
1435 
1490 
1495 
1500 
1505 
1510 
1515 





BOARDS . 


. DATA 


INPUT FOP 


thermele: 






1 EMPTY 
j . (ko.) 


1 HOPIZ 
| DIP14 | 


VERT 
B I P 1 4 


| H0RI2 1 
1 DIP 16 


VERT 1 
DIP 16 j 


HOPIZ 

DIP24 


1 VERT 1 

DIP24 j 



PRINT 

PRINT "|.<ko> | DIP14 | BIP14 j D I P 1 6 | D I P 1 6 | DIP24 | BIP24 | PLACE ft 
KEYS" 

PRINT "TEAR HERE THEN PRESS CON 



1 hopi: 

[ BIP40 



| VERT j FLAT | 
| DIP40 [ 16 PIN | 



FLAT | FLAT | FLAT | 

>4 PIN | 40 PIN | 64 PIN [ " 



PRINT LINO) 

PAUSE 

PRINT "| < k 3 ) j HOPi: I VERT | FLAT I FLAT | FLAT | FLAT J PLACE 
KEYS " 

PRINT "| 

PRINT L I N( 4 > 

PRINTER IS 16 

! NOTE THAT THE KEY# IS THE CODE FOR TYPE OF ELEMENT STORED IN Ityped) 

ON KEY #0 GOTO Emptv 
ON KEY #1 GOTO Hdipl4 
OH KEY #2 GOTO Vd»pl4 
ON KEY #3 GOTO Hdipl© 

ON KEY #4 GOTO Vdiplo 
ON l-EY #5 GOTO Hdip24 
ON KEY #6 GOTO Vdip24 
ON KEY #7 GOTO 1435 
OH KEY #3 GOTO 1435 
ON KEY #9 GOTO Hdip40 
ON KEY #10 GOTO Vdip40 
ON KEY #11 GOTO FI it 16 
ON KEY #12 GOTO FI at 24 
ON KEY #13 GOTO FI at 40 
ON KEY #14 GOTO FI at 64 
IF 0 1 de = F I THEN RETURN 
GRAPHICS 

FOR 1=1 TO Nreg 
IF I type K I > = 1S THEN Ne.tr 
Loop: MOVE X<I),Y<I> 

POINTER M< I ) , YCI ) ,2 

! THIS IS AN ENDLESS LOOP TO ALLOW FOR USE OF THE PRIORITY INTTERUPT KEYS 
GOTO Loop 

Next r : IF Peg cr,a"'0 THEN RETURN 

IF I < >Nreg~ THEN 1475 
IF OldeOPI THEN WAIT lOO + Nregd.S 
POINTER 0,0,2 
EMIT GRAPHICS 
RETURN 
NEYT I 

Empty; GOSUB Erasenwm 
1 1 ype < I ) = Np i n ( I > =We < I • =Le< I > =0 
GOTO Next r 

Hdi p 1 4 : GOSUB Erajenurfi 
Np i n < I > = 1 4 

IF Itype< I )=2 THEN LB IP PI/ ’2 
LABEL USING "K" ; " D IP 1 4" 

We<I>=We=.25^25.4 ! ACTUAL WIDTH IN mm 



1520 Led >=Le=. 725-25. 4 ! ACTUAL LENGTH IN mm 

1525 IF d t y p e d ) = 1 > OP < 1 1 ype < I ) =0 > THEN 1545 

1530 Temp_dim = Wed> ! HERE LENGTH AND WIDTH TRANSPOSED SUCH THAT LENGTH 

1535 We<I J=We=Le<. I > ! IS DEFINED TO BE IN THE DIRECTION OF AIR-FLOW 

1540 Le< I >=L*=Temp dim 
1545 IF I typed >*0~THEN Ityped> = l 

1550 Drawdi p : ! THIS SECTION WILL DRAW ftLL THE DIP FIGURES ON THE CRT IF NpirvdS 



< 6 ) 
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BOARDS .. DATA INPUT FOR THEPMELS',-: 



1555 
1 560 
1565 
1570 
1575 
1530 
1535 
1 5 ?0 
1535 
1600 
1605 
1610 
1615 
1620 
1625 
1630 
1635 
1640 
1645 
1650 
1655 
1660 
1665 
1670 
1675 
1630 
1635 
1690 
1695 
1700 
1705 
1710 
1715 
1720 
1725 
17 30 
1735 
1740 
1745 
1750 
1755 
1760 
1765 
1778 
1775 
1730 
1735 
1790 
1795 
1800 



L*=L*/Sca 



MOVE >{< I > -L i' 2, V ■ I --We 2 
DRAW X< I .• -Le 2 , Y < I 1 + W« 2 
DRAW XCD+Le 2,Y(I'HU 2 
DRAW D+Le ■ 2. YU .'-W* 2 
DRAW XCD-Le 2,YvI>-We 2 
LD I R 0 
GOTO Near 



Vd i p 1 4 : I * ype < I > =2 
GOTO H d i p 1 4 
Hdipl6:G0SUB Eraienum 
Npin(D = 16 

IF I tops I >=4 THEM LD I R PI - 2 
LABEL U5IHG "K"; “ D I P 1 S '* 

Ue< I )=Wa=. 25*25. 4 ! WIDTH IN mm 

Let I >=Le=. 325*25. 4 ! LENGTH IN mm 

IF I type -I >04 THEN It...p*<I> = 3 
IF 1 1 ype ■ I > 4 THEN Drawdip 

T emp_d i m=We CD 
We < I ) =Wf=Lf i I ) 

Le(I>=Ls=Te mp_d i m 
GOTO Drawdip 
Vd i p 16: It ype ■ D =4 
GOTO Hd i p 1 6 
Hdip24:G0SUB Erasenum 
Up i n< I; =24 

IF Itypecl '=6 THEN LDIP PD2 
LABEL USING " K " ; "DIP 24" 

Wet I )=We = . 55*25. 4 



Let I ) =Le= 1 . 25*25. 4 
IF 1 1 y p e 1 1 ) < ; 6 THEN I typed ; = 5 
IF I type CIK,** THEN Drawdip 
T emp_d i r.i = We ■ I 1 
We < D =We = Le t I ) 



! End of Draw_dip 



! End of Hd i p 1 6 



Let I > =L e = Te rr> p_o i m 

GOTO Drawdip ! Ef'd of Hdip24 

Vd i p24 : It ype < l ;• =6 
GOTO H d i p 2 4 

Hdi p40: IF t\ MOD Nxr = c:- OR < I t ype < I + 1 > = 1 O > AND < Re g_c ha=0 > THEN Loop 
G03UB Erase num 

I = I h= I + 1 

GO SUB Erasenura 
1 = 1-1 

Nptnf I)=Npin-:ih.)=20 ! THIS IS FOR EACH HALF OF THE 40 PINS 

I I ype< I > = I type 1 1 h ) =9 
MOVE X( D+Lr/2, Y< D 
LABEL USING " K " : " D I P 4 0 " 

Wet I ■ = We Uh >=We = . 6*25. 4 

Let D=Le < Ih>=Le = 25. 4 ' THIS IS FOP EACH HALF IN EACH REGION 



1305 We=W*'3ca 
1310 L e = L e ' 3 c a 

1815 MOVE X ( I > + L r - 2 - L e- , Y t I > + W e ■' 2 
1320 DRAW XU >+Lr/2-Le, Y<D-We '2 
1325 DRAW X< D+Lr/2 + Le, Y< I >-We 2 
1830 DRAW X < I > L r / 2 + L e , Y U > *■ W e ' 2 
1335 DRAW X< I )+Lr. 2-Le, Y( I)+WfU 
1S40 I=Ih 
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1845 

1 S50 

1855 

I860 

1885 

1S?0 

1 S75 

1SS0 

1335 

1890 

1895 

1900 

1905 

1910 

1915 

1920 

1925 

1930 

1935 

1940 

1945 

1950 

1955 

1 960 

1965 

1970 

1975 

1930 

1935 

1998 

1995 

2000 

2005 

2010 

2015 

2020 

2025 

2030 

2035 

2040 

2045 

2050 

2055 

2060 

2065 

2070 

2075 

2080 

2035 

2090 

2095 

2100 

2105 

2110 

2115 

2120 

2125 

2130 



BOARDS .. DATA INPUT FOP THEPMELE:: 



GOTO Mixtr 

Vdi p40: IF I+N;r'Nr*g THEM 1425 
GOSUB Er a-enufii 
I = I h = I +Nxr 
GOSUB Eraserium 
I = I -N <r 

Np 1 n< I > =Np i n < Ih j =20 ! THI3 IS FOR THE HALF OF THE 40 PINS IN EACH REGION 

ItypeC I )=Ityps< Ih>=10 
LDIR PI '2 

MOVE X < I > , YC I > -Hr/2 
LABEL USING " K ** ; " B I P 4 0 " 

LDIR 0 

Ue < I > = Uf> Ih)=We=. 6^25. 4 

LeC I>=L£CIh.'=L* = 25. 4 ! THIS IS FOR HALF 

We=We- Sc a 
Le*Le/Sc a 

MOVE X(I)tW« 2, Y< I>-Hr^2+L* 

DRAW X< I >-He 2, Y< I >-Hr ' 2+Le 
DRAW XC I )-M«/2,ya)-Hr^-Le 
DRAW XC I > + W e ' 2 , Y I)-Hr^2-Lc 
DRAW X < I > + We 2, YC I >-Hr/2+Le 
LDIR O 

T emp_d i m = Le C I > ! HE 5 E TRANSPOSE Ls AND U* TO MAINTAIN LENGTH DEFF WITH AIR 
Le< I>=Le< Ih>=lU ■: I> 

U?C I >=W* C Ih) = Teri 1 p_di f,-, 

GOTO N*;;ir 
FI at 16: It ype Cl =1 1 
GOSUB Erasenum 
Np l n C I > = 1 6 

LABEL USING M K";"16" 

L«=Ue=W*' I >=L* C I .‘=4.57 
GOTO Dr&wdip 
F 1 at 24 : It yp* C I .' = 1 2 
GOSUB Eraienum 
NpinC I>=24 

LABEL USING ,, K";"24" 

L‘? = W= = We , i:=Le<I)=?.75 
GOTO Dra"dip 
FI at 40: It yp * < I > = 13 
GOSUB Era-inum 
NpinC I>=40 

LABEL USING " K " ; "40" 
l_e=We=W*C I > = L f c I ’ = 12. 19 
GOTO Draudip 
F 1 at64 : I t ype 1 1=14 
GOSUB Eraisnurt 
NpinC I >=64 

LABEL USING " K" ; " FLAT64 " 

L**We=We< I >=Ls < I > = 13. 29 

GOTO Drawdip ‘ END OF COMPONENT DEFF INI T I ONS 

! 

Er as *num : ! HERE WE MOVE THE REGION NUMBER TO THE UPPER LEFT CORNER 
CSIZE 3 

MOVE X < I > , Y < I ) 

PEN -1 

LABEL USING "K";l 
CSIZE 2.5^Sca 
IF R*’5_cha=0 THEN PEN 1 
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BOARDS PniM INPUT FOR t n e k m e l e . . 



2135 MOVE X<I J-.4*Lr.Y<I>-t.4*Hr 
2140 LABEL USING "K" ; I 
2145 MOVE X(P.YU) 

2150 CSI2E 3. 4/Sc a 

2155 RETURN ! END OF ERflSENUM 

2160 ! 

2165 Bd_chaf : GOSUB Error 

2170 Bd_cha: « THIS ROUTINE ALLOWS FOR CHANGE OF TYPE FOR SPECIFIED COMPONENTS 
2175 ON ERROR GOTO £d_ch*e 
2 1 SO S -I = " C H A N G E COMPONENTS" 

2185 GOSUB Pagehead 

2190 PRINT " This section allows correct i ons to the components on the circuit 
board. " 

2185 PRINT "However; you mav not change the circuit board itself. You will firs 
t remove" 

2208 PRINT "all unwanted components (up to 1 O ) by entering the region number on 
the picture." 

2205 PRINT "When removing any component which requires two spaces, use the lower 
region num . " 

2210 PRINT "Entering a 0 or pressing CONT with no ertm will allow program to 
continue" 

2215 PRINT "with the section that allows deffimtion of components in all empty 
spaces. " 

2220 PRINT 

2225 PRINT "Remember to allow space for the larger components that require two r 
eg i ons . " 

2230 ltt=0 

2235 Again: Itt=Itt+l 
2240 flnsIs'O" 

2245 INPUT "ENTER THE REGION NUMBER TO REMOVE COMPONENT < DEFAULT = ) 0 = ) NONE ) . " , An 

si 

2250 IF Ansf="£ACKJJP" T HEN RETURN 
2255 Reg_cha=VAL<Ans *) 

2260 IF Reg_ch a=0 THEN Add_comp 

2265 IF <Reg_cha>0> AND < Reg_c ha-. Nreg+ 1 > THEN 2238 
2270 GOSUB Err i r, 

2275 GOTO Again 

2280 Reg_cha f : 1 1 1 >=Req_cha 

22S5 GOSUB Erase_comp 

2290 Reg_cha=0 

2295 IF 1 1 1 < = 1 6 THEN Again 

2300 S$ = " MAX I MUM NUMBER OF REGIONS" 

2305 GOSUB Page he ad 

2310 PRINT " There are 10 changes to accomplish at this point, and to avoid pro 
blems you" 

2315 PRINT "must redefine the board." 

2320 PRINT LINO) , SPA<15> , "Press CONT " 

2325 GOSUB Add_comp 
2330 GOTO Ed_c h a 

2335 Erase comp!! THIS SECTION REMOVE UNWANTED COMPONENT TO ALLOW CHANGE 
2340 GRAPHICS 
2345 PEN -1 
2350 I=Peg_cha 

2355 IF I type <Reg_cha)=0 THEN 2365 

2360 ON ItypeO) GOSUB Hd i p 1 4 . Vd i p 1 4 . Hd i p 1 6 . Vd i p 1 6 , Hd i p24 , Vd i p24 , Er ro 1 de , Err o 1 de 

, Hd i p40 , Vd i p48 , F 1 at 1 6, F 1 at. 24, FI at40,Fl at 64 

2365 PEN 1 

2370 GOSUB Empty 



< 9 > 
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BOARDS 



data input for thermelex 



2375 IF I < 2 THEM 24*0 

23S'0 IF < It ype< 1-1 >=9 > AMD Cl t ype C I >=9 ) THEM 2390 

23S5 GOTO 2400 

2350 1=1-1 

2355 GO SUB Empty 

2400 IF Kllxr+1 THEM 2420 

2405 IF I type < I -H;>r ><>10 THEM 2420 

2410 I = I -Mxr 

2415 GOSUB Empty 

2420 WAIT 500 

2425 EXIT GRAPHICS 

2430 RETURN 

2435 Add comp: i THIS SECTION ALLOWS ADDITION OF COMPONENTS TO EMPTY SPACES. 

2440 IF Ttt<2 THEN 2475 

2445 IF 01de=PI THEN GOSUB Graphkey 

2450 FOR 1=1 TO Nr eg 

2455 IF ItypeUX'O THEN 2470 

2460 POINTER XCI>,YcI>,2 

2465 GOTO 2460 

2470 NEXT I 

2475 FOR 1=1 TO 15 

2480 OFF KEY #1 

2435 NEXT I 

2450 EXIT GRAPHICS 

2495 RETURN ! END OF ADD COMP 

2500 ! 

2505 Case — type : GOSUB Error 

2510 Case~typ: ' THIS SECTION ALLOWS THE USER TO SPECIFY THE CASE TYPE FOR DIPs 
2515 ON ERROR GOTO Case type 

2520 IF < 0 1 de = P I > AND <I»t<2> THEN GOTO Case__cna 

2525 ! BUT FIRST MUST TELL THE USER TO CHANGE THE KEY-CODE OVERLAY 
2530 31= " CHANGE OF THE SPECIAL FUNCTION KEYS" 

2535 GOSUB Pagehead 

2540 PRINT LINU0,',SPA< lO' , "REMOVE THE KEY CODE OVERLAY FOP COMPONENT INPUT" 

2545 PRINT L I N < 1 ; , SPA 1 0 , " REPLACE THE GENERAL PURPOSE KEY CODE OVERLAY" 

2550 DI3P SPA C 25 , "PRESS CONT WHEN READY" 

2555 PAUSE 

25S0 S f = " T Y P E OF DIP PACKAGE" 

2565 GOSUB Page head 

2570 PRINT " Dual inline packages CDIPs) generally come in either an inject io 
n molded" 

2575 PRINT "plastic case or a ceramic sandwich case. The style of case construe 
t i on af f {•cts" 

2530 PRINT "the heat transfer and must be known for the thermal model." 

2585 PRINT 

2590 PRINT "The methods of specifying the types of cases are listed below." 

2595 PRINT TAB < 15) ; "THESE APPLY TO THE DIPs ONLY" . 

2600 PRINT L I N C 2 ) 

2605 PRINT TABC10);" 1. ALL CERAMIC" 

2610 PRINT TABOO;-; "2. ALL PLASTIC" 

2615 PRINT TAB ' 10> ; "3. MAJORITY CERAMIC t USER SPECIFY WHICH APE PLASTIC)" 

2620 PRINT TAB(10>;"4. MAJORITY PLASTIC CUSEP SPECIFY WHICH APE CERAMIC'-" 

2625 PRINT TAB< 10) ; "5. USER SPECIFY CASE STYLE FOR EACH DIP COMPONENT" 

2630 Aria 1 = " 1 " 

2635 INPUT "ENTER YOUR CHOICE < 1 . 2 , 3 , 4 , 5 ) " , An ■ 1 
2640 IF AnsK /"BACKJJP" THEN 2670 

2645 PRINT L INC 20 SPA' 20 ), "UNABLE TO BACK UP WITHOUT LOSING THE BOARD " , L I H C 5 ) 
2650 BEEP 
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BOARDS 



DATA INPUT FOR THEFMELEX 



2655 DISP "Prsii CONT when ready to continue" 

2660 PAUSE 

2665 GOTO Case_typ 

2670 Ans=VAL ( Ans d > 

2675 IF (Ans<l> OF-: < Ans '5> THEN Case_type 
2680 OH Ana GOSUB C ere , P 1 as , Me e r e , Mp 1 as , Anv 

2635 Caae_cha: Sd= "COR PECT I QNS TO CASE STYLE FOP DIP PACPAGES" 

2690 GOSUB Page he ad 

2695 PRINT " On r he screen below is a line that represents, the case type For 
all regions." 

2700 PRINT "Those regions that contain a plastic cased component, are shown in in 
verse vidio" 

2705 PRINT " < 1 ) while all others are shown in normal video < 1 . " 

2710 PRINT 

2715 PRINT "To make corrections enter the region # that is incorrect , if no chang 
es no entry." , LI NCI } 

2720 FOR 1=1 TO Nreg 

2725 IF I MOD 10=1 THEN PRINT SPfl(lO); 

27 30 IF * Cases [ I , I ] = " 0" '• AND U; = 10> THEN PRINT I ,* 3 P ft I > ; 

2735 IF < Cased C I , I 3 = "0 ’ ■ AND < I 1 y' THEN PRINT SF ft( 1 > : I 1 SPA ■ 1 > ; 

2740 IF (Cases C I , I ] = " 1 " > AND <I>=10) THEN PRINT CHPd < 129); I ; C HR d ( 1 23 > ; SPA ( 1 ) 

2745 IF (Cased! I , I3*"l " ) hND <I<10> THEN PRINT C HRd ( 1 29 • ; " " ; I ; CHRd ^ 1 28 > ; SP A 

<i>; 

2750 IF I MOD 1O = 0 THEN PRINT L I N (1 ■ 

2755 NEXT I 

2760 PRINT 

2765 Ansd="NULL" 

2770 INPUT "ENTER REGION # TO CHANGE OP PRESS CONT WITH HO ENTRY IF OK", fins d 
2775 IF ( Ansd= " NULL" > OR ( Ans d= " BACKJJP " > THEN RETURN 
2730 Peg 1 =VAL < Ansd > 

2735 IF ( R e g 1 - 0 ; AND t Pe g 1 ( =Nr eg > THEN 2300 
2790 GOSUB Err in 
2795 GOTO Lass c h a 

230C Cased ! Peg 1 , Reg 1 2 = VALd ( I NT ( ( VAL < Cased ! Pegl , Reg 1 ] > + 1 > MOD 2>> 

2805 GOTO Case_cha ! END OF Case_cha 

2810 ' 

2815 Cere: 1 ALL CERAMIC CASE? 

2820 FOR 1=1 TO Nreg 
2825 Cased! 1 , 1 3 = " 0 " 

2830 NEXT I 

2835 RETURN !£ND OF CERE 

2840 ! 

2845 PI asl! ALL PLASTIC CASES 
2850 FOP 1=1 TO Nreg 
2855 Cased! I , I ]=" 1" 

2860 NEXT I 

2865 RETURN * END OF PI as 

2870 I 

2875 Mcere: Sd= " MA JOR I T Y CERAMIC ... SPECIFY WHICH APE PLASTIC" 

2830 GOSUB Page he ad 
2335 FOR 1=1 TO Nreg 
2890 Cased! I , n = "0" 

2895 NEXT I 
2900 Ansd= " NULL " 

2905 INPUT "ENTER THE NUMBER OF PLASTIC COMPONENTS " , Ans d 
2910 IF Ansd=“ BACKUP" THEN Case_typ 
2915 Ans = VAL (Ansd > 
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BOARDS 



DATA INPUT FOR THEPNELEX 



2920 IF (fldi tn AMD < Ans<Nreg+l > THEM 2935 

2925 GOSUB Err in 

2930 GOTO Me ere 

2935 IF Ans=Nreg THEM P 1 as 

2940 FOR 1=1 TO Ana 

2945 AnaS = 11 MULL" 

2950 INPUT “ENTER THE REGION NUMBER THAT CONTAINS A PLASTIC COMPONENT ", fins t 

2955 IF finsl*" BACK JJP " THEM Metre 

2960 Reel =VAL •' Ansi > 

2965 IF <Regl>0.> AND ■ Fegl ‘ Mreg+1 > THEN 2930 

2970 GOSUE Emin 

2975 GOTO 2950 

2930 CasefCPegl,Pegl3=" 1" 

2985 NEXT I 

2990 RETURN ! END OF Me ere 

2995 ! 

3000 Mpl as: S-f = " MAJOR I TV PLASTIC ... SPECIFY WHICH ARE CERAMIC" 

3005 GOSUB Pig* he ad 
3010 FOR 1=1 TO Nr eg 

3015 CaseSC I , I 3 = " 1" 

302O NEXT I 

3025 Ana I =“ NULL" 

3330 INPUT "ENTER THE NUMBER OF CERAMIC COMPONENTS ", fins $ 

3035 IF fins** "BACKUP" THEM Case ? yp 
3040 A n a = I N T < V A L 1 A n a > 

3045 IF < fins , O'* AND • Ana-' Mr e g* 1 ■ THEM 3069 

3050 GOSUE Err in 

3055 GOTO Mpl aa 

3O60 IF Ana=Mreg THEM Cere 

3065 FOR 1=1 TO Ana 

307O Ans**" NULL" 

3075 INPUT "ENTER THE REGION NUMBER THAT CONTAINS A CERAMIC COMPONENT ", Ana-F 

3030 IF flns •$ = " BACK_UP " THEM Up las 

3035 IF <Regl>0> AND -.RegK Mreg+ 1> THEN 3100 

3090 GOSUB Emin 

3095 GOTO 3075 

3100 Case* CRegl ,Regl 3="0" 

3105 NEXT I 

3110 RETURN ! END OF Mpl as 

3115 ! 

3120 Any: 1 EACH REGION MUST HAVE IT'S DIP SPECIFIED 
3125 FOR 1=1 TO Nreg 

3130 IF < I type UK1) OP - 1 t vc« C I > > 1 O > THEN 3160 

3135 DISP "IDENTIFY CASE TYPE FOP REGION I: VENTER 1 FOR PLASTIC AND 0 FO 

R OTHEP>"; 

3140 INPUT CaaeIC 1,13 

3145 IF CCaseSC 1 , 1 3 = " 1 " > OP < Case ft 1 , 1 3 = "0" ) THEN 3160 

3150 GOSUB Emin 

3155 GOTO 3135 

3160 - NEXT I 

3165 RETURN ! END OF Any 

3170 ! 

3175 Wha*._nou: ! THIS SECTION IS THE FINAL SECTION 
3180 S T = " U H A T NOW ?" 

3185 GOSUB Pagehead 

3190 PRINT " You have completed one cycle through BOARDS and ha'-'e the follow 
i ng opt i cna : " 

3195 PRINT L I N < 1 ) , T A E 1 1 0 > , " 1 . FEPFCPM THERMAL ANALYSIS ON "jNarae*!" USING THERM 
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< 12 > 



BOARDS . . DATA INPUT FOP: THERMELEK; 



320G PRINT LIN( 1 > , TAB< 10) , "2. INPUT ANOTHER CIRCUIT BOARD DESCRIPTION FROM KEYED 
ARD." 

3205 PRINT LIN<1),TAB<10>, "3. INPUT ANOTHER CIRCUIT BOARD DESCRIPTION FROM MASS 
STORAGE. “ 

3210 PRINT LIN(l), TAB<10>, "4. TERMINATE." 

3215 ! 

3220 Ans$="l" 

3225 INPUT "YOUR CHOICE • 1 , 2 , 3 , 4 > " , Ans f 
3230 IF Ans $ * " BACK_UR " THEN RETURN 
3235 Ans=ABS(INT VflL 1 finis')) ) 

3240 IF A ns \ 5 THEN 3255 
3245 GOSUB Errin 
3250 GOTO 3220 

3255 ON Ans GOTO Therm 1 , F edo , Pe do , Ends 

3260 ThermljDISP "WORKING LOADING THERML TO PERFORM ANALYSIS OF Maine* 

3265 LOAD "THERML " , 35 
3270 STOP 
3275 ! 

3280 Redo: Map=Ans 
3285 RETURN 
3280 ! 

3295 EnderPRINT PAGE , L I N< 20 '■ , SPA v 1 0 > , " NORMAL TERMINATION " 

3300 DISP "WORKING LOADING STANDARD KEY DEFF I M I T I QMS " 

3305 GCLEhR 

3310 LOAD KEY "STDKEY" 

3315 END ! END OF What_now 

3320 ! 

3325 P age he ad : 1 THIS SECTION PLACES DESIRED HEADING ON A BLANK CRT 
3330 PRINTER IS 16 

3335 PRINT PAGE , TAB ■, 34-LEN ( S S '* .2 *** " ; CHR* < 1 32 > ; S* ! CHP* ' 1 28 > ; " + *■*". LIN<2) 

3340 RETURN ! END OF PAGEHEAD 

3345 ! 

3350 Errin: i THIS SECTION ALERTS THE USER TO AN ATTEMPT TO INPUT BAD DATA 
3355 BEEP 

3360 DISP “INPUT OUT OF RANGE TRY AGAIN" 

3365 WAIT 15O0 
3370 BEEP 
3375 WAIT 1000 
3380 BEEP 
3385 RETURN 

3330 ! END OF Errin 

3335 ! 

3400 Leads e: GOSUB Error 

3405 Leads: ! THIS SECTION ALLOWS THE INPUT OF THE NUMBER OF LEADS BETWEEN REGION 
S 

3410 ON ERROR GOTO Le adse 

34 15 $•! = " ELECTRICAL LEADS 0° OTHER CONDUCTION Pfl’HS" 

3420 GOSUB Page he ad 

3425 PPINT " All circuit boards have thermally conductive materials applied t 
o their " 

3430 PRINT "surface either in the form of electrical leads dr as thermal conduct 
ion paths." 

3435 PRINT "Heat will be transferred through these materials between the regions 

3440 PRINT "if there are conduction paths (rails) present, they will have a much 
1 arger " 

3445 PRINT "effect than the electrical leads and the electrical leads will be ne 
g 1 e c t e d " 
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BOARDS 



DATA INPUT FOP THERMELE7 



3450 PRINT "in the thermal model." 

3455 AnsX="N" 

3460 INPUT "ARE THERE METAL CONDUCTION PATHS PRESENT ( N or Y)",fln-I 

3465 IF A n s X = " B A C 1 UP" THEN RETURN 

3470 IF UFCX(AnsX[l", 1 3 > = " N " THEN Lead*_size 

3475 IF UPCXkAnsXC 1 , 1 ] > ="Y" THEN 3490 

3430 COSUB Errin 

3485 GOTO Leads 

3490 Sx = "CONBUCT I ON F AILS" 

3495 A 1 X = " CPU ' s " 

3500 GOSUB Pagehead 

3505 PRINT " Thermal conduction rails are much wider arid thicker than the ele 
c t r i c a 1 " 

3510 PRINT "leads plated onto the circuit board. The thickness is generally uni 
Form , " 

3515 PRINT “however; the rails are typically not of uniform width. It is theref 
ore" 

3520 PRINT "necessary to define a concept of Conduction Path Units (CPU's). The 
rail wi dth" 

3525 PRINT "will be modeled in terms of an i nt e qe r number of CPU's. You will 
be asked to" 

3530 PRINT "input the thickness (mm) of a CPU and the width will automatically d 
e f au 1 1 to" 

3535 PRINT ".1 mm such that a rail of width 1.7 mm can be modeled as 17 CPU's." 
3540 PRINT 

3545 PRINT " For you r convenience the thermal conductivities (watts* M-C) of t 
hree alloys" 

3550 PRINT "commonly used as material for conduction rai 1 s are given below:", LIN 

<1 > 

3555 PRINT " A 1 5052= 133.2 U/M-C AL 6101= 215.7 U M-C Cu 113 = 333 

.8 U/M-C" 

3560 PRINT 
3565 0!de=0 
3570 Ans X = " 1 " 

3575 INPUT "ENTER THE THICKNESS CF THE THERMAL CONDUCTION PATHS or CPU's (mm)", A 
nsX 

3530 IF Ans X=" BACK UP" THEN Leads 
3535 T! =ABS(VAL< AnsX O 

3590 PRINT "THICKNESS OF CPU's = ";T!;"mm" 

3595 AnsX=" . 1 " 

3680 ! INPUT "ENTER THE WIDTH OF THE CPU's 'mm)", AnsX 
3605 IF AnsX = " BACK JJF‘ " THEN 3570 
3610 Ml =AES(VftLtflnsX) ) 

3615 A 1 = - HI *• T 1 ! #»*####*# NOTE THAT THIS AREA IS NEGATIVE FOP CPU s 

3620 PRINT "WIDTH OF CPU's = ";WlJ"»m AREA FOR HEAT TRANSFER FER CPU =" 

; - A 1 ; " mm'" > 2 " 

3625 AnsX=" 133.2" 

3630 INPUT "ENTER THE THERMAL CONDUCTIVITY OF THE CONDUCTION PATHS ( U* M-C ) " , Ans X 
3635 IF AnsX= "BACK UP" THEN 3570 
3640 K1 =ABS( VAL( AnsX > ) 

3645 PRINT " THERMAL CONDUCTIVITY = " ; K 1 ; "W at t s/fl-De a C" 

3650 WAIT 15C0 

3655 GOTO Leads enter ! END OF CPU 

3660 ! 

3665 Leads s i ze : 3X= "NUMBER OF LEADS" 

3670 GOSUB Page he ad 
3675 A 1 X = " LEADS " 

3680 PRINT " This section of the program allows the entry of the size and con 
due t i v i t y of " 
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BOARD 1 : 



DATA INPUT FOP THFPMELEY 



3635 PRINT "leads <traces> on the surface of r.he circuit board. These act as bo 
t h" 

3690 PRINT "electrical and t hernial conductors between the regions. When enter in 
g the width" 

3695 PRINT "figure an average width for all the leads." 

3700 PRINT LINC1>," The thermal conductivity will default to pure copper <334 
Wat t s-'N-De g C>" 

3705 PRINT "with no entry." 

3710 01 de=0 

3715 flns-f = " . 0 33" 

3720 INPUT "ENTEP THE THICKNESS OF THE LEADS •. TRACES? ON THE CIRCUIT BOARD (mm) 
" , AnsS 

3725 IF AnsI="BACKJJP" THEN Leads 
3730 T 1 =ABS VALCflns X ■ > 

3735 PRINT L I M< 1 > . "LEAD THICKNESS =" ; T1 ; "mm" 

3740 A n s 4 = " 1 " 

3745 INPUT "ENTER THE AVERAGE WIDTH OF THE LEADS < TRACES > ON THE CIRCUIT BOAP 
D < mm ) " , Ans I 

3750 IF Ansi=" BRCKJJP" THEN 3715 
3755 W 1 = A B $ ( V A L H n s S > ) 

3760 A 1 = W 1 * T 1 

3765 PRINT " AVG WIDTH OF LEADS = " ;Wl;"mm AREA FOR HEAT TRANSFER PER LEAD 

= " ; A 1 ; " m m A 2 " 

3770 Ans f = " 384 " 

3775 INPUT "ENTER THE CONDUCTIVITY OF THE TRACES ''DEFAULT COPPER = 384U at t. s-' M-d 
egK ) " , Ans* 

3730 K 1 =ABS ( VAL < Ansf > ) 

3735 PRINT "CONDUCTIVITY OF THE LEADS = " ; K 1 J " W at t I ' M-C " 

3790 WAIT 1500 

3795 Le ad s_en ter; 3* ■" NUMBERS OF ":iAlf&" BETWEEN REGIONS" 

3SO0 GOSUB Pagehead 

3305 PRINT " This section of the program allows the entry of the number of th 
e " J A 1 S 

3310 PRINT "that cross each of the internal region boundaries. If these are on 
both sides" 

3815 PRINT "of the board add both numbers together. The flashing cross will mow 
e to the" 

3320 PRINT " appropr i at e location on the screen but the graphics picture will not 
remain, " 

3825 PRINT "If a hard copy of the graphics is needed for a guide, recall that k3 
will" 



3830 PRINT "provide one. The picture will return to the screen while the number 
of " ; A 1 f 

3835 PRINT "is labeled and the flashing cursor will move to the ne t location to 
be ent e-red. " 

3840 PRINT L I N i' 1 > , " DO NOT USE kO < BRCKJJP > WHILE THE NUMBERS OF " ; M 1 X ; " APE BE 
ING ENTERED!" 

3345 PRINT L I N ( 1 > , " There will be correction opportunities later." 

3850 DISP "PRESS CONT WHEN READY TO START WITH REGION #1" 

3855 PAUSE 

3360 GOTO Leads_in 

3365 Le ads_l abe 1 : ! HERE THE NUMBERS OF LEADS OP CPU 'S IS WRITTEN ON THE SCREEN 

3870 CSIZE 2.5 

3875 GRAPHICS 

3880 LABEL USING "K";N1 

3885 PEN -l 

3390 LDIR -PI 

3395 LABEL USING "K" ;N1 
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BOARDS .. DATA IHPUT FOR T HERMEL E. : 



PI THEM WAIT 500 



3900 PEN 1 
3905 LDIP 0 
3910 IF 01 de- 
3915 RETURN 

3920 Leads oldlCSIZE 2.5 
3925 S •>*“ NUMBER'S OF LEAD'S" 

3930 A1*="LEAPS" 

3935 IF AKO THEN Ai*="CFU s" 

3940 GOSUB Pagehead 
3945 PRINT “ This section of 
ng " 

3950 PRINT "boundaries of regions b 
that " 

3955 PRINT " represent the informations stored 
to make" 

3960 PRINT 
equi red. " 

3965 PRINT 

3970 ! PRINT "It' you do not des’re 
enter any " 

3975 ! PRINT "number before you press CONT 
3930 DISP "PRESS CONT WHEN READY TO VIEW DATA " 

3935 Leads_i n : FOP 1=1 TO Nreg 

3990 J=1 ! HERE IS BOTTOM OF THE REGION 

3995 GRAPHICS 

4000 IF I + Nxr< =Nr eg THEN 4025 

40O5 NI<J,I)«H1*0 ' HERE IS BOTTOM OF THE BOARD WHERE 

4010 MOVE '<• I >, Y< I 4 5 ♦Hr 

4015 GO'iUB Le ad s label 
4020 GOTO 4165 



the progr am displays the 



placing numbers on i 



that information and re-store 



y i e w t he 



4 025 IF < I type d >< >10 > OR U t ype < I +Nxr ) , 1 0 ; THEN 4 100 
4030 MOVE xa ) + . 32-Lr, Y< I >-.45-Hr 
4035 IF 01 de< /PI THEN 4050 
4040 N 1 = N 1 C J , I ) 

4045 GOTO 403O 

4050 POINTER X< I '• + . 32-Lr, Y< I >-. 45-Hr, 2 
4055 WAIT 2000 

4060 DISP "ENTER THE NUMBER OF " ; A 1 * ; " BETWEEN REGIONS 
; " AND "*,CHR$< 132); I+N*r : CHR* C 123) J 
4065 INPUT HI 
4070 N1*ABS<N1> 

4 075 N1 < J*2, I+N*r)*Nl < J , I )=Nl 
4080 GOSUB Le ads_l abe 1 

4035 MOVE X < I + Nxr > + . 32 ■’■Lr , Y< I +M/r ) + . 42*Hr 

4090 GOSUB Leads label 

4095 GOTO 4165 

4100 MOVE X<I),Y<I)-.45-Hr 

4105 IF OldeOPI THEN 4120 

4110 N 1 = N 1 < J , I > 

4115 GOTO 4150 

4120 POINTER X<I>,Y<I)-.45*Hr,2 
4125 WAIT 20O0 

4130 DISP "ENTER THE NUMBER OF " ; A 1 S ; " BETWEEN REGIONS 
AND " J CHR* < 132> ; I + Nxr ; CHP* •! 123/ ; 

4135 INPUT HI 
4140 H 1 = A £ 3 < H 1 > 

4145 Hl( J, I ;=N1 < J + 2, I +Nxr>=Nl 
4150 GOSUB Le ads_l abe 1 



< 16 ) 



numbers of leads cross 
he sides of the region 
You will be all owe 
in on mass storage if 

rri’ng number of leads, 



CONNECTORS WILL GO 



" ; CHP*< 132; ; I : CHR*< 123 



" ; c HRs a 32 > ; i : chr* < 12 s 
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BOARDS .. DATA INPUT FOP THERMELEX 



4155 MOVE X(I+hKr),Y<I+Nxr)+.42tHr 
4160 COSUB Li ads_l abe 1 

4165 J=2 ! HERE IS FOR RIGHT SIDE OF REGION 

4170 IF I MOD NxrOO THEN 4185 

4 175 HH2,I)»N1»U ! HERE IS THE RIGHT SIDE OF THE BOARD 

4180 MOVE I> + .42«rLr, Y-. I> 

4135 GOSUB Leads laoel 
4130 GOTO 4340 

4185 IF <Ityp*U><>9> OR < I type (1 + 1)09) THEN 4270 
4200 MOVE X U > + . 42*Lr , Y U > - . 32-*Hr 
4205 IF OldeOPI THEN 4220 
4210 N 1 =N 1 < J , I > 

4215 GOTO 4250 

4220 POINTER XU > + . 42*Lr , Y < I) - . 32*Hr , 2 
4225 WAIT 2000 

4230 DISP "ENTER THE NUMBER OF ' ; fl 1 f ; " BETWEEN REGIONS " ; CHPf < 1 32 > ; I I CHRl 1 1 23 > 
and " ; CHRs 1 1 32 > : i -*■ i : chrs < 1 28 >: 

4235 INPUT N 1 
4240 HI = H B S < N 1 ) 

4245 HI < J , I ) =N 1 < J -*• 2 , I + 1 > =N 1 

4250 GOSUB Leads label 

4255 MOVE X U ♦ 1 )' - . 42*Lr , Y < I ♦ 1 > - . 32 »>Hr 

4260 GOSUB Lsads_lab*l 

4265 GOTO 4335 

4270 MOVE X<I) + .42*Lr, YU) 

4275 IF 01 de< -PI THEN 4230 
4280 N 1 = H 1 ( J , I ) 

4235 GOTO 4315 

4230 POINTER X U ) * . 4 2<*Lr , YU ; , 2 
4235 WAIT 2000 

4300 DISP "ENTER THE NUMBER OF ' ; A 1 T ; " BETWEEN REGIONS " ; C HR t < 1 32 ) ; I J CHRT <1 2S .■ 
; " and “ ;CHPS' 132 >; iu; C hrs U2S u 
4305 INPUT HI 
4310 Nl=AESOn) 

4315 GOSUB L> ads_l abe 1 
4320 HI < J + 2, IM > = N1 ■: J, I 1 = N 1 
4325 MOVE X (I * 1 >- . 42*Lr , Y U > 

4330 COSUB Le ads_l abe 1 

4335 J=3 ! HERE IS FOR TOPS OF EACH REGION 

4340 IF I >Nxr THEN 4360 

4345 N1 <J, I >=N1 =0 

4350 MOVE XU),Y<I; + .42*Hr 

4355 GOSUB Leads_lao.il 

4360 J =4 ! HERE IS FOR LEFT SIDE OF REGION OR BOARD 

4365 IF <I-1> MOD N.<r=0 THEN 4375 

4370 GOTO 4330 

4375 N 1 < J , I ) = N 1 = 0 

4330 MOVE X<I>-.42*Lr,Y(I) 

4335 GOSUB Leads_label 
4330 NEXT I 
4335 GRAPHICS 
4400 MOVE 160,130 
4405 LOPG 6 
4410 CSI2E 3 

4415 LABEL USING "K"; " # OF "fcAl-f 
4420 DRAW X(Hxr>, Y< N.' r > +Hr-'2 
4425 POINTER 0,0,0 
4430 WAIT 5000 
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( 17 ) 



BOARDS 



DATA INPUT FOR THEPMELEX 



4405 EXIT GRAPHICS 
4440 RETURN 

4445 Leads_chae: GOSUB Error 

4450 Leads”cha: S*="CORPECTIONS TO NUMBER OF "iflU 

4455 ON ERROR GOTO Leads_chae 

4460 Check =0 

4465 G03U3 Pagehead 

44?0 PRINT " V'iju may now make correct iorii to the numbers of " ; A 1 f ; " crossing 
the boundar i es. “ 

4475 PRINT "Recall that k3 will provide a hard copy of the graphics if needed 
while k 2 " 

4430 PRINT "will return the graphics picture to the screen for visual checking. 
4485 PRINT 

4 4 SO PRINT "In responce to the prompts below, INPUT Reg #, Reg 4, CORRECT # OF 
■;fil*,LIN<l; 

4495 PRINT "IF THERE ARE NO CHANGES PRESS COMT KEY WITH NO INPUT." 

4500 LORG 5 

4505 IF Check >0 THEN EEEP 

4510 IF Check. 0 THEN PRINT LIIK3) , " CHECK OVER THE FIGURE , YOU HAVE MADE AN IL 
LOGICAL CHOICE OF REGIONS." 

4515 Check=Check+t 
4520 Peg 1 = 0 

4525 DISP "REG 4 , PEG 4 , CORRECT 4 OF ";flli; 

4530 INPUT Re g 1 , F'eg2 , He wnT 
4535 IF Reg 1 =0 THEN RETURN 

4540 IF ( Peg 1 >Nr eg > OR t F'e g2 ' Nr e a ■ THEN 4465 
4545 IF R e g 2 R e g 1 THEN 4565 
4550 Dum=Peg2 
4555 Reg2=Reg 1 
4 560 R e g 1 = D u r» 

4565 J=0 

4570 IF < R e g 1 ♦ 1 = R e g 1 .> AND (Regl MOD Nxr ' >0 ■ THEN J = 2 
4575 IF Regl+Nxr=Reg2 THEN .J=l" 

45S0 IF J<>0 THEN 4655 

4535 EXIT GRAPHICS 

4530 SS= " ERROR IN CORRECTIONS" 

4595 GOSUB Pagehead 

4600 PRINT LIN(5>, SPAC5 ' , " THOSE TWO REGIONS DO NOT CONNECT TRY AGAIN" 

4605 BEEP 
4610 WAIT 2500 
4615 GOTO Leads_cha 

4620 ! HERE IS THE CORRECTION SCHEME 

4625 Leads_er ase : PEN “1 

4630 GOSUB Le ads_l abe 1 

4635 N 1 =Newn 1 

4640 REN 1 

4645 GOSUB Leads_label 
4650 RETURN 

4655 IF J=2 THEN 4745 

4650 IF < Itype <Pegl >< >10 OR ■. I type (Rsg2> . > 1 0 > THEN 4705 
4665 MOVE X < P eg 1 > + . 32*Lr , Y< Reg 1 > - . 4S*Hr 
4670 N1=N1 < J, Regl > 

4675 GOSUB Le ads_er as e 

4630 MOVE X<Peg2.' + .32*Lr , YtReg2> + .42*Hr 
4635 N 1 = N 1 C J + 2 , R e g 2 ' 

4690 GOSUB Le ads_er as e 

4695 N1 <1 , Regl >=N1 < 1 , Peg2>=Newnl 
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< 13 > 



BOARDS 



path ihput for t hef'mele; i 



470O 


GOTO Leads change 




4705 


MOVE XCReql ^ , Y( F'eq 1 . 45*Hr 




4710 


N1=N1 <J,Regl) 




4715 


COSUB Leads erase 




4720 


MOVE X < Reg27, Yv Reg2 ) + . 42*Hr 




4725 


N 1 =N 1 < J*2, Peg2> 




4730 


GOSUB Leads erase 




4735 


N 1 < J,Regl .'=Nl <J*2, Reg2>»Newnl 




4740 


GOTO Leads cha 




4745 


! HERE WE APE TO CORRECT THE RIGHT SIDE 


OF REGION: 


4750 


IF < 1 1 y pe < Re g 1 X >9 '• OR < 1 1 ype v Peg2 ) < >9 ) 


THEN 4795 


4755 


MOVE XcPegl - ' + .42-»Lr,Y(Redl)-.32*Hr 




4760 


N 1 = N 1 < J , R e g 1 ) 




4765 


GOSUB Leads_erase 




4770 


MOVE X ( Re g2~- . 42-Lr, Y<.Peg2>-. 32*Hr 




4775 


Ml =N1 < J + 2, Reg2 ) 




4730 


GOSUB Leads erase 




4785 


N 1 < ■ J , P e g 1 > = N 1 < J + 2 , Peg2 ) =Ne wn 1 




4790 


GOTO Leads change 




4795 


MOVE X\Pegl ;+. 42*Lr , YtPegl ) 




4300 


N 1 = N 1 <J, Regl > 




4805 


GOSUB Leads erase 




4810 


MOVE X * Peg2 ) - . 42*Lr , Y<Pe g2) 




4815 


N1=N1( Jt-2, Reg2 




4320 


GOSUB L e ads erase 




4825 


HI < J, Regl )=N1 f J+2 , Reg2 **Newnl 




4330 


GOTO Leads cha 





4335 Temp i n : ! THIS SECTION INPUTS THE KNOWN TEMPERATURES OP POWERS OF ELEMENTS 

4840 EMIT GRAPHICS 

4345 $!=" TEMPERATURES OR POWER LEVELS OF COMPONENTS" 

4350 GOSUB Pagehead 

4355 PRINT " The thermal model used by this program assumes each component to 
be a heat" 

4330 PRINT "source for which the user specifies either the maximum junction temp 
erature or" 



4365 PRINT 
junct i on" 
4870 PRINT 
c al cul ated 
4375 PRINT 
c t i on" 

4380 PRINT 
4835 PRINT 
4890 PRINT 
del and " 
4895 PRINT 
ever , " 

4900 PRINT 
he same as 



the rate of heat generation within that component. When the maximum 
temperature is specified, the maximum stead 1 ,' state power levels are 
When the rate of power di si pat ion is specified, the steady state jun 
temperature- are calculated." 

The program uses the component, surface temperature in the thermal mo 

therefore requires a case to junction thermal resistance C R j c > ; how 

if zero is specified then the surface temperature is assumed to be t 



4905 PRINT "the junction temperature. When components span two regions, enter h 
a 1 f the " 

4910 PRINT "component power for each region." 

4915 PRINT LIN<1)," All entries must be in Deo C or Watts and deg C-'Watt s. " 
4920 IF 0 1 d e O P I THEN 4945' 

4925 DISP "PRESS CONT WHEN READY TO VIEW THE DATA FROM " ; CHPf < 1 32 ' ; Name S J CHP S < 
123);" FOP POSSIBLE CHANGES" 

4930 PAUSE 

4935 GOTO Ternp_cha 
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BOARDS 



DATA INPUT F 0 R THERMELEX 



4940 ! 

4945 flnsl 3 "l" 

4950 MAT P o w = Z E R 
4955 MAT Tj=ZER 
4960 MAT R j_c = Z£R 

49S5 INPUT "DO VOU DESIRE TO SPECIFY POWER LEVELS (1) OR TEMPS (2) ^",Ansf 
4970 IF AnsS="BACh UP" THEN RETURN 
4975 Ans=VAL<AnfS>~ 

49S0 IF < fln»= 1 ) OP tAns=2> THEN 4995 

49S5 GOSUB Err in 

4990 GOTO Tempi n 

4995 IF Ans = 1 THEN Pcw__in 

50O0 T emp_i n : T em_f o 1 =0 1 HERE INPUT TEMPEFATUFES WILL SOLVE FOR POWERS LATER 

5005 sx= "I nput of junction temperatures " 

5010 GOSUB P age he ad 

5015 PRINT " You ire now entering junction * e mper atures- ',deg C) and junction 
to ciie" 

5020 PRINT "thermal resistance (deg C/Wat t >. Do no* ufe kO (Back-up) option w 
hi 1 e entering" 

5025 PRINT "the data." 

5030 PRINT LINO) 

5035 PRINT " I Tj<I) P j c < I > 

5040 FOP 1=1 TO Nr eg 

5045 IF I type < I '*0 THEN 5060 

5050 DISP "FOP ELEMENT IN F EG I ON " ; I ; " ENTER T junc , F'j_c " 5 

5055 INPUT T j < I ■ , R .j ■ I ) 

5060 PRINT SF'ftv 1', i7TjCI),Pj c(I) 

5065 Tj(I/=Tj(I)+273 

5070 NEXT I 

5075 GOTO Tempj:hi ! END OF TEMP_I N 

50S0 Pow_i n : T em_ = o 1 = 1 ! HERE INPUT POWER LEVELS WILL SOLVE FOR TEMPERATURES 

5035 S f = " INPUT OF COMPONENT POWER LEVELS" 

5090 GOSUB Page he ad 

5095 PRINT " You are now entering component power d i if- i p at i on (Wat*) and jun 
c t i on to" 

5100 PRINT "cafe thermal resistances (deg C'Watt). Do not use k9 <Back_up whi 
le entering" 

5105 PRINT "the data." 

5110 PRINT LINO) 

5115 PRINT " I Pow( I ) Rj c < I ) " 

5120 FOR 1=1 TO Hr eg 

5125 IF 1 1 ype < I > = 0 THEN 5140 

5130 DISP "FOP ELEMENT IN REGION # " ,* I ; " ENTER Pow< I ) , P j_c < I > " ; 

5135 INPUT P o w ( I ) , R j c(I) 

5140 PRINT SP A < 1 > , I , Row < I ) , R j c(I) 

5145 T j ( I ) =273 

5150 NEXT I 

5155 GOTO Temp c ha ! END OF POW_lN 

5160 ! 

5165 Ternp ch-ae! GOSUB Error 

5170 Temp~cha: ! THIS SECTION ALLOWS CHANGES TO THE TEMPS OR POWER LEVELS 
5175 SS="DATA FOR "tPietlbl* 

5180 GOSUB Fagehead 

5185 PRINT " The data lifted below are the c u r r e n t. values for the variables 
spec i f i ed. " 

5190 PRINT " REGION 4 Tjunc (DeqC) Power (Watt 3 ) Re as e- i ' U -'C ) " 

5195 FOP 1=1 TO Nr eg 

5200 PRINT TAB ( 2 ) , I ; TAB', 17 ■ , T j ■ I ) -273 5 TAB ( 33 ) , Row’ I J TAB < 43 > , P j_c ' I > 
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BOARD' 



. DATA INPUT FOR THEPMELEX 



5205 NEXT I 

5210 PRINT L I N < 2 ) 

5215 IF Sk i p = P I THEN GOTO 5360 

5220 PRINT "Us? the DISPLAY up-arrow OR down- arrow to move the data list as requ 
i red" 

5225 Check =0 
5230 flna**"0" 

5235 INPUT "ANY CHANGES “> INPUT PEG # IF YES OR PRESS CONT FOR NO CHANGES ", Ans 

t 

5240 IF Ans T = " BACK UP" THEN Temp in 
5245 R e g 1 = I N T < V AL < fin s 1 ^ ) 

5250 IF (Regl >-l ) AND (Reo KNrsg* 1 > THEN 5265 
5255 GOSUB Err in 
5260 GOTO 5230 

5265 IF < Peg 1 =y > AND (Ch«ck*0) THEN Hard 
5270 IF <■ P e g 1 = 0 > AND <Check<,0> THEN Temp_cha 
5275 DISP "TO CHANGE VALUES IN REGION #“;P«qi; 

52S0 Check = 1 

5235 IF T em_so 1 = 1 THEN DISP "ENTER Power , P.i c"J 

5290 IF Tem_so1=0 THEN DISP "ENTER Tjunc , Pj~c H i 

5295 IF T em_s o 1 = 1 THEN INPUT Pow < Pe g 1 ’ , Rj c < Peg 1 1 
530O IF Tem_sol=0 THEN INPUT Tj > Pe g I ■ , R j_c • Pe g I > 

5305 IF Tem_*ol=0 THEN T j • Pe g 1 '• =T j Pe g 1 > +273 

5310 PRINT LINC 1 ; , TfiB< 2 > , Fegl : TAB < 1 7 ) j T j ( Re q 1 > -273 : T AE 1 33 :■ , Row < Peg l.'- ; TAB 1 ' 43 > , P j 

_c<Regl);" *•-* CHANGE 

5315 GOTO 5230 ! END OF TEMP_CHA 

5320 Hard: ! THIS SECTION PRINTS A COPY OF THE INPUT DATA IF REQUESTED 
5325 A n s 1 = " N " 

5330 INPUT "DO YOU DESIRE A PRINTED COPY OF THE DATA ABOVE < N or Y>", final 

5335 IF Ansi* "BACK UP" THEN Temp_cha 

5340 IF fins 1= " N " THEN RETURN 

5345 PRINTER IS 0 

5350 Sk i p=P I 

5355 GOTO 5190 

5360 PRINT TAP < 39-LEN • "THE ABOVE DATA IS FOP "3.P i c t. 1 b 1 1 > • 2 ) , " THE ABOVE DATA IS F 
OP ";Ptctlbll,LIM<l> 

5365 PRINT "BOARD LENGTH (defined along air f 1 owi*“ : Bdl ; "mrn HEIGHT ="; 

Bdh ; "mm" ,LIN<1> 

5370 PRINT "BOARD THICKNESS*" ; Th i ck_b; "mm CONDUCTIVITY =" 

;Kb; "Watts/M-K" , LIN<2) 

5375 PRINT "THE BOARD MODEL ASSUMES “ ; A 1 X J " AS CONDUCTION P-THS WITH AN AREA OF 
";ABS<A1 >; " mm A 2 " , L I N I > 

53S0 PRINT "THERMAL CONDUCTIVITY OF THE " ; A 1 1 ; " *“;Kl;" Watt.a^M-C" 

5335 PRINT " 

" , LIN<2) 

5390 PRINTER IS 16 

5395 RETURN . ! END OF Hard 

5400 ! 

5405 Stowe! GOSUB Error 

5410 Stow: ! THIS SECTION PLACES THE DESCRIPTION OF THE CIRCUIT BOARD ON TAPE 

5415 ON ERROR GOTO Stowe 

5420 31= " RECORD BOARD DESCRIPTION ON MASS STORAGE" 

5425 GOSUB Pagehead 

5430 PRINT " You may record all the data concerning the circuit board on any 
aval 1 ab 1 e " 

5435 PRINT "mass storage device. This allows any user t * retrieve * he descripti 
on at some" 

5440 PRINT "later time without » he need to input all the details. This option o 
ccurs both" 
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['AT A INPUT FOR 1 HEP MEL EX 



3445 PRINT "before an d after the thermal analysis. Enter desired data file na 
me below," 

5430 PRINT "be sure to include the mass storage unit, specifier if not the defaul 
t . " 

5455 PRINT "For example : T 1 4 , : F8 , Y12 etc." 

5460 Ansf="Y" 

5455 INPUT "HO YOU DESIRE TO RECORD THE DESCRIPTION DATA <Y or N)‘‘" f Ans* 

5470 IF Ane* = " BACK JJP" THEN RETURN 
5430 IF lJF'C*<Ans*Cl , 1 J; = "N" THEN RETURN 
5435 IF 0 1 de< >F I THEN 5505 

5430 PRINT L I N ( 5 ) , "THE PRESENT DESCRIPTIVE TITLE FOR THAT BOARD IS " J C H R * ■; 1 3 2 ; P 
i ct 1 bl * ; CHR* k! 123) 

5435 EDIT "CHANGE THE TITLE OP PRESS CClNT WITH NO ENTRY FOR NO CHANGE. ", Pi C ♦ 1 b 

1 s 

5500 IF Pi ct 1 bl *=" BACKJJP" THEN Stow 
5505 Ansi = tlamet 

5510 EDIT "UNDER WHAT NAME DO YOU DESIRE TO STORE THE DATA (change below';-'" , Name 

* 

5515 IF Nar..e*=" BACK JJP" THEN 5430 
5520 IF Namef =AnsX THEN 5560 
5525 Ans f = " Y " 

5530 DISP "DOES A DATA FILE WITH AT LEAST " J 4O*Nreg+60O; " BYTES EXIST UNDER THAT 
NAME ( Y or N>?" ; 

5535 INPUT Ans * 

5540 IF Ans * = " BACK JJP " THEN 3305 
5545 IF UPC* (Ans* C 1 . 1 ] ) = "Y" THEN 5360 

5550 DISP "WORKING CREATING DATA FILE FOR "JName*;" THAT I S " ; 4 0 * H re q + 6 O O ; "BYTE 
S IN SIZE" 

5555 CREATE Name* , 1 , 40-*Nreg + 6O0 
5560 ASSIGN #5 TO Name* 

5565 PRINT LINO) , "WORKING WRITING BOARD DESCRIPTION OF " ; P i c t 1 b 1 * ; " ON MASS S 
TOPAGE" 

5570 DISP 
5575 READ *3, 1 

55S0 PRINT #5; Name* , Pi c 1 1 b 1 * , C ase* , Bd 1 , Bdh, Sc a, H\r , Hyr , Mreg, Thi ck_b, Kb, A1 , K1 , Ten 
_so 1 

5535 FOR 1=1 TO Nreg 

5590 PRINT #5; Itype-'I;.X< I > , Y< I > , T j < I > , Povm I),Rj_c(I) 

5595 FOR J=1 TO 4 
5600 PRINT #5 I N1 ( J, I ) 

5605 NEXT J 
5610 NEXT I 
5615 PRINT #5 ; END 
5620 ASSIGN #5 TO * 

5625 PRINTER IS 0 

5630 PRINT L I N ( 2 ) , " THE CIRCUIT LISTED BELOW IS STORED UNDER THE FILE NAME " ; C H F: * 
< 132) ; Name* ,* CHP* ( 123) 

5635 PRINT L I N ■' 2 ). SPA • 25 P i c 1 1 b 1 *, L UK 2 > 

5640 PRINT SPA (25 ) , "SAVE FOR YOUR RECORDS" 

5645 PRINT " 

, L I N <. 2 > 

5650 PRINTER IS 16 
5655 DISP 
5660 RETURN 
5665 Error! I 
5670 EEEP 
5675 WAIT 300 

5630 IF E R R H = 5 6 THEN Err_name 
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BOARDS 



DATA INPUT FOR ThERMELEX 



5635 PRINT LIN(20),$PA<10>, "ERROR NUMBER " ; EFRN ; " HAS GCUPFED IN L I HE" ; ERRL ; “ . PR 

ESS CONT WHEN READY - 

5690 DISP 

5695 BEEP 

5700 PAUSE 

5705 RETURN 

5710 ! 

5715 Err_nifrte:i THIS SECTION FOR IMPROPER FILE NAME 
5720 PRINTER IS 1 6 
5725 PRINT PAGE 

5730 M s u s I = " DEFAULT HA 3 S j T 0 R A G E " 

5735 FOR 1=2 TO LEhKN anie X :■ 

5740 IF Name* C I , I ] = " : " THEN 5755 

5745 NEXT I 

5750 GOTO 5770 

5755 Msu = X = Nar„e X [ 13 

5760 CAT Ms us X 

5765 GOTO 5775 

5770 CAT 

5775 PRINT L I N < 2 ' , CHPf < 1 32 > : N am* X C 1 , I - 1 3 : CHRX Cl 23 > ; " is NOT on "jMsusx;" with th 
at spelling. " 

5730 PRINT L I N •' 1 > , " CHECK OVER THE DIRECTORY ABOVE FOR CORRECT NAME 

5735 DISP “PRESS CONT WHEN READY" 

5730 PAUSE 
5795 RETURN 
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10 

20 



THERML THERMAL ANALYSIS PROGRAM FOR 



40 
50 
60 
70 
SO 
90 
100 
110 
120 
130 
140 
150 
160 
170 
ISO 
190 
200 
210 
220 
230 
240 
250 
260 
270 
0) , 
230 



* thermelex 

A SYSTEM OF PROGRAMS FOR THE HP 9345 
TO PERFORM THERMAL ANALYSIS OF ELECTRONIC CIRCUIT BOARDS 

*<*»^t«** + t*'*<>*t4*t** + -»*****'t***t**>**»**r*t + *'‘i***<**** + '»t**t 
* 

* 

* PREPARED AT THE NAVAL POSTGRADUATE SCHOOL MONTEREY, CA 

* R. A. FOLTZ LCDR USN 
HI I I THERML 



JUNE 1930 



OPTION FASE 1 
PRINTER IS 16 

DIM SIC 70] , T --pe I ■ O: 14 - , Hsus I [22] , Y1 bl IC253 

COM Map. Name! .Pic* 1 b 1 I C 5 O ] , C we ft 503 . Bd 1 , Bdh, Sc a, Thi ck_b , Tem_sol , 
COM SHORT Xv50 > . Y ■ 50) , T j <50 > , Row 50> , Pj_c v50> 

COM INTEGER 1 t vpe <- 50) , M 1 1 4,50), N\r , Nyr , Nre g 
INTEGER N p i n < 1 5 ) 

SHORT Ae < 5 0 > , ?e_c end '.50) , Re_convt 50 > , Rt or. e_b < 50. * , Pt op_r < 50) , P 1 < 4 , 
Pt>_conv<50’* 

SHORT Wide set 1 4 . , Le ns e t cl 4 ,< 



290 SHORT Le< 50 > , Me 50 ■ . He <. 50 ' , Te < 50) , T ai r ( 50 > , A ■' 50 , 50 ^ , B < 50 ) , Tb • 50) 

300 M ap = I NT ■' M ap ) 

310 IF < M ap = 0 ) OP THEN Oldpict 

320 St art_over : ON Map GOTO 0 1 dp i ct , Thermal , Therm al 

330 GOTO Oldpict 

340 Oldpict 5 : GOSUB Error 

350 Oldpict: ! THIS SECTION RETRIEVES A BOARD DESCRIPTION OFF A MASS STORAGE DEV 

ICE FOR THE PURPOSE OF DEBUGGING THERMAL 
360 ON ERROR GOTO Oldpict# 

370 0 1 de=P I 

380 Si = " BOA c D DESCRIPTION FROM MASS STORAGE DEVICE" 

390 GOSUB Ragehead 

400 PRINT " You have chosen to input the circuit, board description in T H E R M 
L d i rec t 1 y " 

410 PRINT "from a mass storage device. This program in the T H E P M E L E X System Mi 
1 1 NOT " 

420 PRINT "allow graphical data checking and while faster, there is the c nances 
that the" 



430 


PRINT 


"data is incorrect. If you decide 


t h at 


i t u o u 1 d b e 


b e 1 1 * r t, o check t 


he data" 










440 


PRINT 


"press KO 'Back up) and BOARDS wil 


1 be 


loaded from 


the DEFAULT n. a 


ss s 


t or age 


un it." 








450 


PRINT 










4 60 


PRINT 


" The d at a file c o n t a i n i ng t he 


boar d 


descr i pt i on 


must have been st 


ored 


by " 










470 


PRINT 


"this program. Enter the data til 


e name below, be 


sure to include t. h 


e m a 


s s " 










430 


PRINT 


"storage unit specifier if needed. 




: T14, : FS, : vi 


2 etc ) " 


4 30 


PRINT 


L I N 1 1 1 , " Do NOT use quotes" 









< 1 ) 
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THERML .. THEFHAL ANAL YS I S FOP THEPMELEX 



500 EDIT "UNDER WHAT FILE NR ME IS THE CIRCUIT BORRD FILED Cchange or enter bale 
W) " l N 3. Til c * 

510 IF Nam#**"BACK UP" THEN 530 
520 GOTO 550 
530 Map* 1 

540 LORD "BOARDS" , 1 
550 ASSIGN #1 TO Name* 

560 READ #1,1 

570 READ *1 ; Name*, Pi ct 1 bl *, Case*, Bdl , Bdh,Sca, N • r,Nyr, Nreg, Thi ck_b, Kb. A 1 , K1 , Tem_ 
sol 

580 FOP 1=1 TO Nr eg 

530 READ 4 1 ; 1 1 ype • I J X < I > , Y (I > . T j < I > , Pow < I 1 , R j c < I > 

600 FOR J=1 TO 4 
610 READ 41 ;HKJ, I > 

620 ne::t j 

630 NEXT I ! END OF OLDPICT 

640 ! 

650 GOTO Thermal 

660 Thermals: G 0 S U B E r r o r 

670 Thermal: 1 THIS SECTION IS THE MAIN INPUT AND CALLING ROUTINE 

6oO ON ERROR GOTO Thermals 
690 SS=" THERMAL ANALYSIS OF" 

700 COSUB Pagehead 

710 PRINT T A B v. 3 7 - L EN < P i ♦ 1 b 1 T > * 2 ; CHP* < 1 33 > I P i ct lbl *; CHP* t. 123> ; LIH<2> 

720 PRINT " This section assumes you have completely and correct 1 y described 
the board" 

730 PRINT "itself. You will be asked questions concerning ONLY the environment 

740 PRINT " The first questions are concerned with the cooling air supply. 
Recall the" 

750 PRINT "direction of air flow on the graphics picture is assumed to be left 
to right." 

76 O PRINT "The clearance between the boards is used to determine the velocity o 
f the " 

770 PRINT "cooling air.” ,LIN<2> 

780 Ans * = " 20 " 

790 Map = y 

800 INPUT "ENTER THE INLET TEMPERATURE OF THE COOLING AIR ''deg C ) " , Hns* 

810 IF Ans * = " BACK JJP " THEN Start_over 
820 Tai r =VhL «■ Ans * > 

830 PRINT " INLET AIR TEMP T ai r*" : T ai r ; "deq C" 

340 Tair*Tair +273 ! ALL CALCULATIONS DONE IN ABSOLUTE TE. P 

850 Ans*=" . 00O5" 

360 INPUT "ENTER THE AIR SUPPLY PER BOARD (H A 3/SEC> ", Ansi 
870 IF Ans *=" BACK JJP" THEN Thermal 
880 F ai r*VAL < Ans * ) 

890 PRINT " FLOW FATE OF A I R= " ; F ai r ; " H'WSec " 

900 Ans* = " 15. 24 " 

910 INPUT "ENTER THE DISTANCE FROM THE FACE OF THE BOARD TO THE NEXT OBJECT mm 
>" , Ans* 

920 IF Ans* *" BACK JJP" THEN 780 
930 Zb=VAL < Ans* ) 

940 PRINT " BOARD SPAC ING*" ; Zb; "mm" 

950 2b=Zb+ . OO 1 
960 A n s * = " N U L L " 

970 INPUT "ALL OK >. . PRESS CONT ANY ENTRY FOLLOWED BY CONT WILL ALLOW PEENTP 
Y OF H L L " , H n s * 

980 IF A ns *=" HULL " THEN 1000 



< 2 ) 
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THERHL .. THERMAL ANALYSIS FOR THEPNELEX 



Typically two or three itterations result in a maximum difference 



990 GOTO Thermal 

1300 S$= " CONVERGENCE CRITERIA" 

1910 GOSUB P age he ad 

1020 PRINT " Convergence is indicated by successive itterations 'hit result in 
e 1 ement " 

1030 PRINT "temperatures that differ only by some small amount. Each element te 
nper at ure " 

1040 PRINT "is compared to that obtained in » he previous itteration and if the 1 
argest " 

1050 PRINT "difference is less than a maximum specified error, results are print 

ed. " 

106O PRINT 
1070 PRINT 
on the" 

1030 PRINT "order of one degree centigrade when solving for temper atures and f ou 
r or five" 

1090 PRINT "itterations will result in a maximum difference on the order of .1 W 
at t when " 

1100 PRINT “solving for powers." 

1110 IF T em_so 1 * 1 THEN AnsS=".5" 

1120 IF Tem_so 1 =0 THEN Ansi="l" 

1130 IF T em_so 1 = 1 THEN INPUT "ENTER THE MAXIMUM DIFFERENCE BETWEEN ITTERATIONS ( 
deg C) DEFAULT*. 5" , Ans-f 

1140 IF Tem_soJ=0 THEN INPUT "ENTER THE MAXIMUM PERCENT CHANGE BETWEEN ITTERATIO 

NS (Watt) DEFAULT* 1 \ An s 1 

1150 IF Ans $ = " BACK _UP " THEN Thermal 

1160 Errmax=ABS(VALt Ans* ) > 

1170 GOSUB Database ! THE FOLLOWING LINES FORM THE MAIN CALLING ROUTINE 

1180 GOSUB Calc O 

1190 GOSUB Calcl 

1280 GOSUB C a 1 _ a i r 

1210 GOSUB Calc t 

1220 GOSUB C a 1 c 2 

1230 GOSUB Debug 

1240 GOSUB Solve 

1250 GOSUB Units 

1260 IF Bomb* 1 THEN 1230 

1270 GOSUB Output 

1280 GOSUB What _now 

1290 IF Map=0 THEN GOTO Thermal 

1300 GOTO 01 dpi ct 

1310 ! 

1320 ! END OF THE MAIN CONTROL SECTION OF THE PROGRAM THEPML 
1330 ! 

1340 ! ********* DATABASE* ♦****-*• 

1350 Database: ! IN THIS SECTION MANY OF THE CONSTANTS USED IN THE CALCULATIONS 
1360 ! ARE READ IN FROM THE DATA LINES BELOW 

1370 DISP "WORKING ON NON-CHANGING PARAMETERS" 

1330 PRINT PAGE 

•X-3ECTI0NAL AREA FOP PINS < M 2 ) 

! SURFACE AREA FOR PINS <M«2> 

! VOL COEFF OF EX PA . AIR (1-deg K) AT 300 deg K 
! SPECIFIC HEAT OF AIR 'WATT-SECrKg-deg K ) 

! AVG DISTANCE FROM BOTTOM OF DIF TO BOARD 
! EMMISTVITY OF THE SURFACE OF BOARD 
! EMM I S I V I TY OF THE DIP SURFACE 
(GRAVITY <MrSec A 2) 

(KINEMATIC VISCOSITY AIR <Mt'2rSec> 



1390 Axpin=4.3E-7 
1400 Aspir.= lE-5 
1410 Bet a=3 . 33E-3 
1420 Cpa* 1 . 006E 3 
1430 D i s= . 00 1 
1440 Epsb= . 3 
1450 Epse=.9 
1460 G = 9 . 31 
1470 Gnu* 1 . 634E-5 



( 3 ) 



119 



THSRML 



THEPHAL ANALYSIS FOP THERHELEX 



Hec = . 003 
Hep®. 005 
Ks®59 
Kai r = . 026 
Kpi n=334 
Lpi n®. 0025 
Mu® 1 . 933E-5 
Pr = , 703 
Rho® 1 . 1 774 
S i g=5. 67E-3 
It t=0 

RESTORE 1630 
FOR 1=1 TO 14 
READ W i deset < I ; 

NEXT I 

DATA . 25, . 725. .25, .325, .55, 1.25,0,0, .6, 1, . 13, .35, .43, .72 ! ALL IN INCHES 

FOP 1=1 TO 14 ! THIS READS THE CASE LENGTHS C 0P EACH TYPE 

READ Lsnsft < I > 

NEXT I 

DATA .725, .25, .325, .25, 1, .55,0,0, 1, .6, . 13, .35, .43, .72 ! ALL IN INCHES 

MAT Widese*. = >. 25. 4 >-Wi dsset ! CONVERT TO mm 

MAT Lenset=(25. 4'^Lenset 

FOR I®1 TO 14 ! THIS READS THE # OF PINS FOR EACH TYPE 

READ Np i n ( I) 

NEXT I 

DATA 14, 14, 16, 16,24,24.0,0,20,20, 16,24,40,64 

FOP 1=0 TO 14 ! THIS ^EADS THE CASE LABELS FOR EACH TYPE 

READ T^pel(I) 

NEXT I 

DATA EMPTY, DIP 14, DIP 14, DIP 16, DIP 16, DIP 24. DIP 24 , NULL , NULL , D I P 40, DIP 4 
0, FLAT16,FLAT24, FLAT40, FLAT64 
1730 FOR 1=1 TO N reg 
1790 IF Itype<I)*0 THEN 1330 
18O0 Le(I) = L?n-et., Itvpe<I>> 

1810 Ue< I > = Wi deset < Ityp*< I > > 

1820 GOTO 1S40 
1830 We< I >=Le< I > = 0 
1840 NEXT I 

1850 IF Tern sol=l THEN MAT Te = <300-> 1 INITIAL GUESS FOR TEMP CASE = 27 deg C 

I860 IF Ter., sol =0 THEN MAT ?ow='.25> ! INITIAL GUESS FOR POWER = .25 WatU 

1870 MAT Le=< . 001 >*Le 

1830 MAT Ue-C.001)4Ue 

1830 MAT X*C.001>*X 

1900 MAT Y = < . 00 1 ) *‘i 

1910 MAT Ae=Le.We 

1920 MAT He = ( He ■: > ! ASSUME ALL CERAMIC MAKE CORRECTIONS AS NEEDED IN Oil C0 

1930 Bdl«Bdl*.031 

1940 Bdh=Bdh* . 00 l 

1950 Areg=Bd1 +Bdh/Nreg 

I960 Hr=Bdh'Nyr 

1970 Lr=Bdl/Nxr 

1930 Thick b=Thick_b*. 001 

1990 A 1 =A 1 1 £-6 

2000 RETURN ! END OF DATABASE 

2010 ! 

2020 Cal 0e : GOSUB Error- 

2030 Cal c O : ! ON ERROR GO TO CflLCOE 

2040 Pout ot =N i>->g = He avg=W* a‘;g=0 



! HEIGHT OF CERAMIC 
! HEIGHT OF PLASTIC 



DIP PACKAGE *♦« 
DIP PACKAGE 



! THERMAL COND OF STEEL 
•THERMAL COND OF AIP 
! THERMAL COND OF PINS 
! AVG LENGHT OF PINS 
! DYNAMIC VISCOSITY 300 
I PRANDTL NUMBER 
(DENSITY OF AIR AT 30Q 
•STEFFAN-BOLTZMAN 



< WATTS-' M-deg ID 
< WATTSxM-deg K> 
< WATTS^M-deg K> 
< M ; 

< Kg^ M-3ec ) 

< AT SOOdeg K) 

<. Kgs M^S > 

(. W/'M /s 2-deg K^4) 



300 deg K 



1 THIS READS THE CASE WIDTHS FOR EACH TYPE 



1480 

1490 

1500 

1510 

1520 

1530 

1540 

1550 

1560 

1570 

1530 

1590 

1600 

1610 

1620 

1630 

1640 

1650 

1660 

1670 

1630 

1690 

1700 

1710 

1720 

1730 

1740 

1750 

1760 

1770 
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THERMAL ANALYSIS FOP THERMELEX 



2050 FOR 1*1 JO Hreg 

2060 IF Case* C I , I ]*•• l" THEN He<I>=Hep 
2070 IF Ityp«<I>«0 THEN 2120 
2030 Pout ot = Pout ot +Poui, I 

2030 U«aug*U«avg+<:u«< IW. 001 Nr*g ! RVG WIDTH” INCLUDING PINS 
2100 H«avg*H®avg+<H#< I > + Lpl nVHrsg ! AVG HEIGHT INCLUDING PINS 

2110 Navg = Navg* l/'Hyr ! HUMBER OF ELEMENTS IN AN AVG CROSS SECTION 

2120 NEXT I 



! ROUGHNESS FACTOR 



TURBULENT AT 1000 DUE 
TO MANY TRIPS OF COMPS 
LAMINAR FOR iGz>.05 



! END OF CRLCO 



2130 Aai r =Bdh*Zb-N aug- We aug*H« aug 
2140 Peri m=2»Edh+2 *N avg-^He aug + 2-^Zb 
2150 Fr=l + 5*<F«rim-2-*Bdh-2-*-Zb)'P*rifri 
2160 Dh=4*Aai r/Feri m 
2170 Vai r = Fai r/ Aai r 
2130 R®*Vai r*Dh '"Gn*.4 

2130 IF R e- > 1 0 0 O THEN Hbf =. 023*Kai r/Dh»R®''. 8 
2200 ! 

2210 IF Re< = 10OO THEN Hbf = 5 . 40-*Kai r* Fr-'Bh 
2220 ! 

2230 Rl_hor=Lr/<Kl *ABS<A1 >> 

2240 R I ”ver»Hr/ (Kl+ABS^Al >) 

2250 Rb~hor=Lr. / ■' Kb*Hr*Th i ctc_b) 

2260 Rb~v«r = Hr/(Kb*Lr-<-Thi ck_b) 

2270 RETURN 
2230 ! 

2230 Cal cl*: GOSUB Error 

2300 Cal c 1 : ! THIS SECTION COMPUTES SOME OF THE NON-CHANGING PARAMETERS 

2310 ON ERROR GOTO Calci-e 
2320 FOR 1=1 TO Nreg 

2330 1 BELOW HERE WE CALCULATE THE CONVECTIVE LOSSES FOP EACH BOARD REGION 
2340 IF X<I)/Dh>10 THEN 2330 ! OUTSIDE OF THE DEVELOPMENT REGION 

2350 Gz=Re + Pr*D*v'X< I > 

2360 Hb=.664*Kair,'<l. l*Dh'»-3QP'.Gz*< l+7.3*SQfi<Pr'Gz>>'Pr>*-Fr ! Eq 13.43 KNUTZEN 

4 KATZ LAMINAR IN DEVELOPMENT REGION 

2370 IF 1/G=>.05 THEN 2330 

2330 GOTO 2400 

2330 Hb=Hbf 

2400 Aregt ot = 2*Rreg-Ae < I • 

2410 Rb_c ony( I > = 1 ■ ' < Hb-Ar-s gr ot > 

2420 ! NOW WE GET TO THE ELEMENTS ON THE BOARD 

2430 IF Itype<I>=0 THEN 2640 ! IF NO ELEMENT THEN SET VERY HIGH RESISTANCE 

2440 IF 1 1 yp* < I ) < 1 1 THEN 2430 

2450 Dis=. 1*DU 

2468 Lp i n= . 1 *Lp i fi 

2470 Asp i n= . I -Asp i n 

2430 Rpi ns_cond=Lpi rv'CKpi mrR.x.pi rt+Npi n( Itype < I ) ) ) 

24 30 Rgap_c ond=D i s^< Kai r ( I ) ) 

2500 IF ItypedXll THEN 2540 
2510 D i s= 1 0*D 1 s 
2520 Lpin=10*Lpin 
2530 Asp i n= 1 0-Rsp 1 n 



2540 Re_cond< I >=Rpi ns c ond*Rg ap_c ond/ < Rp i ns cond+Rgap cond) 
2550 Abl as=Ae< I >+2*L*“l )-H«= + Aspi n«Mpin< Ityp*' - I > > 

2560 Heb 1 as =Ht 

2570 IF I < 6 THEN Hb^ I .*»Hb 

2530 Rebl as=l/(Hebl as*Rbl as> 

2530 Hest ag = . 57-Kai r«*Pr* . 4*SQR < V ai r/(H-: < I >*Gnu; > 

2600 fistag=2*We< I ) *He < I ; 

2610 R«stag=l/<H*it ag-Rst ag ) 
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2620 

2630 

2640 

2650 

2650 

2670 

2680 

2680 

2700 

2710 

2720 

2730 

2740 

2750 

2760 

2770 

2780 

2790 

2800 

2810 

2820 

2830 

2840 

2850 

2860 

2870 

2830 

2890 
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3190 
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Re cony<I) = Peb1as*Restag''<Peb1as + Restaa) 

GOTO 2660 
Re_corw< I )=1E30 
Re condc I )= 1E30 

! BELOW HERE CALCULATE THE BOARD CONDUCT I VE RESISTANCE 
J=l ! HERE IS BOTTOM OF REGION 
IF I *• N x r > = N r e g THEN 2720 
IF NK J, I >«0 THEN 2740 

R1 <J, I >*R1 vsr/NI < J, I)^Rb ver/(R1 wer/HKJ, I)+Rb ver) 

GOTO 2750 

R 1 < J , I ) = ft b_v e r '2 E - 4 O 
GOTO 2750 
R 1 < J , I ) *Rb_ver 

J=2 ! HERE IS RIGHT SIDE OF REGION 

IF I MOD Nxr*0 THEN 28O0 
IF N 1 < J , I i=0 THEN 2320 

R 1 < J , I)«R1 hor/Hl <J, I ) *Rb_hor < R 1 _hor/N 1 <J, I)+Rb_hor> 

GOTO 2830 

R 1 ( J , I >=Pb_hor/2E-40 
GOTO 2830 
R 1 ( J , I ) = P'b_hor 

J = 3 ) HERE IS THE TOP OF THE REGION 

IF I < = Nxr THEN 2330 
IF NK J, I ;=0 THEN 2900 

R 1 < J , I ) = R 1 ver- NKJ, I wRb **er^ R1_o*r/Nl <J, I>+Rb_ver> 

GOTO 2910 

RKJ, I ) =Rb ver/2E-40 
GOTO 2910 
R K J , I ) =Rb_uer 

J = 4 ! HERE IS THE LEFT SIDE OF THE REGION 

IF <I-1) MOD N ^ r = O THEN 2960 
IF NKJ,I)=0 THEN 2930 

R1(J,I)=R1 hor^'N 1 < -J , I > »Pb hor/CRl hor^Nl < J, I ) +Pb_hor > 

GOTO 3000 

R 1 C J, I >»Pb_hor/2E-40 
GOTO 3O00 
R 1 (J, I ; = Rb ho r 
GOTO 3010 

Rb_c onu <I) = lKHfc*Aregtot) 

NEXT I 

RETURN ! END OF CflLCl 

i 

Calc teZGOSUB Error 

Ca1c”tZ! THIS SECTION DETERMINES SURFACE TEMP FROM JUNCTION TEMP AND Rj_c 
ON ERROR GOTO Calc te 
IF Tem_S 0 l=l THEN RETURN 
FOR 1=1 TO Nr eg 
Te< I ) = Tj< I )-Pow( I>*R j_c< I ) 

IF T*<IXTair<I>vi. l THEN Te < I > -Tat r < I >■* l . I 
NEXT I 

RETURN ! END OF Ca1c_t 

Cal c2e ! GOSUB Error 

C a 1 c 2 *. ! THIS SECTION CONTAINS THOSE PARAMETERS WHICH CHANGE WITH TEMPS 

ON ERROR GOTO Ca1c2e 

DISP "WORKING ON CHANGING PARAMETERS" 

FOR 1=1 TO Nr eg 

IF 1 1 y p e < I > < / 0 THEN 3230 



122 



THERML .. THERMRL ANALYSIS FOR THERMELEX 



3200 Rt op r < I > = l E 5 0 
3210 Rgaprad=lE50 
3220 GQT0""*3250 

3230 Rgap_r ad = f Epi ■f-Epsb-Epae-'-Epab > / < 4*3 i g*Epse*Epab*R's<I^*T£<I>''-3> 

3240 Rt op r < I ■ = < < 1-Epsb>' < Ep3b~Rreg) + 10/Rreg+ < 1-Epse > '< R« ■. I >*Epse > ) / <4*Si g*Te< I ) 
-3) 

3250 Rtote_b< I ) =Rgap_r ad*Ps_c ond < I )/ <Rgap_rad+Re_cord< I > > 

3260 ReU<I)=Re cony(I)tStote_bU) !!!! MUST BE SUM FOR RLGEERR I C PERSONS 

3270 NEXT I 

3230 RETURN ! END OF CRLC2 

3230 ! 

3300 Calc ai re: GOSUB Error 

3310 Calc~air:' THIS SECTION FINDS THE RIR TEMPS FOR ERCH REGION BRSED ON POWER 
3320 ON ERROR GOTO Calc_air* 

3330 C f r=Cp a* r ii r - P h o 
3340 Pout ot =0 
3350 FOR 1=1 TO Nreg 

3360 Poutot =Powtot *Pow< I > 

3370 Tair<I>=Tair+.5*Pou<I>*Nyr/Cfr 

3330 IF <1-1 > MOD Nxr=0 THEN 3420 

3390 FOR K= I - 1 TO I-<I-1> MOD Nxr STEP -1 

3400 Tai r v I > = Tai r c I > +Pow< k) *Nyr/Cf r 

3410 NEXT K 

3420 NEXT I 

3430 Tout =Tai r+Poutot /Cf r 

3440 RETURN ! END OF CRLC_RIR 

3458 ■ 

3460 So 1 yee : GOSUB Error 

3470 Solve: ! THIS SECTION SOLVES THE PROBLEM 

3480 ON EFPOP GOTO Solve* 

3490 DISP "WORKING ON SETTING UP THE MRTPIX" 

3500 Itt=Itt+l 

3510 BEEP 

3520 WRIT 300 

3530 BEEP 

3540 GOSUB 3*t_up 

3550 GOSUB Elu 

3560 Err=0 

3570 Tm ax = 0 

3530 Pmin=100O 

3590 IF T fi»_so 1 = 1 THEN 3730 

3600 FOR 1=1 TO Nreg 

3610 Pneu=(T?( I >-B<I ) >/Rtote_b< I > + < T e < I >-T ai r < I ) >/R*_c on-; < O + < Te < I ) - B < I ) >/Rt 

op_r < I ) 

3620 IF Pneu<0 THEN Pneu=Pow< I >/2 

3630 IF PnewKPmin THEN Frmri=Fneu 

3640 IF ftBS< <Pow< I>-Pnew)/Pow<I> )>Err THEN Err = RBS < < Pow < I )-Pneu) / Pou < I > > 

3650 Pou < I > =Pneu 

3660 NEXT I 

3670 Tro ax = 500 
3630 MRT Tb-B 

3690 IF Err.Errmi 1O0 THEN 3340 

3700 GOSUB Calc t 

3710 GOSUB Calc”air 

3720 GOTO 3920 

3730 FOR 1=1 TO Nreg 

374 0 Tneu = < Pou < I ■ *Rt ot e_b* I ) *Pe_c onv < I > + B< I '> 4Pe_c©r>" < I > + T ai r < I > -Pt . ot. e_b < I > ) 

/Rele(I) 



< 7 > 
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3750 IF ABS<Te<I>-Tnew>'Err THEN Err= ABS < Te < I > -Tnew > 

3760 Te < I > = T new 

3770 Tj<I>=Te< I) + Pow( I > •» R j c < I) 

3730 IF Tj<I)>Tmax THEN Tmax»Tj<I> 

3790 NEXT I 

3800 Pm l n=0 
3810 NAT Tb = B 
3820 MAT A = ZER 

3830 IF Err>Errma;f THEN 3920 

3840 DISP 

3850 FOR C=1 TO 4 

3860 BEEP 

3870 WRIT 90 

3830 BEEP 

3890 WRIT 150 

3900 NEXT C 

3910 RETURN 

3920 GOSUB Temp_print 

3930 IF < 1 1 1 >30-' OR Err >20*0 ) OR <Tout>400> THEN GOTO Bomb 
3940 GOSUB C al c 2 

3950 GOTO Solve ! END OF SOLVE 

3960 ! 

3970 Set_up: ! THIS SECTION SETS UP THE Nreg 3IP1UL EQUATIONS IS MATRIX FORM 
3980 FOR 1=1 TO Nreg 

3990 IF Tem_io 1 =0 THEN B ■ I > =Te <1 > /Pt ot e_b < I ) + T ai r < I > xRb_c onv < 1 ) 

4OO0 IF Ism sol =1 THEN B I ) = < Pow ( I > +R>?_c ortv < I ) +T ai r < I > > ✓ Re 1 e < I ) + T ai r < I ) /Rb^c onu < 
I) 

4010 Ledge=Re dge = Tedgs = Bedge = 1 
4020 IF (1-1) MOD N x r = 0 THEN Ledge=2 
4030 IF I MOD Nxr = 0 THEN Pedge=2 
4040 IF K*Nxr- THEN Tedge=2 
4050 IF I+Nxr=Nreg THEN Bedg*=2 
4060 IF <I-1) MOD N- r=0 THEN 4030 
4070 fl<I, I-1> = -1/R1 i'4, I) 



4030 IF I MOD N x r = O THEN 4100 
4090 R< I, If 1 > = -l PK2, I ) 

4100 IF I < = N v r THEN 4120 
4110 A< I, I-Nvr > = - 1 / P 1 * 3 , I > 

4120 IF I +Nxr >=Nreg THEN 4140 
4 130 fl< I , I+Nxr )=-l PI a. I) 

4140 R< I , I )=<Bedge=l >/Pl < 1 , I >* < Redge= 1 > /R 1 <2, I > + <T*dg*»l > - P 1 < 3 , n + < Ledge* 1 > /P 1 < 4 



, I) + l/'R»ct e_b ‘ I ) + 1 • Rb c or»v < I > - < T ern_so 1 = 1 > *Pe_-: omn I < Re 1 e < I ) *Rt ot e_b I > > 

4150 NEXT I 

4160 RETURN ! END OF SETJJP 

4170 • 

4130 E1'.j: ! THIS SECTION PERFORMS R LU DECOMPOSITION OF THE 'fi'MRTRIX 

4190 DISP "WORKING ON ITERATION NUMBER ";Itt 

4200 Nm 1 = Nr e g- 1 

4210 FOR K= 1 TO Nml 

4220 Kp 1 =K+ 1 

4230 FOR I =Kp 1 TO Nreg 

4240 G*-A<I,K>'A<K,K> 

4250 A < I , K > =G 

4260 FOR J=Kpl TO Nreg 

4270 A<I, J)=A<I,J)+G»A<K,J> 

4280 NEXT J 

4290 NEXT I 

4300 NEXT K 



< 8 ) 
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4319 Solver: ' THIS SECTION SOLVES THE NEW MATRIX AND PLACES THE ANSWERS 

4320 Np 1 =Nr eg* 1 ' INTO MATRIX B TO BE PASSED BACK TO SOLVE 

4330 FOR K = 1 TO Nml 

4340 Kp 1 = K + 1 

4350 FOR I=Kpl TO Nreg 

4360 B<I>*B<I>+A<I,IO-*B<K> 

4370 NEXT I 

4330 NEXT V 

4390 B<Nreg>=B<Nreg' ' A < Nreg , Nreg ' 

4400 FOR K=2 TO Nreg 

4410 I =Np 1 -K 

4420 J 1 * I ♦ 1 

4430 FOR J=J1 TO Nreg 

4440 bu)=B(I;-A'I,j!'»bj) 

4450 NEXT J 

4460 B< I )=B< I >'A< I , I > 

4470 NEXT K 

4430 RETURN ! END OF SOLVER 

4490 T emp_pr i nt : ! THIS SECTION USED FOR INTERMEDIATE OUTPUT 
4500 EXIT GRAPHICS 

4510 IF <T emprt = 1 ) AND fRns=PI> THEN PRINTER IS 8 

4520 PRINT SPA< 10' , "DATA FOR " J P i c t 1 b 1 J S " # " ; 1 1 1 J " ITTERflT ION" 

4530 PRINT LIN<2) 

4540 FIXED 2 

4550 PRINT " PEG 4 Tease 'DeqC' Tjunc <DegC> Treg CBeqC) Pou 

<U) R 2 -C >U/C' ■* 

4560 FOR 1=1 TO Nreg 

4 570 PRINT I J TABk 12) , Te< I ) -273 ; TAB < 26 > . T j < I ) -273 ; TAB C 4 1 ,' , Tb < I > -273 ; TAB < 5 2 > , 

Pou< I ) ; TAE<62) , P j_c < I > 

4580 NEXT I 

4590 PRINT LIN<2> 

4600 PRINT "BOARD THICKNESS*" I 1 0 0 O e T h i c k _b ; "mm AND CONDUCTIVITY =" ; Kb J "Wat 

t s/M-K" 

4610 PRINT 

4620 GOTO 4580 ! PRINT " I RB_CONV RE_C0NV TE Hb" 

4630 FOR 1=1 TO 5 

4640 PRINT USING 4660 ; I , F b_c orw < I ) , R*_c onv U > . Te a > -273 , Hb < I ) 

4650 NEXT I 

4660 IMAGE DD, 4X, 4 ■ 4D. 4D, 4X) 

4670 PRINT 
4630 FIXED 4 

4690 PRINT "COOLING AIP FLOW OF " : F a i r ; " M 3 per SEC VEL= " ; V-ai r ! " M. Sec 
Vai r*39^3; “ FT/S > " , L I N < 1 > 

4700 PRINT "INLET AIR TEMP= " ; T ai r-273 ; " deg C OUTLET AIR TEMF=" ; Tout-273; "de 
g C" 

4710 PRINT LIN<2>. "LARGEST DIFFERENCE BETWEEN ITTE RATIO NS = " ; Err ; " * * — * " 

4720 PRINT " 

“ , L I Mi 2 ' 

4730 PRINTER IS 15 

4740 STANDARD 

4750 RETURN 

4760 Debug: flrt;l="N"! TEMP DEBUG FOR RESISTANCES 
4770 RETURN ! 

4730 INPUT "DO YOU WISH TO HAVE A LIST OF ALL THE RESISTANCES PRINTED OJO or YES 
>?" , A ns* 

4790 IF UPCKAnsI C 1 , 13 > = "N" THEN RETURN 
4300 PRINTER IS 8 
4310 FIXED 5 
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4820 PRINT " PEG# RE CCNH RE CONV PTOTE B RTOP R RB 

CONV " 

4330 FOR 1*1 TO Nreg 

4340 ! PRINT TfiBt2> . I , TAB ( 4> , Re_cond< I > , TflB<4> ; R* c onv< I) , TAB C 4 ) , Rt. ot e_b ( I) ; TAB < 
5);Rtop_r(I);TftB(4'.;Rflea) 

4350 PRINT TAEC 2 ) , I J P<? condtDjP? conv(I);Rto»? b<I>;Pt-op r<I>;Rb conv( I) 

4330 NEXT I 
4870 PRINTER IS 13 
4330 STANDARD 
4330 RETURN 

4900 Units!! THIS SECTION CONVERTS FROM METER TO mm 

4910 Bd 1 = 1 000-Ed 1 ; HERE CORRECT UNITS FOP USE IN GRAPHICS AND TO REWORK ANALYSIS 

4920 Bdh* 1 000-Bdh 

4930 MAT X=<10O0^X 

4940 MAT Y = < 1 0 0 0 ) -* Y 

4950 MAT L**<1O00>*L# 

4 960 MAT We*<100o>*We 
4970 Thi ck_b*1000*Thi ck_b 
4930 A 1 =A 1 * 1 Eo 

4990 RETURN ! END OF Units 

5000 • 

5010 Outputs : GOSUB Error 

5020 Output:! THIS SECTION OUTPUTS TO GRAPHICS ON A BLANK BOARD 
5030 ON ERROR GOTO Outputs 
5040 PLOTTER IS "GRAPHICS" 

5050 GRAPHICS 
5060 MSCALE 0,10 
5070 CSI2E 2 
5030 LORG 5 
5090 MOVE 13, 130 

5100 LABEL " sir flow > " 

5110 CSIZE 3 

5120 Bdh=BdhxSc a 

5130 Bd 1 *Bd 1 /'Sc a 

5140 IF Sc a* 1 THEN 5220 

5150 LORG 5 

5160 LDIR PI/2 

5170 MOVE 175,70 

5130 IF Sc a> 1 THEN LABEL "THIS PICTURE IS 1 ' " &VALf (Sea '• Z< " SIZE" 

5190 IF Sc a< 1 THEN LABEL "THIS PICTURE IS 2X SIZE" 

5200 LORG 5 
5210 LDIR 0 

5220 MOVE 1 00-LEhK "OUTPUT DATA FOP " &P i c t 1 b 1 S > ^2 , 1 35 
5230 LABEL USING "K"; "OUTPUT DATA FOR "tPIctlbl* 

5240 MOVE 0,0 
5250 DRAW O , Bdh 
5260 DRAW Bd 1 , Bdh 
5270 DRAW Bd 1 , 8 
5230 DRAW 0,0 
5290 LINE TYPE 3 

5300 Lr=Bdl/Nxr ! LR = LENGTH OF EACH REGION 

5310 Hr=Bdh/Nur ! HR = HEIGHT OF EACH REGION 

5320 FOR 1=1 TO Nxr-1 

5330 MOVE I *Lr , O 

5340 DRAW I *■ L r , B d h 

5350 NEXT I 

5360 FOR 1=1 TO Nyr-1 

5379 MOVE O, I *Hr 



< 10 ) 
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5380 DRAW Bdl , 
5390 NEXT I 
5400 LINE TYPE 1 
5410 LORG 8 
5420 Nr =0 
5430 CSIZE 2.4 



5440 


FOR 1*1 TO Nreg 








5450 


MOVE X < I> - . 


35*Lr,Y(I)+.35+Hr 






5460 


LABEL USING 


"K" ; VALi( I > 






5470 


MOVE X < I ) + . 3 L r , Y < I > + 


25*Hr 






5480 


LABEL USING 


" K " ; T m p e t 


1 1 y p e v I > > 






5490 


IF 1 1 y p <5 < I > 


s 0 THEN 5570 






5500 


IF Tj<I>-27: 


3< . 95* (To, a 


<-273) THEN LABEL IJ'i 


SING 5610; T i ( I )-273 


5510 


IF Tj(I)-27 


*. 95*(Tfi..ax-273> THEN LAEEL USING 5600 J T j < I ) -273 


5520 


IF P oij < I ) =1.05*Pmin THEN LABEL 


USING 56] 


:-0 ; Pom (15 


5530 


IF Pow(I)>l 


0 5 * P m i n THEN LABEL USING 5620 ; Row ( I > 


5540 


LABEL USING 


5610; Tet I) -273 






5550 


GOTO 5590 










5560 


! 










5570 


LAEEL USING 


" K " ; " " 








5580 


LABEL USING 


5610; Tb< I ) -273 






5590 


NEXT I 










5600 


IMAGE D D D . 


, " c " 








5610 


IMAGE DDD.D," C* 










5620 


IMAGE D.DDD," W' 










5630 


IMAGE D . D D D , " W" 








5640 


FIXED 4 










5650 


DUMP GRAPHICS 










5660 


PRINTER IS O 










5670 


PRINT SPA < 20 > , " 










5630 


PRINT S P A < 2 O ; , " 








" 


5690 


PPINT 6PAC20 ■ , " 


TYPE 




EMPTY 


" 


5700 


PRINT SPA (20)," 


T j unc 


or 




" 


5710 


PRINT SP A ( 20 > , " 


POW 




Tboard 


" 


5720 


PRINT SP A ( 20 ) , " 


Tc ase 






" 


5730 


PRINT SPA(20 ) , " 








" 


5740 


PRINT L I N < 2 ; 










5750 


PRINT SPA ( 14; ; ' 


1 FLOW PATE VELOCITY 


Tin Tout " 


5760 


PPINT USING 


5300; "COClInG AIR "; 


! Fai r; “ M ~ 3 s S e •: " ; Vai r; "MxSec " ; T ai r-273; " 


deg C 


Tout -273 deg < 










5770 


STANDAPD 










5780 


IF Tern *ol=0 THEN PPINT USING 58 1 0 LARGEST CHANGE IN POWER BETWEEN ITTERAT 


IONS 


#“ ; 1 1 t-l ; " AND ; l» t ; " = ' 


' ; Err 






5790 


IF Ten, so1=l THEN PPINT USING 5310; "LARGEST CHANGE IN TEMP BETWEEN ITTERATI 


0N3 #" ; Itt-1 ; " AND 4" ; 


; 1 1 1 ; " = " ; 


i Err 






5800 


IMAGE 1 1A, 2X, . 40 , 7 A , 2X , DD. 


, 4D, 5A, 2X. 3D. D, 5A. 2: 


s, 3D. D, 5A 


5810 


IMAGE < , /’45A , DD , r 


j A , DD , A , DD . 


, 4D 






5820 


! MOPE OF THE CIRCUIT DESCRIPTION CAN BE PREi 


iENTED here 


5830 


PRINT LIN< 1) , "CIRCUIT BOARD DESCRIPTION IS STORED UNDER " J CHPf < 132 > : N am** ; 


CHR$< 123) 










5840 


PRINT " 












" , L I N ( 2 ) 










5350 


PRINTER IS 16 










5360 


EXIT GRAPHICS' 










5370 


RETURN 








! END OF OUTPUT 



5380 ! 

5890 What_noue .* GOSUB Error 

5900 What now: ! THIS SECTION PRESENTS 



THE VARIOUS OPTIONS AVAILABLE AND DIRECTS 
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5910 ! PROGRAM CONTROL AS NEEDED 

5920 ON ERROR GOTO Mhat_nowe 
5930 S*="WHAT NOW ?" 

5940 GOSUB Pagehead 

5950 PRINT TAB< 15 ) , 11 You ha."* completed the thermal analysis of 
5960 PRINT TAB ( 37-LEN v P i c t 1 b 1 * > ; C H R * ( 129 > : Pi a 1 bl * ; CHR* (123) 

5970 PRINT TAB < 1 7 « , “ t. he options available are listed below?" 

5980 PRINT LIN< 1 ) , TAB' 5) , " 1 . CHANGE COOLING AIR PARAMETERS AND PE-hNALYZE ";CHR 
1(132); Name*; CHR* k 123) J “ . “ 

5990 PRINT LIN< 1 > , TAB 1 5) , “2. SENSITIVITY ANALYSIS FOR CHANGES IN COOLING AIR PA 
RRMETERS OF " ; CHR* (132) ; Names ; CHP*< 123)1 " . " 

6000 PRINT LIN< 1 ),TAB<5> . "3. MAKE CHANGES TO BOARD DESCRIPTION IN FILE NAME " J C 
HR*< 132) ; Name* ; CHR* ■ 128) ; " . " 

6010 PRINT LlhKl), TAB'S), "4. °ETP I EVE A NEW BOARD DESCRIPTION FROM MASS STORAGE 



6020 PRINT LIN< 1 ) , TAB<5) , "5. INPUT A NEW BOARD DESCRIPTION FROM THE KEYBOARD." 
6030 PRINT LlN< 1 ) , TAB (5) , "6. TERMINATE SESSION" 

6040 Map=Bomb=0 
6050 Ans* = ,1 l" 

6060 INPUT "YOUP CHOICE FROM ABOVE <1 , 2 , 3 f 4 , 5 . 6 > ? " , Ana* 

6070 IF Ans * = " BACK-UP " THEN RETURN 
6080 Ans= I NT < VAL ( Ana * > ) 

6090 IF < A n s > 0 ) AND ’.Ans<=6) THEN 6120 
6100 GOSUB Errin 
6110 GOTO What_now 

6120 ON Ans GOSUB Pe do , 3e ns i , Ch ange , Get _new , K ev_ne w , Term i n at e 
6130 IF Ar.s=l THEN RETURN 
6140 GOTO What_now 

6150 Redo*. RETUP N ! THIS WILL ALLOW RESTART OF THIS PROGRAM WITH SAME BOARD 
6160 ! 

6170 Sens i e : GOSUB Error 

6180 Sen-si :» THIS SECTION PRODUCES PLOTS OF OUTPUT VS COOLING AIR PARAMETERS 
6190 ON ERROR GOTO Sensie 
6200 3-1 = " SENS I T I VTY ANALYSIS" 

6210 GOSUB Page he ad 

6220 PRINT " This section allows you to investigate the effects of variations 
in the air flow rate. " ; 

6230 IF Tem_sol=l THEN PRINT "Plots of Maximum Junction Temperature vs Flow Rate 
of the air are produced" 

6240 IF T*m_sol=0 THEN PRINT "Plots of Minumum Power vs Flow rate of the air are 
produced. " 

6250 PRINT " You specify the maximum flow rate per board <M A 3rsec) and five s 
e par ate M 

6260 PRINT "analyses are performed and the results plotted. NOTE: selecting a rn 
ax i mum flow" 

6270 PRINT "that is evenly divisible by five (5,10,30) will result in better loo 
king axes. ",LIN(1> 

6280 PRINT "These plots may be produced on either the screen (with hard copy via 
Key3) or " 

6290 PRINT "on a peripheral plotter such as the Hp 9872A. In addition a printed 
t abu 1 ar " 

5300 PRINT "output of the results at each of the airflow rates may be produced." 

6310 PRINT 

6320 Ans*= " NULL “ 

6330 INPUT "WHAT IS THE UPPER LIMIT ON THE AIR FLOW RATE FOR THE SENSITIVITY ANA 
LYSIS’", Ans* 

6349 IF Ans* = " BhCK_UP " THEN What_r.ow 

6350 IF Ans* ="NULL " THEN 6328 
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6360 Ans=ABS C VAL < Rns-$ ) > 

6370 Xmax=Ans 
6380 Xm i n=0 

6390 Xstep=<Xniax-Xmin>^5 
6400 Ymin=0 
6410 Mult=10000O 

6420 IF Xmax >= . 00 i THEM Mult-10000 
6430 IF Xmax > = . 0 l THEM MultMOOO 
6440 IF Xmax >= . I THEM Mult-100 
6450 AnsS=“N“ 

6460 INPUT "BO YOU WISH TO HAVE TABULAR RESULTS FOR EACH OF THE AIR FLOWS? <N or 
Y> “ , Ansf 

6470 IF Ans-f = "BAC K_U P “ THEN Sensi 
64S0 Temprt=0 

6490 IF UPC*<An*J’‘- ,, Y" THEN Temprt=l 
6500 GOSUB Dat iDai« 

6510 Map=i 

6520 IF Tem_sol=l THEN 6550 

6530 FOR Fair»X»ax TO Xstep STEP -Xstep 

6540 GOTO 6560 

6550 FOR Fair* Xstep TO Xmax STEP Xstep 

6560 GOSUB C a 1 c 0 

6570 GOSUB Cal cl 

6530 GOSUB Cal c_t 

6590 GOSUB Cal c 2 

6600 GOSUB Calc air 

6610 GOSUB Solve 

6620 IF T e m_io 1 = 1 THEN Y=Tmax-273 
6630 IF T ?m_i o 1 = 0 THEN Y*Pfnin 
6640 IF Tempr t = 1 THEN PRINTER IS 0 
6650 IF Temprt-1 THEN GOSUB Temp_pr i nt 
6660 IF Map = 1 THEN GOSUB Plot. 

6670 IF Aris $ = “ BACKJJP" THEN 6450 

6630 Ans* It t -Map-0 

6690 GRAPHICS 

6700 PEN 1 

6710 MOVE Fair.Y 

6720 LORG 5 

6730 LABEL USING " k " ; " * " 

6740 PEN u 
6750 NEXT Fair 
6760 EXIT GRAPHICS 
6770 DUMP GRAPHICS 
6780 PRINTER IS 0 
6790 PRINT LINC2) 

6300 PRINTER IS 16 
6810 GOSUB Units 

6820 GOTO What_now ! END OF SEHSI 

6830 ! 

6340 Change:! THIS SECTION LOADS BOARDS TO ALLOW CHANGES TO THE CIRCUIT DESCPIPT 
6850 M ap = 1 
6360 PRINT FACE 

6370 DISP "WORKING LOADING BOARDS " 

6330 LOAD "BOARDS", 1 

6890 Cet_new: ! THIS SECTION EXPLAINS THE TWO WAYS TO GET A NEW BOARD DESCRIPTION 
6900 Si*"" INPUT OF NEW BOARD DESCRIPTION" 

6910 GOSUB Pagshead 

6920 PRINT " The new board description may be read in from mass storage in tg 
o different" 
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6930 PRINT "programs. THERML (the program in core now) will not allow visual 
checking or" 

6940 PRINT "modification of descriptive data for the circuit board. The data is 
read in'* 

6950 PRINT "at your direction but only questions concerning the environment exte 
rior to the" 



6960 PRINT "circuit board are asked. The second method of input from mass stora 
ge involes" 

6970 PRINT "the program BOARDS which allows both data checking through graphic 
s and " 

6990 PRINT "opportunities to modify the descriptive data." 

6930 PRINT 

7000 PRINT LIN< 1 ) , TAB (5) , " 1 . READ A NEW BOARD DESCRIPTION USING THERML . (no c 
orr ec t i ons ) “ 



7010 PRINT L IN < 1 ) , TAB 5 > , "2. FEAD A NEW BOARD DESCRIPTION USING BOAPDS . <allo 
ws corrections)” 

7020 PRINT L I N <. 2 ) , " REMEMBER THE THERMELEX SYSTEM MUST BE IN THE DEFAULT MASS ST 
ORAGE DEVICE" 

7030 Ansi® " 1 " 

7040 INPUT "YOUR CHOICE FROM ABOVE < 1 , 2)?" , Ansi 
7050 IF Ansi® "BACK-UP " THEN What_now 
7060 Ans=VAL( Ansi) 

7070 IF < Ans >0 ) AND <Ans:3:- THEN 7100 

7080 GOSUB Errm 

7090 GOTO Get_new 

7100 IF Ans = 2 THEN 7130 

7110 Map® 1 

7120 RETURN 

7130 ! 

7140 Map =2 

7150 DISP "WORKING LOADING BOARDS" 

7160 LOAD "BOAPDS" , 1 
7170 ! 

7180 Key_new: ! THIS SECTION LOADS BOARDS WITH THE INTENT TO INPUT NEW BOAPD DESC 
7130 Map =3 

7200 DISP "WORKING LOADING BOARDS" 

7210 LOAD "BOARDS", 1 
7220 ! 

7230 Terminate: ! THIS SECTION TERMINATES THE SESSION 

7240 GCLEAR 

7250 DISP "WORKING LOADING STANDARD KEY DEFF IN I T IONS" 

7260 LOAD KEY "STDKEY" 

7279 PRINT LINC20),SPA-: 15), "NORMAL TERMINATION" 

7280 PRINT LIN<2),SPA\15>, " THANK YOU 
7290 DISP 

7300 END ! END OF TERMINATE 

7310 ! 

7320 Pageheadl! THIS ROUTINE PLACES THE PAGE HEADINGS FOR THE INSTRUCTIONS 
7330 PRINT PAGE, TAB < 34-LEN < Si > •' 2 > , "*+ " ; CHPi<. 132) ; Si; CHF:t< 123) ; " *4",LIM<2> 

7340 RETURN 
7350 ! 

7360 Errir.: ! THIS SECTION ALERTS THE USER TO AN ATTEMPT TO INPUT BAD DATA 
7370 BEEP 

73S0 DISP "*>* + 4* INPUT OUT OF RANGE TRY AGAIN" 

7390 WAIT 1500 
7400 BEEP 

7410 RETURN ! END OF ERR I N 

7420 ! 
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7430 Error: ! THIS SECTION IS THE ERROR TRAPPING ROUTINE FOR THE ENTIRE PROGRAM 

7440 ! PROGRAM FLOW RESUMES AT THE TOP OF THE SECTION IN WHICH THE EPPOR 

7450 ! OCCUPED AFTER THE USER PRESSES COHT 

7460 EXIT GRAPHICS 

7470 PRINTER IS 16 

7430 PRINT PAGE 

7490 BEEP 

750O WAIT 30O 

7510 BEEP 

7520 IF EPFN=56 THEN Er rename 

7530 PRINT L I N v 20 > , SPA ' 10), " EPPOR NUMBER " ; ERRN J " HAS OCUPRED IN LINE" ; ERRLJ " . PR 
ESS CONT WHEN READY" 

7540 DISP 
7550 BEEP 
7560 PAUSE 
7570 RETURN 
7530 ! 

7590 Err name.*! THIS SECTION FOR IMPROPER FILE NAME 
7600 PRINTER IS 16 
7610 PRINT PAGE 

7620 Msus*=" DEFAULT MASS STORAGE " 

7630 FOR 1 = 2 TO LEN f N amef > 

7640 IF N&m*«CI, I] = " : '' THEN 7670 
7650 NEXT I 
7660 GOTO 770O 
7670 Msus*=N am*£ C I ] 

7680 CAT Msusl 
7690 GOTO 7710 
7700 CAT 

7710 PRINT LIN<l,',"FiW Name " ; CHP$ < 132 > ; Name f C 1 . I - 1 3 ; CHR5 C 1 23 > ; ” is- NOT on " ; CH 
Rs< 132) ; Msus* ; CHPf<123) : " with That spelling " 

7720 PRINT LIN< 1) , "CHECK OVER THE D I PEC TORY " ABOVE FOR CORRECT NAME OR SPELLING.. 

7730 DISP "PRESS CONT WHEN READY" 

7740 PAUSE 

7750 RETURN ! END OF ERROR 

7760 Plot. elGOSUB Error 

7770 Plot: ! PLOTTING ROUTINE FOP THE AXES 

7780 ON ERROR GOTO Plot* 

7790 PLOTTER IS "GRAPHICS" 

7890 IF Te«i_sol=0 THEN 7330 
7810 Ystep=10 
7820 FOR 1=0 TO 7 

7830 IF Tmax-273>50+I-*25 THEN Ys t * p = 1 5+ 1 *5 
7340 NEXT I 

7850 PI t 1 bl *="Tjunc vs Air Flow" 

7860 Y1 bl $= " June Temp <deg C>" 

7370 GOTO 7940 
7830 Yst € p= . 1 
7890 FOR 1=1 TO 10 

7900 IF Pm i n > . 5* I THEN Y s t € p = . 1 * < I + 1 ) 

7910 NEXT I 

7920 P 1 1 1 b 1 * = '* Prn l n vs Air Flow" 

7930 YlbU="Min Powe r 'Ccmp < Watts)" 

7940 Y m ax = 5 * Y s t * p 

7950 GRAPHICS 

7960 LOCATE 15,120,10.95 

7970 SCALE Xm i n , Xm a , Vm i n , Yrnax 
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7980 AXES Xst ep, Ystep, Xmi n, Ymi n 
7990 CSIZE 3 
3000 LOIR 0 
3010 LORG & 

3020 FOR Xpos=Xmin TO Xm ax STEP Xstep 
3030 MOVE Xpos, Vmin-. l*Vit ep 
3040 LABEL USING "K ,, ;Xpo«*Mwlt 
3050 NEXT Xpos 

3060 MOVE 2.5-Xstep,-.4-Ystep 

8070 LABEL USING " K " ; " A l r-F 1 ow/ Board ( " S.VAL* < 1 •' Mu 1 t > & " M^S'SecV* 

8030 LORG 3 

8090 FOR Ypos=Ymin TO Ymax STEP Ystep 

8100 MOVE Xm i n- . 1 *Xst ep , Ypos 

8110 LABEL USING "k" : Ypos 

8120 NEXT Ypos 

8130 LINE TYPE 1 

3140 LORG 4 

8150 LDIR PI/2 

8160 MOVE -.4+Xstep,2.5*Yst€p 
8170 LABEL USING " K " I Y 1 b 1 $ 

8180 MOVE Xmi n + 2. 5*Xst*p f 5. l^Ysr ep 
8190 CSIZE 4 
8200 LDIR O 

3210 LABEL USING " K " 1 P i c t 1 b 1 f 
8220 LABEL USING " K " ; P 1 1 l b 1 f 

8230 RETURN ! END OF PLOT 

8240 ! 

8250 Bomb:! THIS SECTION DELIVERS MESSAGE TO THE USER OF FAILURE TO CONVERGE 
8260 WAIT 2000 
8270 BEEP 

3230 PRINT PAGE, LI N< 10), "UNABLE TO ACHIEVE CONVERGENCE DUE TO NUMERICAL INSTABIL 
ITIES" 

3290 PRINT LINO). "I SUGGEST A CHANGE IN EITHER THE INSTALLATION PARAMETERS OP” 
83O0 PRINT LIN' 1 • , “THE CIRCUIT BOARD PARAMETERS " 

3310 PRINT LIN<2) . “THERMELEX PREDICTS TEMPERATURES MUCH MUCH BETTER THAN POWER L 
EVELS" 

8320 PRINT LINC 1 ) , "TRY SPECIFYING THE COMPONENT POWER LEVELS . " , LIN < 3 > 

3330 DISP "Press CONT when ready to return to option list" 

8340 BEEP 
8350 PAUSE 
8360 DISP 
8370 Eomb= 1 

8380 RETURN ! END OF BOMB 
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